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1 . 0 SUMMARY 


' / i* * 


Data acquired on the transmission of upstream-generated, low frequency 
noise transmission through aircraft engine-type turbines during NAS3-19435 
showed the existing theory to be inadequate. This program, NAS3-20027, 
was directed towards improvement of the theory and evolution of a working tool 
to predict the low frequency noise transmission through turbines. 

A comprehensive analytical study was performed to define the improved 
theory. Two approaches were utilized in the study: the existing, actuator- 

disk analysis and a new, finite-chord analysis. The frequency dependence 
was preserved through the latter, finite-element treatment of nozzle and 
rotor blades. However, it reproduced the results of the actuator-disk 
analysis for low frequencies and indicated that the simpler actuator-disk 
modeling was valid for frequencies as high as 0.4 to 0.5 of the blade 
passing frequency. This encompasses the entire frequency range of interest 
for combustor noise. 

The existing, actuator-disk analysis and the new, finite-chord analysis 
both utilized an Isolated blade row assumption. The effect of interaction 
with adjacent blade rows (multistaging) was added to the actuator— disk 
analysis, and due consideration was given to spinning modes by modeling 
these as equivalent plane waves. A frequency inverse energy distribution 
corresponding to the asymmetric sound introduction in the NAS3-19435 tests 
was specified for the multiple modes. The resulting analysis was compared 
with the bathtub spectrum specified in the experimental investigation and 
was found to be in very good agreement with the midfrequency floor, that is, 
nominally the 200 to 1200 Hz region. The lobe encountered in the data for 
this region was Indeed shown to be causec by the first spinning mode cut- 
on. Subsequent cut-ons were found to be responsible for much of the data 
scatter noted previously about the nominal floor. In fact, the floor was 
found to extend beyond the frequency range first specified to about 2000- 
2500 Hz. The increase in transmission loss constituting the high frequency 
end of the bathtub spectrum was attributed to the diffraction by the blades. 

The effect of area variations was studied and found to be responsible 
for a spurious increase in the transmission loss data at the very low 
frequencies (below 200 Hz). The apparent increase was caused by location of 
the downstream sensors at a pressure cancellation point in the turbine rig. 

The interaction of the acoustic waves with blade passage shocks was 
analyzed and found to be a very weak, second-order effect. 

The new theory was validated by comparison of the predicted and observed 
trends for the floor as a function of the pressure ratio and speed. 



Two computer programs Incorporating the new theory were written, and 
the program listings are provided in the Appendices, along with user 
instructions. One program is for unchoked turbines and uses an exact 
solutxon method. The other uses an iterative solution and is a generalized 
procedure for any combination of choked and unchoked rows. 

The programs were exercised parametrically and charts constructed to 
approxxmately predict the low frequency noise transfer for single and 
multistage turbines. The transmission loss through the exhaust nozzle was 
found to merit consideration also, and was separately defined. 

Recommendations were made for continuing work and include: 

• Coupling of the turbine and exhaust nozzle wave systems. 

• Completion of the modular prediction method for combustor noise. 
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2.0 INTRODUCTION 
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Studies of advanced aircraft propulsion systems indicate that combustor 
noise is a potential contributor to overall systems noise. This is especially 
true for propulsion systems with reduced fan and jet noise either due to 
cycle selection (for example, high bypass and turboshaft engines), or 
through incorporation of advanced acoustic treatment and/or mixed-flow 
exhaust systems as proposed for the Energy Efficient Engine. There has 
also been much speculation (see Reference 1) that "core,’* "tailpipe," 
or "excess" noise, all of which are generic terms for internally generated 
low frequency noise, constitute a floor in-flight for turbojet engines, 
such as used on the Concorde, or for low bypass engines that might be 
proposed for American AST application. 

Accurate prediction of the different components is an important element 
of systems noise analysis. While General Electric’s Unified Line combustor 
noise prediction method (Reference 2) has been found to be a reasonably ac- 
curate predictor of far-field levels for current engines, there is some 
question about adequacy for engines employing advanced combustors and 
turbines. The Unified Line method consists of a semiempirical correlation 
of engine data and makes no attempt to separate the individual elements. 
Recognizing that the problem is a great deal more complex than a black- 
box approach can cope with. General Electric has been engaged in defining an 
alternative, modular approach to combustor noise prediction under NASA and 
FAA sponsorship. The different modules consist of: 


• Noise generation at the source 

• Transmission through downstream turbine blade rows 

• Transmission through the exhaust nozzle 

• Propagation through the jet stream(s). 

The acoustic characteristics of combustors at the source have been 
researched both experimentally and analytically in recent years (see 
References 1-6). Also, the investigation is continuing most actively at 
the NASA Lewis Research Center and at General Electric under NASA contract 
(NAS3-19736) . The latter involves measurement of the source character- 
istics of an advanced, low emission combustor installed in an engine and 
the associated turbine transmission loss. 

The salient features of low frequency noise transmission through tur- 
bines were determined on a component basis during an earlier NASA contract 
(Reference 7). Comparison of the data with an actuator-disk, isolated- 
blade-row, analytical model (Reference 3) showed the existing theory needed 
improvement. This program contained specific tasks to alleviate the short- 
comings in the existing theory and to formulate an alternative theory free 
of the limitation associated with actuator-disk models. 
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The desired program goals were to: 


1. Define an improved, validated theory for predicting the acoustic 

ransfer function for low frequency noise propagating through 
aircraft engine turbines. 

2. Provide working charts to predict the transfer of low frequency 
noise through single and multistage turbines. 



3.0 THEORY 




3 . 1 BACKGROUND 


An analysis was performed previously (Reference 8) which examined the 
transmission and attenuation of sound waves through a turbine row on the basis 
that both the pitch and chord length of the turbine row were infinitesimally 
small compared to the wavelength of the sound impinging on it. In this 
limit, the turbine row may be modeled as an actuator disk which creates an 
abrupt discontinuity of the flow on either side of it (Figure 1). By employ- 
ing conservation of mass flow and energy flux normal to the blade row, and by 
using the Kutta condition, the attenuation of a sound wave was calculated. 

This analysis was valid only for subsonic flow throughout but was later 
extended to include supersonic relative exit flow under NAS3-18551 (Reference 
3) . One of the key features of the new analysis was replacement of the Kutta 
condition by a choked-flow relationship. The analyses were programmed and 
exercised in a parametric study of the NASA Core, single-stage, high pressure 
turbine. The results are shown in Figure 2 in the form of the predicted atten- 
uation for the plane— wave case as a function of the turbine stage pressure 
ratio with percent design speed as a parameter. The attenuation for the 
nozzle and rotor are shown separately and then summed to provide a stage 
attenuation. The supersonic and subsonic regimes are demarcated, and there 
is little discernible deviation going from one to the other. The predicted 
attenuation apparently increases slightly with pressure ratio over the sub- 
sonic range, remains flat in the transonic regime, and then decreases as the 
Mach number increases to well above unity. 

An obvious problem with this analytical model was the loss of frequency 
content due to the actuator-disk assumption. Also, the upper frequency limit 
on the model was undefined. A second, more subtle problem was the “isolated 
blade row” assumption: that is, the use of anechoic terminations both up- 
stream and downstream of the blade row in question. The effect of adjoining 
blade rows or discontinuities was not addressed. 

An experimental investigation of low frequency noise through aircraft 
engine-type turbines was conducted under NAS3-19435, and the results are 
reported in Reference 7. The data from these scale-model-sized turbines were 
compared with the theory, and discrepancies between theory and data were 
noted. The experimentally determined transmission loss indicated a frequency 
dependence below 100 Hz and above 1500 Hz, increasing in both cases from a 
fairly constant value of attenuation in between. For single-stage turbines, 
the attenuation associated with this “bathtub” floor was found to correspond 
closely to the transmission loss predicted by the actuator-disk analysis. 
However, for a three-stage configuration, the attenuation was overpredicted 
by six to seven dB. An earlier check (Reference 3) of the analysis indicated 
that the attenuation for a six-stage arrangement was overpredicted by 20 dB. 
This clearly indicates that the attenuation for a multistage configuration 
cannot simply be obtained by summing up the attenuation for each individual 
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Figure 1. Geometry of Wave Incident on Stage Element. 
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blade row. The Interactive effects of adjacent blade rows must be given due 
consideration. The interactive effect is integrated into the theory in the 
multistaging" analysis in Section 3.3. 

An analytical model utilizing a finite-chord-airfoil model, in order to pre 
serve the frequency, is described in Section 3.2. The effect of acoustic wave 
interaction with the weak shock waves encountered in the flow passages is 
explored separately. The influence of abrupt area variations is examined in 
an attempt to discern associated frequency dependence, particularly effects 
which would influence the data obtained in NAS3-19435; that is, to note 
trends introduced by the unique facility used to obtain these data. 

These data are compared with the theory in Section 3.4. A computer pro- 
gram incorporating the analysis is presented, along with operating instructions 
and sample printout. A simple, first-cut method of predicting turbine attenu- 
ation for preliminary design use is described. The method is the result of a 
parametric exercise of the analytical prediction program for a number of exist- 
ing aircraft engine turbines. These include turbofans, turbojets, and turbo- 
shafts. The final section consists of conclusions and recommendations for 
future work. 


3.2 FINITE-CHORD ANALYSIS 


The basic idea adopted to consider the effect of finite-chord length 
(and finite, transverse pitch) is illustrated in Figure 3. The process of 
transmission of sound waves across the turbine blade row is "dismantled" into 
an incidence problem, a "passage" problem, and an "emission" problem. In 
other words, as in Reference 9, we assume: the incident sound wave first 

excites duct waveguide modes as if the blade row was a semi-infinite row of 
flat plates; secondly, these duct waveguide modes propagate through the tur- 
bine row as if it were a doubly infinite passage of varying area and a 
straight axis; finally, they reradiate plane waves on the emitted side as if 
the blade row was again a semi— inf inite blade row of flat plates. 

The above idealization considers, to a reasonable extent, the physics of 
the blade row; except, the curvature of the row is not being accounted for in 
the passage problem (though the curvature of the blade row is accounted for 
in treating the two semi-infinite blade rows corresponding to the incidence 
and emission problems as of different stagger angle). With "t" and "M " 
denoting the normal pitch at the inlet and inlet Mach number to the blade 
row, if the frequency of excitation in Hz is below [a A - Mn^/2t] only the 
lowest duct waveguide mode of all the duct waveguide modes excited will be 
propagating. Under these circumstances, Cummins (in Reference 10) has shoxm 
experimentally (with no flow) that even curved bends with 180® turning pro- 
duce very little t ransmiss ion loss. We will restrict the analysis to fre- 
quencies below [a A - Mn‘^/2t] (a is the speed of sound at the inlet to the 
row). The effects of variable area and variable Mach number in the passage 
problem are accounted for. 
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Figure 3. Dismantling of Transmission Process. 



The incidence and emission problems are largely a matter of applying 
the results of Reference 9 and, hence, will not be discussed further here. 
The transmission of the lowest duct wave-guide mode through the variable 
area, variable Mach number, but straight passage region is discussed next. 



The equations governing the propagation of the lowest mode may be writ- 
ten in terms of a nondimens ional acoustic pressure <Kp'/yp) and nondimen- 
sional velocity v(u'/U) as: 


U 



U 


^ _ 

dx 


u) (j) = 0 


(continuity equation) 


(la) 


(-7 - l) f - 3 “)♦ 


+ [2 ^ - j w]v = 0 

(x-momentum equation) (lb) 


In (la) and (lb) , a time dependence for all quantities of type exp (-j oj t) 
is assumed, U(x) denotes the steady, average, axial velocity in the nozzle; 
M(x) the associated steady, average, axial Mach number; and y the specific 
heat ratio . The above equations are given in References 11 to 14. 


p(x) is the average, steady, static-pressure distribution in the nozzle. 


10 



ORIGII'JAL PAGE 13 
POOR QUALITY 


For convenience of the computational scheme to be used, we first intro- 
duce (}) and (<() + v), rather than ({> and v as the dependent variables. We thus 
rewrite (la) and (lb) as: 


U 


dx 


(v + (f)) = j 0 ) 4) 


and 


(2a) 




dx 


= [(y + 1) 


dx 


- 2 j 0) ] (f 


- [2 


dx 


2 w 


(v + 4>) 


(2b) 


Secondly, it will prove useful to choose the independent variable as 
x’ = (L - x) where L is the length from U = 0 to U = a*, where a is the 
sonic velocity at the throat for the equivalent "linear** nozzle (following 
Reference 11) • Equations (2a) and (2b) become: 


U (v + (())= -j 0 ) 


^ {[(T + 1) £r+ 2 j -H 

(1 - M^) 


- |2 £r + j + <^)} 


(3a) 


(3b) 


The **linear** nozzle approximation assumes that U(x) varies linearly from the 
inlet to the outlet, that is: As stated in Reference 13, is a 

**suprisingly satisfactory approximation for conventional nozzles,** We next 
nondimensionalize (3a) and (3b) by using a* as a velocity scale and L as the 
length scale; (3a) and (3b) become: 


(1 - 5) ^ (V + <))) = - j n 4- (^a) 

(1 - 5) ^ = {[-(y + 1) + 2 j n]4> 
dC 

m2 

+ [2 - j n] (v + 4>)} y- (^b) 

(1 - M^) 
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where n = wL/a* and ? = x'/L 

Now, for the linear nozzle, m 2 / Q - m 2 ) may be shown to be 

2(1 - 

(y + 1) C(2 - O 


SO that we have to Integrate the pair: 




< 5 a) 


^ ^ 2(1 - ^) _ . 

de (y + 1 ) 5(2 - 5 ) + 1) + 2j n] 


+ (2 - j n) (<j) + v) } 


(5b) 


J^<"m 1” (with 

0 < Mi < Mf < 1), the Initial and final values of 5 are: 

1/2 


5i = 1 - 


(y + 1) 

L 2 + M 2 _ 1) 


(6a) 


?f = 1 - 


(y + 1) 

2 + M^^(y - 1) 


1/2 


(6b) 


Note that if 0 < Mi < < 1, then 0 < 5i < 5^ < 1. 

Suppose we start the integration near the nozzle throat at 5 = 5^ . As^ 
there is only a transmitted wave in the nozzle; hence, we may show that if 

La ? ♦'^l' Equations ( 5 a) 

and ( 5 b) can be integrated by a Runge Kutta fourth-order scheme from 5 = 5,- 

to 5 - 5 f with the above initial values for <J> and ((|> + v) at 5 = 5 ^ If the 
terminal values of (f. and (^ + v) at 5 = 5f are known, by use of impedance 
relations for forward and reflected waves, (f>inc. at 5 = 5f may be shown to 


1 . 

2 


{ (1 - M^)(() + M^(ij) +v) } 
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where (fi and (<J> + v) are the computed values at C = Cf (where the Mach 
number Is M^) . 

The above describes the essence of the computation scheme that was 
adopted in the present study. Mesh size was normally taken as the smaller of 
( 5 f ~ 5 i )/100 or Tr/ 20 n so that it was the smaller of one-hundredth of the ^ 
(nondimens ional) nozzle length or one-fortieth of a wavelength (based on a , 
the speed of sound at the throat). However, for ^ small or 5 close to unity, 
the derivatives d({i/d^ and d/d^ (v + <j>) can be quite large 5 hence, the mesh 
size was reduced to one-eighth times the lesser of ^ or (1 - O times the 
usual step size for ^ < 0,125 or 5 > 0.875. 

The inputs are M 2 , and a frequency parameter taken here as f = 
[frequency in radians/sec^ x actual nozzle curved length/speed of sound at 
stagnation conditions. Then n may be shown to be: 

n = f /(y + 1) /2 / (Cj ” 

The analysis assumes y = 1.4. It calculates the static pressure ratio 
(Pf/Pi) of the steady, ideal flow. Marble, in Reference 13, shows that as 
p 0 we may expect a result for (p’transm. /p' inc.) of 



1 + M^ 
+ M^ 


(1 ^ ^ m /) 

|(i+^”iM2)| 


(Pf/Pp 


This is the result for "compact" nozzles. As q “ we may expect a limit 
from the point of view that the Blokhintsev energy is conserved at high 
frequencies (as pointed out in References 12 and 14, so that (p'transm./ 
p’inc.) would tend to 




1/2 



(1 + 




1/2 


1/2 


M, 


(1 + 


M^) 


times (p^/p^)- 


The results for M^ = 0.05, 0.1, and Mf = 0.95, 0.975, and for f ranging 
from 0.1 to 20, are shown in Figure 4. Notice the figure shows excellent agree- 
ment at high and low values of "f" with the theories of Blokhintsev and 
Marble. 


The analysis described above for the passage problem was coupled to the 
solutions from Reference 9 for the incidence and emission problems to derive 
the complete, though approximate, solution. 


and M 2 are sometimes used to denote M^ and Mj 


respectively in what follows. 
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To check that such a ’’dismantling” process is valid, comparisons were 
made with the present method and with the actuator-disk method for a very low 
excitation frequency. Excellent agreement was obtained between the results 
of the two methods as shown in Figure 5. 

Repeated calculations with the present method showed, however, that up 
to frequencies (f) defined by 


a X - 

1 


the calculated results are rather insensitive to frequency. Since the rotor 
blade passing frequency can be taken as Wp^/t, where Wp is the wheel tip 
velocity of the rotor and t the transverse pitch of the rotor, the above 
indicates that, up to half the rotor blade passing frequency, the results 
are rather insensitive to frequency. It should be pointed out that more 
exact calculations in Reference 15 for flat-plate cascades bear out these 
conclusions* The passage problem does have a frequency dependence, but it 
turns out that, once the initial Mach number (Mj^) to a row exceeds 0,3, the 
frequency dependence is very slight with even the zero and infinite frequency 
limits being within a dB of each other. 

Thus, the most important conclusion of Section 3.1 was that, in fact, 
the actuator-disk model has a high regime of validity; it is valid up to 
roughly (at least) one-half the blade passing frequency. In practical terms, 
for core noise interests which extend to less than one-half the blade pass- 
ing frequency, there is no need to consider any frequency dependence insofar 
as the analysis of the transmission phenomenon is concerned; although, fre- 
quency dependencies may arise in a given experiment due to the fact that 
given source types couple into a duct in a frequency-dependent manner, and 
incidence angles on the blade row may be frequency dependent. 


3.3 MULTISTAGING 

3.3.1 Problem Formulation 


A sound wave incident on a blade row will generally give rise to a 
reflected sound wave, a transmitted sound wave, and a shear (vorticity) wave. 
The latter two are formed downstream of the blade row and propagate in that 
direction. The former is encountered upstream of the blade row and will 
propagate in a direction opposite to that of the incident wave. 

The transmitted wave will, in turn, be responsible for another set of 
three waves on encountering the next blade row. Further, the reflected wave 
from the second blade row interaction will interact with the first blade row 
giving rise to yet three more waves! It is convenient to collect all the 
upstream and downstream waves after a ’’steady state” has been attained such 
that there exist a pair of forward- and backward-propagating waves between 
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each blade row (see Figure 6). Assuming anechoic terminations upstream and 
downstream of the turbine, the incident wave provides the only forward- 
propagating energy upstream, while a transmitted wave contains all the sound 
energy downstream and propagates away from the turbine. 

In addition to the sound waves, there exist vorticity waves at each 
interface. These propagate with the flow and can only exist on the down- 
stream side of each interaction; that is, the vorticity wave between an 
upstream nozzle and a rotor is determined by the interactions at the nozzle. 

Since the wavelengths of interest here are of the order of a foot, while 
the blade chords and spacings are of the order of an inch, an actuator-disk 
analysis is conveniently applicable. Also, the phase differences between 
interfaces are small and can be neglected, considerably simplifying the 
problem. 

A two-dimensional Cartesian coordinate system, fixed with respect to 
each blade row in turn, is used. Hence all quantities assume their relative 
values at each rotating blade row, as distinct from their absolute values. 

In this analysis, the relative inlet Mach number and the axial component of 
the exhaust Mach number are being limited to subsonic values. At any inter- 
face, upstream quantities will be denoted by the subscript n and downstream 
quantities by m. Hence, in a three-stage turbine, n can assume values from 
one to six, and m from two to seven, as is shown in Figure 6. 


3.3.2 Wave Description 

The wave interaction at each interface can be described schematically as 
in Figure 7. The direction of rotation defines the positive y-axis and 
the axial flow direction the positive x-axls. The flow angles are given by a 
and 3 upstream and downstream -of the blade row respectively. Since alternate 
blade rows rotate and are fixed to each blade row in turn, and are not 

equal but are related by the rotor velocity component. Note that for tur- 
bines 3 will generally be negative downstream of a rotor and positive down- 
stream of a nozzle. 

The sign on the wave propagation angles is defined solely by the y- 
component of the velocity, as the x-components are predetermined by the for- 
ward- and backward-propagation terms. Hence all 0*s shown in Figure 7 are 
positive. 

The frequency across any interface is preserved. However, since the 
acoustic velocity varies and the wave number is defined by m/a, upstream and 
downstream wave numbers are related by 


_ £n 

^n ^m 


(7) 
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Figure 7. Schematic of Wave Interation at a Blade Row. 



where a — ambient acoustic velocity 
k - wave number, w/a 
0 ) - circular frequency, Zirf 

pressure perturbation associated with forward- and backward- 
traveling sound waves can be expressed asr 


PPn " exp j 


L 1 + M cos Gp + M sin Op 
n ny 


- 0)t 


( 8 ) 


and 


PB„ = exp j 


k (-X cos 0B + y sin 0 r ) 
n "n 


1 + M Sin 0g 

nx "n ny 


- cat 


(9) 


w^r the amplitude In thn Incident 

wave, that is, the incident wave is given by: 


p^ = exp j 


kj^ (x cos + y sin 0^.) 

1 + Ml cos 0 + Ml sin 0, 

X I -»-v 1 


- tot 


( 10 ) 


The corresponding density and velocity perturbations are given by: 


pf, 


PPm ~ 


m 


a2 

m 


m = 2, 3, ... 7 


( 11 ) 


The primed quantities denote a perturbation value 
state values. ’ 


ns distinct from steady— 


vp^) = (cos 0p^, sin 0p ) 


Pf, 


m 


“ Pm^m 


( 12 ) 




a^ 


n 


= 1 . 2 , 


(13) 
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'^Bn) = (-cos 6 b^, sin 63 ^) 


Pn^r 




(15) 


(uf. 


Vj) = (cos ©I, 


sin 6j) 


Pi 




(16) 


There are no pressure or density perturbations associated with a vorticity 
wave, hence 


P 





0 


(17) 


The velocity perturbations convect with the flow and assume the form: 

^“Qm’ ^Qy) Qm c^Pj (W X + k^ny y - wt> (18) 


where the direction cosines and Kq^ remain to be defined. 

The y-dependence of all the waves is determined by the incident wave: 


kn sin 0 F 

^_n 

1 + cos 0 p^ + sin 0 p^ 


kn sin 

1 - Mnx cos + Mny sin Gg^ 


(19a) 


km sin Gp^^^ 

1 + cos Gpjjj + Mmy sin Gp^ 


(19b) 


km sin Gg^ 

1 - Mnx cos Gg^ + Mn,y sin Gg^ 
kmy 


(19c) 

(19d) 
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After some manipulation, the following expressions can be derived for 6 r 
. 0Bm* 0Fm terms of the "known" 0Fn = Oj): ' 


tan 0 


Bn 




sin 0p 


n 






tan 0B = _3nn Mmx (1 - Mpy) + G^nTCd ~ Gmn M,nv)^ ~ (1 - M^,,) 

a - Gmn JVy)2 _ g2„ 


tan 0F = ^ fV (1 - Gmn M^y) + ^ [(i _ M„.y)2 _ (l - 

“ ^mn Mmy)^ - 


^SZ _ P 
k ~ ^ 


mn 


where 



sin 01 ? 

^ 

1 ^nx ^^n ^ny ^^n 


The quantity kinx is determined using the fact that the vorticlty wave con 
vects with the flow. That is, the wave will appear fixed (free of time 

dependence) in a coordinate frame moving with the fluid. The coordinate 
transformation is given byi 


^ ^ 


a M t 
m mx 


yp = y “ 


a M 
m my 


t 


The exponent in equation (18) becomes ... 




^m ^^mx ^ny (yp 


^m ^my “ ^t} 
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(20) 


(21) 


( 22 ) 


(23) 


(24) 


(25) 
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Since the time dependence must vanish, 

^nx ^mx ^my ^my ** = 0 


^^m - ^^my ^y 
‘Snx " M„ 


^mx 


since = u)/a, 


m 


Therefore ^ux _ 1 ^^my/^m) Mmy 

^n ^mx 

^nx 1 “ ^mn ^my 
or = ^ 

^mx 




-. T 

f 'j.^ A 


PA.'l-'' 


v*i^ 
U ■ 




(26) 


The direction cosines are determined from the fact that the vorticity wave is 
divergence free, so that 


3x 

This requires 


Qi + ^_2i = 0. 


3y 


k u +k v' = 0. 
mx Qm my Qm 


Equation (18) can then be expressed as 


(Uf 


■fm 




) = 


my 


-k 


mx 


'2 2~ 

+k 

mx my 




2 2- 

+k 
mx my - 


Qm j P^mx (27) 


The reflected and transmitted waves always appear on the opposite side 
of the axis from the incident wave. Using the sign convention of Figure 7, 
this means 


> 0 and 0^ > 0 when 0j > 0 


9r < 0 and 0^ < 0 when 0j < 0 


0R - 0T - 0 when 0j = 0 
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Cutoff Angles 


There are two cutoff criteria for each blade row. 

(a) Upstream Cutoff 

On the upstream side of a blade row, the fact that a wave is forward 
propagating implies that 


l®Fnl ^ 90“ + sin~^ 

This condition can alternately be expressed as 


(28a) 


Un + a^ cos > 0 (28b) 

Hence the upstream cutoff angles are determined by using an equality sign in 
expression (28), Waves exceeding |6 f^| cannot be incident on the blade row 
in question as they convect upstream. 

(b) Downstream Cutoff 

On the downstream side of a blade row, a forward— propagating wave 
implies that 

|0Fml < 90“ + sin“^ (29a) 

This gives cutoff angles of: 

tan 0p cut-off = 

” -Mmx (29b) 


This 3lso defines the transmitted wave angle for which the radical in equa- 
tion (22) becomes zero. For angles larger than this cutoff angle, the 
^®*^Tcal becomes negative and the wave decays exponentially. 

Corresponding to the epm of equation (29) are epn, which can be derived 
using equation (19b) 

or tan 6p„, cut-off - ^n ffl Mnx ( l-GnmMn Y)+GnmV (l-GT,niMnv) ( 1-M^y) 

^ ^~^nm^ny ) ~Cjyj, 

^1 l"*^x‘^os0pjjj+MmySin0Fjjj (31) 

equation (29b). 


where 

and is defined by 


Real values of from equation (30) impose further limits on forward- 
propagating waves that are transmitted through any blade row. 
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3.3.3 Matching Conditions 


Mass and energy conservation provide two sets of equations. A third set 
is derived from imposing the Kutta condition at the trailing edge (this is 
for subsonic relative exit flow; for supersonic flow, the choking condition 
is used instead) . 


Subsonic Relative Exhaust Flow 
The linearized equation for mass conservation gives 

[Up' + pu']^ = [Up' + pu^]^ (32) 

where the subscripts indicate evaluation of the quantities in the square 
bracket on the upstream and downstream sides, respectively, of the actuator 
disk. 


The linearized equation for energy conservation along with the adiabatic 
flow relation, p/pY = constant, in a frame of reference fixed to the blade 
yields: 

[^+Uu' + Vv'] = [2-+Uu" + Vv"] 

P n p (33) 

If a stationary or laboratory coordinate system is used, the rotor energy 
must also be included. 


^3.11y , the Kutta condition requires the flow to leave tangent to a 

• Since the unit vector normal to the exit stream is given by 
(“Sin 3 + cos 3- By) , the Kutta condition gives 


[v^ 


(- sin 3 


e + cos 3 

X 


or 


e )] = 0 

y m 


[-u'^ sin 3 + v^ cos 3]^ = 0 (34) 

In general, the quantities both upstream and downstream will consist of 
a forward— propagating sound wave, a backward— propagating sound wave, and a 
vorticity wave. However, upstream of the first blade row there is no vor— 
^Bve (Qj = 0) , and downstream of the last blade row there is no 
backward-traveling sound wave (B2 n+ 1 = 0)* where N is the number of stages in 
the turbine. Since = 1, that leaves 6N unknowns. However, there are 2N 
blade rows with three equations at each blade row. Therefore the problem can 
be solved. 


Application of the matching conditions (32) - (34) to the first blade 
row gives the following equation set which can be expressed in matrix form 
as: 
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cos 



where 



Using 
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Then 



where D 


-1 



Is the inverse of D^, 


that is in 



D, 


gives the 


( 38 ) 


identity matrix 



10 0 
0 10 
0 0 1 


Similarly, for any blade row it can be written: 



( 39 ) 


( 40 ) 
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For the last (2N) blade row, = T and = 0. therefore 


2Nfl 


2N+1 


^2N^ 


2N 

^2N 

^2N 


^^2N ^2N^ ^°2N-1^2N-1^ 




1 

Bl 

0 


TC^l 

TC^2 

TCi3 


’ 1 j 

TC21 

TC 22 

TC23 


«1 

TC31 

TC32 

TC33 


0 1 


transition coefficients relating the transmitted and 
incident perturbations. 

The second row of ( 45 b) shows that 
TC. 


~ TC 


21 


22 


whereupon, it can be seen that 


T = ITC 


(’ 


ICjj 

^11 ■ TC 22 


) 


A computer program to utilize this matrix- inversion technique can be 
found in Appendix A, along with a flow chart and typical output. 


(45 a ) 


(45b) 


(46) 


(47) 
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Supersonic Relative Exhaust Flow 


When the relative flow exiting from a bj.ade row becomes supersonic, the 
Kutta condition is replaced by the choked— flow condition. A discussion of 
application to disturbed flow at a blade row can be found in Reference 16. 

The interaction with the shock that occurs due to the locally supersonic 
conditions is considered separately in Section 3.4. 

Supersonic flow actually implies two separate governing equations — 
one upstream of the blade row and the other downstream. The downstream 
condition is analogous to the Kutta condition in that it determines the 
relative exit angle. The Kutta condition states that the relative flow angle 
leaving the blade row is given by 

3 = cos ^ (d /t) = constant 
o 

where d^ defines the cascade throat and t the blade-to-blade pitch (see 
Figure 8), However, when the critical pressure ratio is exceeded, the flow 
angle for low supersonic Mach numbers is given by: 


(48) 



where A/A* is defined as in the usual sense (Reference 17): 


A 

A* 


i] f 2 + (y-I)M^ 
MJ L Y+1 


exp 


Y+1 

2(y-1) 


-] 


(49) 


The one-dimensional area function defined in (49) is valid only for 
small supersonic Mach numbers because it ignores shocks. The flow turning 
provides the extra area required to pass the flow defined in the throat. 


However, the downstream choking condition and the mass conservation 
equation cannot both be used simultaneously as the former implicitly contains 
the latter and the resulting equations are no longer linearly independent. 


The upstream choking condition requires that the corrected mass flow be 
dependent only on the upstream stagnation parameters (Reference 16). That 
is : 



constant 


or 



constant 


(50) 
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where: m = mass flow rate = pUA 

A = cross-sectional area 
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= stagnation pressure 
Tq = stagnation temperature 
R = gas constant 
Y “ ratio of specific heats 


Using conventional gas dynamic relationships and taking the logarithmic 
differential yields: 


2 h' _ I±1 a' = 0 (51) 

U T Y-1 1 


where: T' 

y 


= temperature perturbation 

■ ^ 


M = absolute flow Mach number 
abs 

After some further simplification and assumption of isentropic flow (see 
Reference 3), the following equation in u’, v' and p’ results: 


^ “abs ♦] 

- [(Y+1) - 0 (52) 

where (|) = absolute flow angle* 

Proceeding as in the subsonic flow case, with equation (52) replacing 
equation (34), the (An) and (Dn) matrices assume the following form: 


= 


(M H-cos0p ) 
nx ^n 


(M -COS0B ) 
nx "n 


k /k 
K 
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. . M (k /k )-M (k /k ) 
[1+M cos(a -ep )] [1-M cos(a +8 b )] nx ny n _ ny — nx — n_ 
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where 


nx n 

*32 • - M^‘=»s%+(y+ 1)M„ cos(»-+eB ) 

o.. (k /k ) 

A == AiL yn n 
33 M (Y+l)Mn 

nx n 
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CT32 

CT33 



(56a> 


(56b) 


Here (CT) is the transition coefficient roatr^fv f nr* 

flow turbine. Equation (56b) can be used to obtain ®^haust 


CT, 


T = F„ = 33 

2N+1 CT„CT_-CT_CT, 
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Unfortunately, most turbine configurations incorporating supersonic 
exhaust flows do so only for the Initial few blade rows. The matrices decouple 
at each subsonic/supersonic interface, and neither (TC) nor (CT) can be 
defined. 

There are several alternative solution methods, including following 
acoustic waves through successive interactions with blade rows. This ap- 
proach is used later on to validate the matrix- Inversion technique for single- 
stage turbines. The implementation, however, becomes quite cumbersome and 
complex for multistage turbines* 

A generalized solution results from the realization of the fact that, 
out of the six amplitudes involved at each blade, two are fully defined at 
the first blade row (F^ = 1 , = 0 ). Guessing at one of the other four 

amplitudes, the other three unknowns can be obtained utilizing the three 
matching condition relationships at the first blade row. Since F2, B^.and 
are now known, Fo, B3, and Q3 can be found by using the relationships at 
the second blade row. Finally, F2N+1 and Q2N+I are calculated. Since an 
anechoic termination is assumed, B2N+1 = 0 . If the computed value of B2N+I 
is not zero, a second Iteration is made through the turbine. Note that this 
guessing routine allows for solutions of nonanechoic terminations. It is 
sufficient to define the relationship between F2N+I and B2N+I to the 

termination. Then the computation loop-escape condition becomes (B2N+I/ 

F2JJ4-1) convergence to the ratio determined by the termination rather than 

B 2 N+I = 0- 

Implementation of this solution routine is made somewhat complex by 
supersonic exhaust blade rows because only two equations are available to 
define downstream quantities. Therefore it is necessary to start a new guess 
at each supersonic blade row. The solution scheme is outlined in Appendix B, 
along with a time-share program listing and typical output. An interesting 
result of supersonic-flow blade rows is that the sound Waves move upstream 
slower than the flow moves downstream, and therefore negative values of the 
backward-traveling wave become possible. 


Validation of Multistaging Approach 

An acoustic wave incident on a multiblade- row vehicle will generally 
give rise to a system of acoustic and vorticity waves which can be evaluated 
in two different ways. The current, multistaging analysis postulates an 
"equilibrium" state solution; wherein, all the reflected and transmitted 
acoustic waves are combined into a pair of forward- and backward-travelling 
acoustic waves in each Interb lade-row space and the associated vorticity into 
a vorticity wave downstream of each blade row. The other approach considers 
each blade row interaction as an isolated blade-row Impingement and then 
follows the resulting reflected, transmitted, and vorticity waves through 
successive interactions with adjoining blade rows. The solution the lim 
of infinite interactions should approach that of the equilibrium model. This 
has been verified for a number of cases ranging from low to high pressure 
ratios, zero and nonzero acoustic wave incidence. Two representative compar- 
isons are provided in Table I for the first stage only of the HLFT IVA low 


35 



Comparison of Successive Interaction 











pressure turbine tested in NAS3-19435. Case I corresponds to the lowest 
pressure ratio tested for the full three-stage turbine, while Case II corre- 
sponds to the highest pressure ratio, at which the first stage was very 
nearly choked at 100% speed. The first column provides the acoustic wave 
amplitudes for the isolated blade-row interaction in which only the trans- 
mitted wave at each blade row is preserved; the reflected and vorticity waves 
are discarded immediately after the interaction. The second column gives the 
amplitudes if the vorticity wave from the first blade-row interaction were 
also preserved and made to interact with the next blade row. The succeeding 
columns contain the amplitudes due to successive interactions of the re- 
flected wave from the second blade row. The last column gives the values 
predicted by the raultistaging computer program. The convergence of the 
successive interaction solution to the multistage values is surprisingly 
rapid, particularly for low pressure ratios. For example, the final trans- 
mitted wave amplitude reaches a value of 99% of the multistage solution after 
only two interactions at the lower pressure ratio. At the higher pressure 
ratio, the transmitted wave amplitude reaches 92% of the multistage value 
after two Interactions, and 96% after four. 


3.3.4 Energy Transmission 


The energy transmitted can be computed using the results of Blokhintsev 
(Reference 18). The energy density e is given by 


• e 



(a+V 


abs 



(58) 


where ^abs absolute flow velocity and 6p the unit vector normal to the 

wave front. Also, the intensity flux vector t is given by 


I = 


e (ae +V 
P 


abs"^ 


(59) 


Only the axial component is of interest here 



I = 

e(ai +V ^ )-e 


X 

p abs X 

or 

I = 

G(a cos9„4U) 


X 

F 

also 

->> 

^abs 

= Ue +(V+W_)e 
X ^ y 

where 

Wr is 

the rotor wheel 


37 



therefore 



where Wj^/a, the blade tip Mach number, 
conditions at exit from the last blade row, 
inlet to the first blade row, i,e., 1=1, 


and the subscript T would denote 
i.e., T = (2N+1) , and I those at 


For a first approximation, 
plane waves propagating between 
Morfey (Reference 19). 


annular spinning modes can be treated as 
infinite p±ates ~ as was demonstrated by 
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A plane wave approximation for m = 3 spinning lobe is provided as an 
example. The annulus is assumed to be cut and straightened out (unwrapped), 
so that the cylindrical walls become a plane sheet. Continuity in the 
circumferential (y) direction requires that the wave pattern be repeated 
every 27 Tr* (mean circumference); that is, 


or 


M( A - r ~) = 27Tr* 
sin 0^ 


. mX 

sin 0^ = 


m 


27rr’^ kr' 


or 


. “1 , m V 


where m = circumferential lobe number 


= 0 , 1 , 2 , . , . y 

k = wave number 




r* = root mean square radius 


_ (tip radius)^ + (hub radius)^ [ 

’ L 2 J 


1/2 


(62) 


l-/hen more than a single dimension is involved, the wave number k is the root 
mean square of the wave numbers associated with each of the dimensions, e.g. 
in the axial and circumferential directions. 


Note that m = o corresponds to a plane wave propagating axially down the 
annulus and is the only cut-on mode for kr* <1. As soon as kr* exceeds one, 
the first pair of spinning modes (one corotating and one counterrotating) 
appear - as was indeed observed during the NAS3-19435 tests. 

Each mode is associated with a different incidence angle, and the corre- 
sponding transmission loss can easily be computed using (61) . The question 
now arises as to the appropriate energy assignment. Equal energy distri- 
bution between all cut-on modes has been frequently postulated in fan noise 
and treatment work. Experimental observations indicate that this is not an 
unreasonable distribution for symmetric sources particularly. The siren tone 
injection into the turbine plenum during the NAS3-19435 tests corresponded 
closely to a point-source placed in an annulus. A simple, no-flow, analyt- 
ical modeling (Appendix C) of the resulting duct coupling can be used to show 
that the energy distribution is given by 


E 

m 


(f2 - £^2, 


1/2 


(63) 


39 



where 


Em ~ energy assignment to m*^^ mode 

f = frequency of interest 

fj. — cut— on frequency for inverse 


An obvious outcome of this frequency inverse dependence is that all the 
available energy is biased towards a mode just cutting-on. But 6t for this 

cut-on, almost ensuring complete reflection at 
the blade row Hence, cut-on should be associated with a sudden increase in 

inconsistent with observations made during 
NAS3-19435, as will be shown in Section 4 . ® 


Once the energy assignment has been made, it is a simple matter to 
compute the summed transmission loss for any given frequency. The computer 
programs in Appendices A and B provide both the individual transmission 
osses for each cut-on mode and the summed transmission loss. 


3.4 SECONDARY EFFECTS 


3.4.1 Duct Termination and Area Changes 

"^e area variations encountered in the turbine tests of NAS3-19435 may 
be modeled as shown in Figure 9. There is a gradual area change from the 
inlet plenum to the inlet casing (Si to S 2 )‘, there are sharp area changes 
associated with each blade row (S 3 and S 5 ), and then there is a sudden expan- 
sion as the exhaust flow dumps into the exhaust plenum (S^ to S 7 ). Each area 
discontinuity is associated with reflected and transmitted waves. The answer 
being sought here is the effect on the transmission loss and, in particular, 
he unique or spurious effects imposed on the data acquired during NAS3-19435. 

The area changes associated with the blade rows are properly accounted 
or in the analyses, but not the associated phase changes over the lengths 
M* 1 : 5 , i 6 » etc. The multistaging analysis, for example, assumes negligible 
change in phase over the interblade- row spacing Since /c = 1.31 cm for 

the turbine of Figure 9, the actual phase change [angle in degrees = (snacinc/ 
wavelength) x 360 ] would be about 1° at 100 Hz and 18° at 2000 Hz which 
represent the J-imits of the frequency range of major interest. Hence the 
assumption would be strictly valid only at the low frequency end. 


The major Impact would appear to be that of the area change at the 
exhaust Plenim. It will be shown that the reflected wave at this termination 
is almost 180 out of phase with the incident wave for low frequencies and 
has an amplitude almost as large, making the duct termination a pressure 
node. Pressure measurements in this region would then indicate inflated 
values for the transmission loss. The exact degree of pressure cancellation 
at a given sensor is a function of the amplitude and phase of the reflected 
wave, the wave number, and the distance to the sensor (£9 or £c + £q). To 
our knowledge, there are no exact solutions available in literature applic- 
able to this particular geometry. However, several approximate methods are 
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f the strip theory modeling by Mani (Reference 20) which 

"he somewhat simpler no-flow models used to analog 

ReSrence®2^ ^ (see. for example. ® 

strator? ‘ analysis is perfectly adequate here - as a demon- 

Assuming, for the moment, a cylindrical duct of radius r discharging 
into the plenum, the ratio of the reflected to incident wave can be written: 


!e _ ^^o - j ^o 

■6 (^o ^ Po^/^6> + j ^o 


lnterfa° reactive components of the Impedance at the 

interface and S is the cross-sectional area. 

In the l^it that (Sy/Sg) is finite, and the wavelength is large com- 
pared to the duct characteristic dimension, R^ = (p^a/S^) and == 0. Then, 

i ^7 - ^6 

rL'lJIts’^Lrobtain^":'' """ """ = 5160 cm2, the following 


High Pressure Turbine 


= - 0.80, ATL = 14 dB 


^c-Stage Low Pressure Turbine 


= - 0.62, ATL = 8.4 dB 


Three-Stage Low Pressure Turbine 


- 0.37, ATL = 4 dB 


SLelL^L^ at thf 1 transmission loss due to pressure 

5160 cm2 t ‘iownf ream sensors. Note that Sy actually varied from 

exhaust duct termination to about 13700 cm2 the scroll 

:?ssLTLse:?^^ th:°trans- 
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Table II. Exhaust Duct Termination Effects. 


Exhaust Duct Duct 
Area, S 5 Height 

Turbine Configuration (cm2) (cm) 


High Pressure 562.5 
1 Stage, Low Pressure 1206.5 
3 Stage, Low Pressure 2387.2 


Note: S 7 - 5160 cm^ at the termination, but increased to about 13700 

cm^ at the exhaust scroll collector. 


3.81 

6.60 

12.45 


Lengths, J., for Sensors (cm) 
Wall 
K, 



Probe 


^10 

4.06 

6.60 

5.84 

3.30 

7.37 

4.83 


For the case of very large (Sy/Sg), the impedance can be considered the 
same as that acting upon a piston mounted in an infinite baffle: 

pa 

R = -^ R (2 kr) 
o S. 

6 ( 66 ) 

P ’a 

X = X (2 kr) 

o S ^ 

D 

where 2 4 ^6 

° (2H4) " (2) (42) ( 6 ) ( 2 ) (42) (62) ( 8 ) 

4.x x^ , X^ . 

^ ” (3^) (5) (32) (52) (7) * • * ^ 


^ = R(2kr) - 1 + i X(2kr) l 

R(2kr) + 1 + j X(2kr) 

For example, using the truncated series representation of Equation 67, 

B (o) 

kr = 0.2 gives / y = 0.99 exp [j (170“)] 

F6(o) 

where (o) mean kx = 0.0 


* See Reference 21, page 146 
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r-cTi the reflected and incident waves provide almost complete can- 
cellatxon at the duct termination. At higher frequencies. ^L cance^Lt Jon 
is not as complete because of changes in both amplitude aAd phaseJ 


kr = 2 gives 


Bg(o) 

F6<°) 


= 0.554 exp [j (107”)] 


The effect at the measuring station can be computed using: 

B6<i.) = Bg(o) exp [j(-k£)] 

(69) 

^6^^^ “ ^6^°^ fj(kJl)] 

7)“nd‘a/+ I MS3-19435 Final Report, Reference 

tLblne rite ’ i: J; *>■' "^8^ pressure 

curoxne tests. It xs obvious that these measurements were very nearlv in 

pressure cancellation region. In contracr rh^ i ^ nearly in the 

1 J in contrast, the low pressure turbine t-ranc- 

?orwMc;T„fr:„:rLr1er'"l2'!9®«^ »au-»nnted Knlltea (K5 and K6) 

the wall and probe sei::L „e givL”n l;ble“i. ‘ 

in apparent tranamlLLriLf SL^Jere^obse^vS^te^ dpurloua Increases 

Tz tLti:rdisJort‘L‘ ..rZizx 

becan^r^t:^^ 

s;si»'L'''”arit“f spctioid^ ?L"c““ 

generated bv tL^Ln ^‘'’'ucture a test to measure wave patterns 

g me ry, rather than measuring real transmission characteristics 

zero net ’effect/^ « Rood first estimate would be to assume J 

plane’iLertef "®or"^S;^l“tr ''hl“ »»ly for no-flow and 

Ss^lo^r"' fn.ifcate“ffL““ll®''Lc:e2: ”"hrtr"L-°' 

conflg^ratloisr'tS arL“IrlatirJn?re “““■foolo" In engine 

would be that at the core nozzle exit THp ^’^®’^®®^sslon loss 

for nozzles such that (2kr) >> 1 Evpti t-h would be a nonzero B 2 N+I 

turbine transmission iL's^^ill J;sufr 

ly large incrp^cp ^ I nuwever, mere will be a comparative- 

overlooked. transmission loss which should not be 
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/j, 


There is also potential for a shift in the transmission loss spectrum 
due to the "gooseneck" sometimes encountered between the high pressure turbine 
exit and low pressure turbine inlet for high bypass turbofans. The gooseneck 
is typical of the CF6 family of engines and Involves a large Increase in the 
mean radius. The modal content of the acoustic energy propagating between 
the two turbines will change, since the first spinning modes will cut-on at a 
lower frequency (cut— on is computed using kr* =1, 2, .... along with a 
Doppler correction for flow). That is, the sudden increase in transmission 
loss characterizing modal cut-on could shift to lower frequencies. 


3.4.2 Shock Interaction 


Since turbine blade passages are not normally designed as converging- 
diverging nozzles, the existence of supersonic flow results in shocks in the 
vicinity of the blade passage— but only at the trailing edge, as illustrated 
in Figure 10(a) (See Figure 21(c) of NACA RM EIK25 for Schlieren photograph 
of such shocks) . 

The interaction of acoustic waves with shocks has been investigated 
analytically by Landau and Lifshitz (Reference 22) for normal shocks and by 
Moore (Reference 23) for oblique shocks. 

In general, the weak disturbance field resulting from shock interaction 
with an acoustic wave can be considered to include two components: 

(a) an unsteady, isentropic, irrotational perturbation satisfying the 
wave equation, i.e., a sound wave 

and (b) a steady (relative to the flow), rotational perturbation of con- 
stant pressure, i.e,, a vorticity wave. 

Strictly speaking, an entropy wave is also created (Reference 16). However, 
the acoustic perturbations are assumed to be small and the shock weak (the 
flow in turbine passages will rarely exceed M = 1.2). Under these circum- 
stances, it would appear that the resulting entropy waves could be neglected. 

As shows in Figure 10(b), Moore discusses the case of a shock overtaking 
a sound wave (Problem A) , and that of a sound wave overtaking a shock (Problem 
B) . The case of interest here corresponds to Problem A in his frame of 
reference. Within the blade passages only the zeroth order mode, an axially 
propagating wave, can be cut-on for the frequency range associated with 
combustor noise. Referring to Figures 10(a) and 10(c), one can see that the 
incidence angle i|; between shock and acoustic wave can then be taken as 
approximately zero. The case of interest here then corresponds to Problem A 
in Moore’s frame of reference with M 1 and = 0. Using the appropriate 
results, the net effect is a weak refraction of the incident sound wave, as 
shown in Figure 11(a) and (b) , The associated vorticity wave occurs at 
4*3 ''' 4>l/2 (approximately parallel to the shock in this case) (Figure 11(c)), 
but the velocity and density effects are very nearly zero, even for Mach 
numbers up to 1.5 (Figure 11(d)). The order of magnitude of the overall 
effect would appear to be much smaller than that resulting from the actuator- 
disk interaction and may be ignored for all practical purposes. 
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• ^ i Trailing Edge 



(a) Shock Pattorns in a Transonic Turbine Cascade 


• Coordinuto Frame Moving with Fluid. 




Sound Wave Overtaking Shock from Behind (Problem 


(b) Shock Interaction with Sound Waves (Reference 6) 


■ Coordinate Frame Moving with 
Fluid 



(c) Shock and Acoustic Interaction 
in a Moving Coordinate Frame 


Figure 10. Shock and Acoustic Interaction. 
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Figure 11. Shock Interaction with Sound Waves (Reference 6). 
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^ • 0 theory/data comparison 


4.1 


BACKGROUND /DATA ACQUISITION 


An experimental investigation of the low frequency noise transfer 
through aircraft engine-type gas turbines was conducted at General Electric 
un er NASA Lewis Research Center sponsorship (NAS3-19435) Details of 
test end the results obtained can be found L Reference ! ?^ese Ltf a^e 

l^tTon It of the previous 

sr»d ^L-fficiftLri- 


The program objectives in NAS3-19435 wer^ ^n n ^ ^i. 

transtasslon loss of sound injected ip”Lr:f he urblnns^funS^n" 

the^r: “O'* oo»pare thLi da^rnith erlstjnr 

theory in order to assess the validity of the theorv Th^ ^ f ^ . 

order to accomplish these objectives is outlined in Figure 12." IWo^turbines 
single-stage, high-pressure turbine (NASA core) and a 

^SrhSs“f^; r°:;?d:rirT:b!L^^n a^: 

acquired at both choked “nd unchoLd'‘coStilnsr”^^®“"''‘°”‘ 

FacllUy (PlgS “^rsolrsourcrco^K^f of f Mgf ■ 

coupled to the inlet plenum through a tranSti^^^oln and^a rl^S 

tra™uiar?i»°p°rth'®"’* acco.pUsM^r^ch;ngf^°“cj!““S 

tions also included one parr"r“S-Lun^^d K^iL'r “hf^""' 
scheLticalirirF^JuJes if and sections are displayed 

ton! „e1e\\s- r^^a^^ijj Jb^f^ L“1dj!S—t!: 
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High Pressure Low Pressure 

Turbine Turbine 

NASA Core | j HLFT-IVA 

1-Stage 1“ and 3-Stage 
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Figure 12, Structure of Acoustic Tests to Accomplish the Program Objectives. 





Table III. High Pressure Turbine Design Characteristics 
(NASA Core Turbine). 


w/r 

Wt. Flow Function, ’-p 

0.81 

Loading, 

£ Ij2 
P 

1.66 

Pressure Ratio (Total) 

1.83 

Speed, N//f 

362 

Stator Vanes 

36 

Rotor Blades 

64 

Radius Ratio 

0.85 

Tip Diameter (Stage Exit), (cm) 

50.8 


Table IV. Low Pressure Turbine Design Characteristics 
(Highly Loaded Fan Turbine, HLFT-IVA). 



— - Stage 



1 

2 

3 

Overall 

Wt. Flow Function, 





’ P 




1.57 

Loading, 

ZU 

3.52 

3.12 

1.60 

2. 70 

P 





Pressure Ratio (Total) 

1.73 

1.81 

1.41 

4.72 

Speed, 









204 

Stator Vanes 

100 

144 

140 


Rotor Blades 

206 

190 

160 

_ 

Radius Ratio 

0.811 

0.735 

0.663 

_ 

Tip Diameter (Stage Exit) (cm) 

63.55 

69.08 

73.18 

- 
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(b) HLFT-IVA Low Pressure Turbine, 1-Stage Build. 


Figure 14. Schematic of Low Pressure Turbine Configurations. 
52 
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Figure 15. Schematic of NASA Core High Pressure Turbine Vehicle. 



















Table V, High Pressure Turbine Test Matrix 
(NASA Core Turbine) 
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Table VI. Low Pressure Turbine Test Matrix 
(HLFT-IVA). 


• Desigr Speed, N/ = 204 

• Flow Function, w/t/P = 1.57 

• Inlet Total Pressure PxQ ~ 275.8 kN/m^ 

• Inlet Total Temperature, Txo = ^22 K 


Pressure Ratio (Pxo/Ps2) 


% Design 
Speed 

Speed 

(rpm) 

Single-Stage Build 

Three-Stage Build 

1.6 

1.9 

2.2 

2.5 

2.0 

3.0 

4.0 

5.2 

50 

2100 

X 

- 

- 

- 

X 


- 

- 

70 

2940 

X 

X 

X 

- 

X 

X 

X 

- 

90 

• 3780 

X 

X 

© 

■ X 

X 

© 

X 

X 

100 

1 

4200 

X 


© 


X 

X 

X 

© 

110 

4615 

X 

X 

X 

X 

X 

© 

X 

X 


O - Repeat Point 
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Data analysis techniques included very high resolution data reduction 
and coherence analysis between upstream and downstream sensors in an effort 
to unmask the siren tones downstream of the turbine. The latter was found 
to be more successful. A typical coherent spectral comparison is shown in 
Figure 16. The figure clearly shows large transmission losses for the 400 
and 800-Hz tones, but a much smaller value for the 1200-Hz tone. The com- 
parison of Figure 16 is on a SPL basis. A more meaningful result was 
obtained by correcting the data frr flow, specific impedance, and area to 
arrive at the corresponding power levels (see Reference 7). The area 
correction assumed zero-th order, radial-mode distribution, that is, constant 
energy distribution from hub to tip. This has been found to be a reasonable 
assumption for low frequency noise measured in an engine core (Reference 2). 

Typical plots of the siren tone attenuation as a function of the tone 
frequency are shown in Figure 17 for the high pressure turbine at design 
point. The spectra display a very distinct, bilobed shape, with large 
increases in attenuation below 100 Hz and above 2000 Hz and a secondary 
peak between 350 and 400 Hz. This secondary peak was found to correspond 
to cut-on of the first spinning mode. How this cut-on increases the trans- 
mission loss has already been discussed in Section 3, 


The data appeared to exhibit a fairly large amount of random scatter 
possibly as a consequence of duct-related phenomena and Interference between 
forward- and backward-propagating acoustic waves. The 2.54-cm axial spacing 
between sensors was found to be Inadequate to separate the two wave systems 
because of the large wavelengths and high broadband "noise” levels. Ultimately, 
the only viable option available was data averaging - use of large samples 
and as many of the sensors as possible. The midlobe, however, remained 
readily discernible, even for the low pressure turbine data where the siren 
frequency corresponding exactly to the first modal cut-on was assiduously 
avoided. Partly because the size of the midfrequency lobe was believed to 
be a consequence of the asymmetric sound injection into the turbines, and 
facilitate comparison with the existing theory at that time 
Uy76), a bathtub spectrum shape was postulated as shown in Figure 18 The 
attenuation spectrum was divided into three distinct regions as shown ’in 
^ frequencies (below 100 Hz), midfrequency floor (200 

to 1200 Hz), and high frequencies (above 1500 Hz). It was hypothesized 
mechanisms Involved at the low and high frequency ends 
wh ch either invalidated the basic theoretical (actuator disk) assumptions 
or involved phenomenon not considered in the analysis. The bathtub floor 
was found to correspond closely to the actuator-disk theory. Coincidently , 
the floor spans the major frequencies of interest for combustor noise. The 
frequency span of the floor could easily be extended to 2000 Hz without any 
loss of generality, as is obvious from Figures 17 and 18. 

transmission loss corresponding to the floor was 
obtained for each turbine operating point by averaging the attenuation 

nrivfn f I r midfrequency region. This technique 

proved to be remarkably successful in collapsing the data and revealing 

trends. The collapse achieved is shown in Figures 19 through 21 for the 
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Sound Pressure Level 


• High Pressure Turbine 

• Siren = 1204 rpm 

• Inlet Temperature = 450 K 

• Inlet Pressure = 389.6 kN/m^ 



0 400 800 1200 1600 2000 

Frequency, Hz 


Figure 16. Comparison of Upstream and Downstream Signals 
Showing Turbine Transmission Loss. 
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• NASA Core High Pressure Turbine 

• Averaged Values from Coherent Spectra 

• Inlet Temperature = 450 K 


• " 2.14, N//T = 100% 



Frequency, Hz 

(a) HPT Cold Design-Point Attenuation Spectra 


• NASA Core High Pressure Turbine 

• Averaged Values from Coherent Spectra 

• Inlet Temperature = 783 K 



Frequency, Hz 

( b) HPT Hot Design— Point Attenuation Spectra 

Figure 17. High Pressure Turbine Design-Point 
Attenuation Spectra, 



Attenuation for Design Point 



Figure 18. Bathtub Spectrum Shape 











PWL 


□ 110% N//T 

O 100% N//t • 1-Stage HLFT-IVA 

A 90% N//T • P,j.^ = 275.8 kN/m^ 



Figure 19, 


Effect of Turbine Pressure Ratio on Attenuation of Single- 
Stage Low Pressure Turbine, 



110% N//f 



Figure 20. Effect of Turbine Pressure Ratio on Attenuation of Three-Stage Low Pressure 
Turbine. 







□ 110% N//f 

• 100% n//t 
A 90% N/v/t 
Q 70% N//T 


• NASA Core 

• = 389.5 kN/m 







turbines tested. The plots* show the floor transmission loss as a function 
of the turbine pressure ratio, with the turbine speed as a separate parameter. 
Clearly, the turbine speed is not a significant variable. The data trends 
did indicate a pressure dependency in that the attenuation increased (very 
slightly) with the pressure ratio for subsonic flows, flattening out, and 
even decreased by a small amount for choked flows. However, the total 
variation observed for any turbine was about 3 dB or less over the entire 
test matrix. The test matrices for these component tests represented far 
greater excursions from design than would be encountered for turbines 
installed in engines. Hence, the data trends would suggest very minor, 
certainly less than 3 dB, changes in midfrequency transmission loss over 
the normal operating range for aircraft engine turbines. 

These figures also show the analytical predictions using the actuator- 
disk theory of Reference 3. The prediction involved two major assumptions 
in addition to the actuator-disk modeling. First, only the plane wave 
propagating axially down the duct was considered (0i = 0); spinning modes 
were ignored because of the low frequency nature of the sound. Secondly, 
it was assumed that the attenuation due to each blade row could be computed 
separately with anechoic terminations both upstream and downstream and the 
individual attenuations were additive in arriving at the attenuation for 
the turbine. Both assumptions were necessary in order to maintain a viable 
mathematical model and extract a solution. Comparison of the predictions 
with the data trends in Figures i9 and 20 left little doubt that the existing 
analysis needed further modification. Figure 19 shows remarkable agreement 
for pressure ratios below choking, but the pronounced dip in predicted 
attenuation above choking was not matched by the data trend, and a 3-dB 
discrepancy resulted. Further, while good agreement was ^ound for single- 
stage turbines in the subsonic flow regime, the three-stage turbine data 
were overpredicted by 3.5 to 7 dB proceeding from the lowest to highest 
turbine pressure ratio tested. The question then became: could the actuator- 

disk theory be modified sufficiently through recognition of higher order 
(spinning) modes, multistaging, etc., to obviate the discrepancies noted 
above and explain the observed frequency spectrum? 


4.2 COMPARISON OF THE DATA WITH THE IMPROVED THEORY 


The predictions used here were generated using the computer programs 
listed in Appendices A and B. The program in Appendix A can be exercised 
only for subsonic turbines, while that in Appendix B is a generalized 
program which can accommodate both choked and unchoked blade rows. However, 
the Appendix A program provides an exact solution and is considerably 
cheaper to execute. 


* 

The values are slightly different from those shown in Reference 7. The ASPL 
to APWL conversion in Reference 7 was made using average values (one for 
each turbine) of the specific impedance and Mach number in order to facili- 
tate data reduction. The data shown here have been corrected using the 
exact values for each different operating point. 
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Figures 22 and 23 provide data comparisons with predictions using the 
program of Appendix A for the single-stage, high pressure 
(NASA Core) turbine. The prediction in Figure 22 uses an equal energy 
distribution and is seen to skim along the bottom of the data points. 

There is an increase in attenuation at each modal cut-on frequency. The 
effect of the first one is most pronounced; suddenly two- thirds of the 
incident energy is transferred into the two new waves that are completely 
reflected. At the second modal cut-on, two-fifths of the incident energy 
IS transferred into the new waves; therefore, the indicated Increase in 
transmission loss is correspondingly smaller. As the number of existing 
modes increases, the effect of subsequent cut-on naturally diminishes. The 
variations in the measured transmission losses, however, are somewhat 
larger than predicted. 

out frequency in Figure 22 and through- 

t NAS3-19435 in order to facilitate comparison with one- third-octave band 
spectra characteristically utilized in the analysis of combustor noise. It 
is more instructive to evaluate these turbine test results on a linear 
frequency scale for current purposes. Such a linear plot is shown in 
Figure 23, along with a prediction made with the frequency inverse energy 

""""sy distribution model bialel 

the available acoustic energy into the highest cut-on mode and, in fact 
assigns all the incident energy to a new mode at the instance of cut-on (f 

propagation angle, 6i = ±90°, almost ensures complete 
reflection and therefore infinite transmission loss. The program assumes a 
more reasonable finite value of 20 dB at this point. The pred^^Ior^rbe 
seen o ^ very close agreement with the measured data, particularly in 
Filled-ir^d^r^ increased transmission loss points due to modal cut-on. 

Filled in data point symbols in the figures denote masking of the downstream 
tone by broadband noise. Hence the actual transmission iLs wafat I^LI as 
much as shown by such a symbol, but it could have been significantly higher. 


that rha ^ ^ between the data and saw-toothed prediction implies 

that the apparent scatter in the data about the "mean" bathtub floor was 
in part, a manifestation of a modal cut-on, due to asymmetric noise injection 
and not a random error in the measurements. The fact that greater scatter 
was observed for the low pressure turbine data (see Figure 24) than for the 
high pressure turbine data provides further verification of this thesis 
Because cut-on occurred earlier in the low pressure turbine as a consequence 
ILI % therefore, be more cut-ons over ^ 

L -t^ associated with the first modal cut- 

on is obvious in the high pressure turbine transmission loss spectlrLl 
conspicuous by its absence from the low pressure turbine data Lly because 
movinrri apparent during the testing and was carefully avoided by 

emphasized hele'^ tLt f requenc.^i. It was recognized then, and is ■ 

phasized here, that the prominence of the cut-on effect in the test data 

was most probably due to tbe method of sound injection into tbe turbines 

losf snectrim sib distribution and a flatter transmission 

loss spectrum such as shown by the solid line in Figure 22. 


64 



Turbine 
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Figure 22. Comparison of Data and Prediction Using Equal Energy Distribution. 




NASA Core Turbine 



Figure 23. Comparison of Data Prediction Using Frequency Inverse Energy Distribut 














! influence of the modal cut-on phenomenon in the test 

Hz or\t t extension of the bathtub floor to 2500 

Hz, or greater, from the original 1200 Hz used in NAS3-19435. The gradual 

frequencies above 3000 Hz coulS bfa^tributed 
3 2 T^P "" f discussed in the finite-chord modeling of Section 

3.2. Tlie increase in transmission loss at the very low frequencies has 

the exhaust nozzle, not the turbines. prooaDxy be due to 

transmission losses represent the bathf-nh p ..v measured 

as discussed earlier In Sectior4 1 I? ! 

mission loss at design point for the NASA core turbine would be 7.2 d? 
low prJss^re'LSl”:. '"iherri's^S It configuration of the 

above choking (see Figure 19) tHp aiff ® decrease in transmission loss 
poratlon of ^he splnnLg^Jd’i; Jifc'-r^tlrSiciLf 

the ;S^g;r«Ln“(s:r"l\“ i^6?\ comparison for 

multistaging analysis wherein the influence of adiacenrbrr^ 
included. The predicted transmission loss if of thrsLp a T 

measured: 10 dB. The data Hn ^na■fr.o^- ^ order as that 

between the 2 0 and s 9 dicate a small Increase, about 3 dB, 

deslgr '•'P P«==ooe ratios, but only for speeds other tian 


theorrin'’FlLrr27!‘" Bofl! compared with the new 

with theory. As in the case of the^”' Inlet data show agreement 

slight Increase In transmission loL iuJ'^reSure r^irr^' ““ 

decrease above choking is reproduced. ^ assure ratio below choking and 

of lhri^L^rfct?rSi:^:ci”:f“:dlaL„1 bit' incorporation 

model provided the critical element ^ '’'P actuator-disk 

the available datl alcmants to successfully explain the trends In 
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Figure 27. Comparison of Theory and Data for the 

Single-Stage High Pressure Turbine, Hot 
and Cold Inlet Flow. 
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5-0 USE OF THE THEORY AS A WORKTNn TOOL 


5.1 CONCEPTUALIZATION 

bi„a "echanlsm behind low frequency noise attenuation by gas tur- 

nere first proposed by*R*°S‘ae pa«“ran^™ °bl' 

d-i1e!“ 

"L— f— 

ft^-s^baiic 

staging. „ultl„ode programs L gr>™ b;io^ °/ ‘"f " 
listings can be found in the apSLe^°" descriptions and 

low freque^rno^'srtL'Ziss^^rt^^ ^he 

design cycle deck. Section sf information available in a preliminary 

mations that lend themselves to "bacrof threnJSjf "‘'t^JrLlculSK^^ 

the sI^onTSecti'vi": this"prog;L?°""^"^“^" -- 


5.2 COMPUTERIZED PREDICTION 

for timLsharrSagr ^ThrSsiJ''f 3^^ Appendices are in t^ORTRAN and written 
input required is shown in TablP VT^ “ Figure 28. The 

absolute flow angle, wheel speed (in thf 038^0^ a^L^atinf hi ^^ionity. 

is a little more difficult to arrive ^®®lgner. Off-desxgn Information 

suggests very little chance In i-Lp t- * a ^J^'^nately the available evidence 
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M0DMLT Flow Chart 


Input: No, of Stages (N) 

Performance Parameters 

(U, 0, P , T^) — (2N + 1) Times 

Tip Radius (R), Hub/Tip Ratio (a) 

•" t " 


Next 1/3 OB Frequency 


F = 50 to 4000 H 


Compute: For a Given Frequency (F) 

No, of Cut- on Modes (NTT) 
Energy Distribution 
[NTT Limited to lOl] 


I 


[ J = 1 to NTT [ - 


Next Mode 


Compute; For a Given Mode (J) 

Incidence Angle (0j) 

Relative Mach Nos. and Flow Angles 

(M ,M ,M ,M ,«,3) 
nx ny mx my n m 

Wave Angles and Ratios 

p /P , k /k , k /k , k /k , k /k 1 

m n nx n ny n my m my m-* 


I 


Test for Propagation and Total Reflection 

V 

Form: Matrix Elements [A. .1, PD 1 

Invert Matrix [D] 

Compute Matrix [TC] = (D^^)(A2fi) (d/)(Aj^) 

Compute Transmitted Wave Amplitude (T) 

Compute Transmission Loss [TL(J)] 

I 


TL = 20dB 


No 


Next Mode 


J = NTT? 
Yes 


Compute: Summed Transmission Loss for All 

Modes TL_ 


i 


I F = 4000 HZ? h 
^ Yes 
pNext Case? 


No Next 1/3 OB Frequency 


Figure 28. Flow Chart - Multistage, Multimode Computer Program. 
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an isolated blade-row mode by specifying zero number of stages; the program 
will then faithfully reproduce the results of the previous published computer 
programs (References 3 and 8) . 

The program starts with the lowest specified frequency and computes the 
number of cut-on modes, the energy distribution, the equivalent plane-wave 
incidence angle, the transmission loss associated with each spinning mode, 
and finally the summed transmission loss for that frequency. The frequency 
distribution specified in the program is the center frequencies for the one- 
|-]^£].^_Qctave bands from 50 to 4000 Hz. However, this can be changed very 
conveniently to any other frequency distribution, for example, the siren tone 
frequencies from NAS3- 19435. 

The most important frequency is the first modal cut-on and this is 
calculated and printed out using the inlet mean radius. The transmission 
loss at this point exhibits a sharp spike, and the subsequent values of 
transmission loss register a significant increase as shown in Figure 29, 

^Ijxch corresponds to the input of Table VII. The transmission loss below the 
cut-on frequency corresponds to the axial plane wave only and is 1.5 dB 
in the example shown. The loss above the cut— on is controlled by the spinning 
modes and levels off at 9.5 dB. 

Equal energy distribution is specified in the two computer programs. 

This can be changed to any other desired energy distribution, including 
frequency inverse, as Indicated in Appendix A. The differences between equal 
energy and frequency energy distributions have already been discussed. The 
latter gives prominent spikes at each cut-on frequency. The height of each 
spike will depend on the assymmetry of the source: a line source giving equal 
values for each spike. A symmetric source, such as an annular combustor 
arrangement, probably will result in rapidly diminishing spikes. Whether 
these spikes will be discernible in broadband combustor noise spectra remains 
to be seen. It may be possible to use very high resolution (narrowband) 
analysis to detect the modal cut— on defects in the transmitted combustor 
noise spectrum in the exhaust nozzle. Also, the cut— on phenomenon could 
diffuse over a wide frequency band due to viscous effects, random flow 
variations, or "soft" duct walls. Cut-on for fan noise has indeed been 
observed to be a diffuse rather than discrete frequency phenomenon. Some 
ciarif ication may b® provided by the results from CF6— 50 tests now proceeding 
under NASA Lewis funding (ECCP III, NAS3-19736). 

It should also be recognized that turbine area, and therefore mean 
radius, generally will increase proceeding downstream. At the same time the 
static temperature will decrease. The cut-on frequency is proportional to 
(a/r*), and therefore will also decrease. This will not only contribute to 
the diffusion of the first cut-on spike, but also will mean a sudden shift 
to a lower frequency in the case of a gooseneck between high and low pressure 
turbines, as found in the CF6 family of engines. The Investigator may prefer 
to use the radius and hub/tip ratio downstream of the gooseneck in the 
computer program instead of the high pressure turbine inlet values. These 
dimensions are used to compute the cut-on frequencies and for no other purpose. 
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6~Stage 



Figure 29. Typical Transmission Loss Spectrum* 








5 . 3 APPROXIMATE ESTIMATION OF THE TRANSMISSION LOSS 


The computer prediction methods were parametrically exercised for a 
number of different aircraft engine turbine systems (Table VIII) in an effort 
to discern trends and simplifications that could be used in a semiempirical 
prediction technique. The net outcome was the prediction spectrum shown in 
Figure 30(a). The transmission loss for frequencies below first cut-on is 
constant, corresponding to 0i = 0. Then at f = f^, the loss increases to 
10 dB, which represents a mean value obtained using equal energy distribution 
for multistage turbines. This value will, of course, be higher in the case 
of an as}nnmetric source. The maximum value indicated by the test data is 
20 dB. For frequencies higher than fc> the transmission loss decreases to a 
value somewhat below a final asymptotic value which is attained with a small 
jump at the second cut-on. 

The first cut-on frequency is clearly the most crucial element here 
because the variation in the flat part of the transmission loss spectrum is 
fairly small from turbine to turbine. Figure 30(b) provides a convenient 
method of estimating this frequency given the mean turbine radius and static 
temperature. The inlet axial Mach n umber is assumed to be 0.3. The Mach 
number correction is actually /l -M^^; higher Mach numbers result in lower 
cut-on frequencies. 

In general, the transmission loss below the first cut-on is very low 
(5 dB or less). Therefore, a small turbine would offer little resistance to 
the transmission of peak combustor noise levels which, it is generally 
accepted, occur near 400 Hz for current engines. For example, the turboshaft 
engine turbine system used in the study (Table VIII) will induce only 3.2-dB 
transmission loss below 1350 Hz because of its size. The predicted trans- 
mission loss is shown in Figure 31. It is interesting to note that an engine 
(core noise) data correlation using combustor source noise parameters col- 
lapsed the available data along two lines as shown in Figure 32 (Reference 
2). Comparison with the component data line suggests much lower overall 
transmission loss for the three turboshaft engines than the turbojet and 
turbofan engines. One of the obvious differences is the exhaust transmission 
loss due to the nozzle, and flow is much lower for turboshafts. The other 
difference is that all three of the turboshafts in Figure 32 were very small 
engines and would have turbine transmission loss spectra similar to that 
shown in Figure 31. 

Table VIII suggests that 9 dB is a good value for the f > 2f^ asymptotic 
part of the transmission loss spectrum. Keeping in mind that the frequency 
range of interest for combustor noise is normally below 2000 Hz, a constant 
value of 9 dB above 2f^ would result in less than ±1.5-dB error for the 
turbines in Table VIII which encompass a very wide range of variables. 

The transmission loss below cut-on is defined by Figure 33. The loss 
actually decreases with pressure ratio and Mach number for multistage turbines. 
The reason for this, and the constant asymptotic value for f > 2fc, lies in 
the influence of upstream blade rows on the reflected upstream— propagating 
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4-stage Turbine System (Turboshaft) 
Equal Energy Distribution 



Figure 31. Transmission Loss Spectrum for a Turboshaft Engine Turbine. 
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Figure 32. Engine Data 
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waves from downstream rows. High pressure ratios cause almost complete re- 
flection at upstream rows and almost total restoration of the upstream- 
propagating energy to the downstream direction. In fact, there can be no 
upstream transmission of acoustic energy through a choked blade row. On the 
other hand, blade rows operating with small pressure drops will permit two- 
way transmission of the acoustic energy. The net effect for turbines having 
two or more stages is increased transmission loss at low pressure ratios. 

Figure 25 is recommended for predicting single-stage turbine transmission loss. 

Finally, the data suggest that the value between f^ < f < Zf^ is 
approximately two-thirds of the difference between the asymptotic and axial 
plane-wave values. That is, if the asymptotic value of the transmission loss 
is 9 dB and that for 0j = 0 is 3 dB, the value between fc and 2fp should be 
taken to be 7 dB. 

The above described approximate method of estimating turbine transmission 
loss is summarized below. The intent of the procedure is to generate a trans- 
mission loss characteristic such as the one shown in Figure (30a) for a 
specific turbine design. 


Turbine Transfer Loss Approximation Procedure 


Determine f^, from f^ = 


A - M ^ 

x_ 

2tt 


where f^ = turbine cut-on frequency 
M = turbine inlet Mach number 

X 

a = turbine inlet speed of sound based on inlet 
static temperature 

r = turbine mean radius 

(Note that Figure (30b) shows calculated values of f^ for 
Mx = 0.3). 

The value of TL (transmission loss) at f^ will be 10 dB 
or more, depending on source symmetry as discussed in the 
first paragraph of Section 5.3. The potential effects of 
a gooseneck transition between high and low pressure tur- 
bines should be considered here, as discussed in Section 5.2. 

At frequencies above f = 2f , TL is 9 dB for a multistage 
turbine, and is determined from Figure 25 for a single stage 
turbine. 


3. At frequencies below f = f , TL is determined from Figure 33 
as a function of the exit fiach number from the first blade 
row. 
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4. 


The constant value of TL in the range above and below 2f(, 
[referring to Figure (30a)] can be estimated to be 2/3 of the 
way between No. 2 and No. 3 values determined above. 

5. A transmission loss spectrum similar to Figure (30a) can now 
be drawn for the specific turbine design being evaluated. 


It is important to remember that this procedure yields only the turbine 
transmission loss. The transmission loss through the exhaust nozzle can also 
be an important consideration for gas turbine engines, particularly turbojets 
and turbofans. A fuller discussion of exhaust nozzle transmission loss can 
be found in Section 2.4 of Reference 2. Briefly, the loss can be modelled as 
a transmission loss, due to flow changes at the exhaust nozzle and through 
the jet(s), and a radiation loss, due to passage of the acoustic wave from a 
duct into open space. Classical analysis of the latter suggests that this 
part is negligible except for nozzles with characteristic dimensions very 
small compared to the acoustic wavelength. This is not usually the case 
except for very low frequencies. The transmission loss part postulates the 
same mechanism, specific impedance, and Mach number discontinuities as used 
i^ bhe turbine blade-row transmission modeling. A closed— form solution can 
be obtained for axial flow and = 0 and is given in Reference 2. A chart 
is shown in Figure 34(a) for the total loss proceeding from inside the 
exhaust nozzle to ambient conditions. The computer prediction programs were 
used here to generate the convenient chart in Figure 34(b) for the transmission 
loss due to changes in the flow through the exhaust nozzle. The effect of 
higher order. modes is included. This chart defines the exhaust nozzle trans- 
mission loss as a function of the temperature ratio across the nozzle. 
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6,0 CONCLUSIONS 


Two theoretical models were presented to describe the transmission of 
low rrofjuency noise through aircraft engine turbines. The somewhat complex, 
finite-chord analysis indicated that the simpler actuator-disk analysis was 
valid for frequencies as high as 0.4 to 0.5 of the blade passing frequency. 

In essence, it meant that the simpler model was adequate over the entire 
frequency range of interest for combustor noise. It was shown that multiple 
blade row and spinning mode considerations also had to be introduced into the 
analysis in order to fully explain the transmission through the blade rows. 
Interaction of acoustic waves with turbine blade passage shocks was found to 
be a very weak, second-order effect. 

The improved theory demonstrated that turbine transmission loss spectrum, 
in the midfrequency range (200-1200 Hz) was indeed flat as postulated by the 
floor of the bathtub spectral fit to the data of NAS 3- 19435. The scatter in 
the data about this floor was found to be due to higher order mode cut-ons and 
a biased energy assignment because of the assymmetric sound injection. Also, 
the flat transmission loss apparently extended to 2000-2500 Hz. Diffraction 
by the turbine blades was probably responsible for increasing the trans- 
mission loss at higher frequencies, giving one end of the bathtub spectrum. 

The other end of the bathtub, the rise at very low frequencies, was shown to 
be a spurious effect introduced by the location of the downstream sensors at 
a pressure cancellation point. 

The theory suggested a step function type of transmission loss spectrum, 
with the jump occurring at the first modal cut— on frequency. The attenuation 
below this frequency was predicted to be constant (with frequency) as would 
correspond to the transmission loss associated with the axial plane wave 
alone. This value would vary from about 5 dB to 1 dB, decreasing with 
pressure ratio for multistage turbines. The transmission loss (prediction) 
for frequencies greater than the cut-off frequency was found to vary betvjeen 
(9 ± 1.5) dB, independent of the number of stages or turbine pressure ratio. 


The first cut-on frequency, which is inversely proportional to size, 
appeared to be a critical element in the transmission loss since combustor 
noise generally peaks in the vicinity of 400 Hz. In particular, small turbo- 
shaft engine turbines would suffer since the step jump to the 9-dB asymptotic 
value is delayed to beyond 1000 Hz. Turbojets and turbofans should exhibit 
higher transmission losses, not only due to earlier cut-on, but also because 
of higher losses at the exhaust nozzle induced by accelerating flow. 

The exhaust nozzle and turbine transmission losses were computed 
separately and independently. It should be a fairly straightforward matter 
to link the two wave systems with due consideration being given to the phase 
change between turbine exhaust and core nozzle exit. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


The analysis performed in this contract has provided two of the four 
modules required for the modular prediction of combustor noise. The work 
proceeding under another NASA Lewis Contract, NAS 3-19736, wherein combustor 
noise is being measured at the source and various locations in a CF6-50 
engine, should further the activity. 
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APPENDIX A 


MATRIX INVERSION COMPUTER PROGRAM 


The transmission loss for an unchoked turbine can be determined exactly 
using the matrix inversion procedure outlined in Section 3.3. The method has 
been programmed for time-share usage and a FORTRAN listing is provided in 
Figure 35. 

A flow chart for the computer program is given in Figure 36. The 
program reads the input parameters and then, for each of the one-third 
octave band center frequencies from 50 to 4000 Hz, calculates the cut-on 
modes (the axial plane wave is always cut-on) . Equal energy is assigned to 
each mode. An equivalent plane wave (see Equation 62, and text, page 39) is 
defined for spinning mode along with a corresponding incidence angle. 

The angles for the forward- and backward-travelling waves are calculated 
at each blade row and, if total reflection occurs or if there is no forward 
propagation, the transmission loss for that mode is set to 20 dB. Otherwise, 
the (A) and (D) matrices are formed, (D) inverted, and (D’^A) computed. The 
transfer coefficient matrix is obtained by multiplying the matrix product for 
all the blade rows and the transmitted wave amplitude T extracted. The 
transmission loss for that mode is then found. When this computation has 
been effected for all the modes at any frequency, the transmitted waves are 
weighted according to the energy distribution and summed to define the trans- 
mission loss at that frequency. 

The frequency and energy distribution can be redefined as required by 
the user. The working frequencies are listed in lines 310 through 330 and 
the energy assignment imposed in line 1570. For example, if frequency 
inverse energy distribution is desired, line 1570 is changed to: 1570 E(J1) = 
l./SQRT (1. - (FRSQ**2)), where FRSQ is the ratio of the cut-on frequency to 
the working frequency. 

The turbine tip radius and hub/tip ratio are input in line 300. 

Normally, the values at the high pressure turbine inlet are used. However, 
the values at intermediate stations, such as the low pressure turbine inlet^ 
may be more advisable in case of large variations in tip dimensions. 

An input sheet is shown in Figure 37. The performance data are stored 
in a data file and the name of this file inserted when requested by the 
program. A typical output is given in Figure 38. The print-out includes the 
input parameters, the cut-on frequency, and the transmission loss for each 
frequency. The number of cut- on modes at each frequency is also shown, and 
the angles for the incident, reflected, and transmitted acoustic waves; the 
amplitudes of the transmitted and first reflected waves, the energy fraction, 
and the transmission loss are provided for each mode. 
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*'*’*'■*'*'*"*'*■'*■* MUUMLl ■*■*■■*■*■*****■*■ 


00 10*frKUiJr! *; 


0O20*#L I ARY /MT I N V » R= ( UL I 3 ) USEHL I B/TUS » R 
0030C ***** FTLENAME MOUMLf ***** 

<■^0 COMMON y CAN GR / P I , TOUEG, TOtfAD 

nObO COMMON '/Cl N PUT/ NSTAGE, lOPT, I PRINT, PTO, TTO,STAGEX ( 5 , 2 1 ) , 

0060& NSTA, TITLE 

00/0 CHARA»TER TITLE*40 

OOdO COMMON i^CUTOFA/ THCL.THCU 

0090 COMMON /CAERO/ V (2 1 ) , MX i 21 ) , MY (2 1 ) , AMQAN (21 ) , RHDRATC 21 ) 

(H00& MACH (21), AS (21) 

Olio REAL KN QK M ( 2 1 ) , M X , ,M Y , MACH 

0120 INTEGER FRET) (20) 

'1130 EOU I VALENCE (KNQKM , AMOAN ) 


0140 

OIdO 

01o0& 

01 /O 

OldO 

0190 


COMMON 

COr/ii'/j[]N 

CO wIhON 
COi - . .jN 


/CAEROI/ U(2] ),PHI(21),VR(21),p(21),T(21) ,GAV(21) 
/cuur / TLOSS(IOO), 1'HI(100),THR(1U0) ,THT( ’ .)), 
Q(100),B(100 ,TK(100) 

/CMATRX/ 0(3, 3,21 ),lC (3, 3,21 ) ,A(3,3,21) ,PR0D(3,3) 
/CAITCH/ CFl , CF2 , CF3 , oUF ( 380 ) 




0200 

0210 

0220 

0230 

240 

02dO 

0260& 

0270 

0280 

0290 

0300 

0310 

0320a 

0330& 

0340 

03bO 

0300 

0370 

0380 

0390 

0400C 

0410 

0420 

0430 

0440 

04b0 

0460 

0470 


Ol MEi)S I . iN SI AGEP (10b) 

OIMEiLSIijN £(100) 

EQUIVal. NCE (STAGEP(l), STAGE/ (1,1 )) 

CHARACTER CFl*l/"/"/,CF2*8, CF3*l/« ; '•/ 

CHAhACTER TrfLE*40, BLANK*40 

REAL MACHN,MACHM,KMYK/-:. KNY<.-M,KMXK'”*,KNXiC.>i,iCMYSAV,KMXSAV, 
KNN,KMM 

EOU I VALENCE ( I B I T5 , b I I'S ) 

GAiA f5l rS/U3777777 1 1111/ , ■v/34037()00000o0/, JPO/(J04007b040007/ 

OAiA m J, SIGMA/16. 8, 0.889/ 

OATA FpEQ/bO,63, 80, 100, 12b, lo0,200, 

2b0, 3 15, 800,630,800, 1000, 

1250, 16C ), 200(1, 2500, 3150,4000/ 

OAiA PI , rUDEG.TORAJ/3. 141592 i ,57.29578, .01 /453292b/ 

OAi'A BLANK/" „/ 

namelist /TiJOISE/ IuPT,PTU, lT0,STAGEP,TirLE,(3AM, IAErJ 
i'An(a)= SIN(X)/Cu (a) 

SET Up namelist INpU'i FILE 

I AERO = 0 

CALL FPARAM(3, JPO) 

PRiN'i," INPUT FILE NAME " 

READ,CF2 

CALL ATI ACH ( 1 , CF 1 , 1 , 0 , STAT, BUF) 

IF( STAT.EQ.O. .UR. STAT. EO. OKA ) GO TO 5 


Figure 35, Program Listing - Matrix Inversion Program. 
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04tiO 
0490 
0500 
0510 
0d20C 
0530 
0540 
O'jbO 
0560 
0570 
0580 
05 yO 
0600 
OolO 
062u 
0o30C 
0640 
0650 
0660 
06/0 
0680 
0690 
0700 
0/10 
0720 
0730 
0 740 
I /50 
0760 
0//0 
0 /80 
0/90 
0800 
08lO 
0820 
0830 
084 ' Jlx 
0d::^(, 
0860 
08/0 
0880 
0890 
0900 
o9lO 
0920 
09JJO 
0940 
0950 
0960& 
09 /L 
0980 


PRINT 1,STAT 

1 Fl3i/MAT(" INPUT FILE STATUS=" , U12 ) 
STOP 


I N I f I AL I ZATI ON *■*■■*•***•***■**'*■•*'■*■*■*■*"*■■* 


o IOf'L = 1 
JO 10 1=1,21 
STAOhX(l,I)= BITS 
oAN(I)= BITS 
10 CONTINUE 

PTO = 14.690 
no = 518.7 
ilfLE = BLANK 


*s 

OF POOR 


hEAO INPUi FILE **** COUNT NP. OF STATIONS 


15 CALL RHAONA{l,TNOISL,"$TuOISE",JEND) 

IF( JENU.EQ.O ) GO TO 400 
00 17 1=1,21 

I F ( STAGEX ( 1 , I ) . EQ . B ITS ) GO TO 1 8 
U ( I ) = S i'AGHX (1,1) * 3 .048 

PHI (I) = STAGEX (2, 1) 

7R ( I ) = STAGEX (3,1) * 3 . 048 

P(I) = STAGEX (4, I) / 6.895 
l'( I) = STAGEX (5, I ) * 1.8 
IF( GA,/1( n.NE.BITS ) GO TO 17 
GAM(I)= GAMX(Td)) 

1/ CONTINUE 
18 NSTA = I-l 

IF( NSTA. EO. 2 9 ) NSTA = 21 
NSTAGE= (iJSTA- )/ 

PRINT 21 , TLiLr-, ; IaGE 

21 FOi64AT(//16X,A4u//32X, 12, •• STAGES"//) 

PRINT 22 

22 F0HMAT(28X, AERO- xHEtJMO PARAMETERS *"// 

2X, "STAGE", 3X, "SfATION " , 3X, "U- FPS" ,3X, "PHI- UEG" , 
3X,"Vk- FPS",2X,"PS- RSIA",2X, "TS- JEG R"/) 

NGlG = 0 
JO 24 1=1, NSTA 

1F( ( (1/2 )'-^2)/I.EQ.O ) NSTG=NSTG+1 
IF( I.EQ.NSfA ) NSlG=lblTS 

PRINT 23,NSfG,I,U(I),PIII(I),VR(I),P(I),T(I) 

23 FUHMAi (4X, I1,7X, I2,F12.3,4F10.3) 

,_4 COiJTINUc 

IF( lAERu.EO.O. ) GO Ti, 2/ 

PRIiKI 25 

25 F0RMA'I(/2X, "S'lAUE", 3X, "STATI ON", 5X , "MX" , 8X , "MY" , 

/ X , " M AC H " , 5 X , " K N QKM " , 7X , " V " / ) 

NSiG = 0 

:/0 26 I = l,NS'iA 
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lui^i 
1020 
1030 
1040 
1013 Od 
1060 
1070C 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
llbO 


n/2)*2)/I.EQ.0 ) NSTG=NSTG+1 
■y .I*HQ-_1 ) GO TO 26 

26 CONTINUE ^ ^ ^ ^ »MACH( I )»KNQKM( I ) , V( I ) 

2/ CONTINUE 

2o FONMAT ( //2X, ^TilETA-I ^ 3X, 'THETA-N%3X, 'THEIA-T', 6X. 
'i 9X, 'B', 9X,'E^, 6X, ^T-LOSS'/) 

CALCULATE AEkO-THEIJNU PANAMETENS 

00 29 I=1,NSTA 

AS(I) = 41.42*SQPT(GAM( I )*T( I )) 

MX(I) = U(I)/AS(I) 

IF( I .EQ. 1 ) (j(] Tt] 29 

AMQAi J ( I ) = bOli: 1 ( GAM ( I ) *T ( I ) / ( UAM -1 ) *1 ( I - 1 ) ) ) 
NJ1UHAT(I)= T(I)vcB(I-i,/(Y(I-i)*^. 'I )) 

29 CONTINUE 


1160 
1170 
1180 
1181 
1 182 
M84 
1 18b 
118o 

1187 

1 188 
1189& 
1190 
1200 
1210 
1220 
1230 
1240 
12bO 
12o0 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340C 
13b0 
1360 
1370 
1380 
1390 
1400 
1410 


ASl = AS(1) 

NMEAN= P0'-V5qhT( ( 1 .+8 IGMA’t^P) /2. ) 

AMEAJj = NMEAN/2.b4 
FPl = 1. 

F hEQCO= ( ( F P 1 * ( as 1 * 1 2 . ) ) / ( 2 . *P I *PMEAi\ ) ) *SQPT ( 
FNEQC1= AINT(FltEQCO) 

PHI NT 30,FPE0C1 
PRINT 28 


-X'U**2) 


30 FUNMAT(//10X, ^-Jr**** FIRST CUT-ON OCCUi^S AT',lX,Fb.O 
IX, ■'HZ ///) ’ 

00 300 L=l,20 

FP = (2.*PI*FPEQ(L) *PMEAN )/(ASl*12. 1 
XjM 1 — mX ( 1 ) 

NTH=FP/SQN T( 1 . -XM1**2) 

1F( NTH.GT.50 ) NTH=bO 
NTT = 2*NTH+1 
FPSQ= o. 

THI(1)= 0. 

JO 32 J=l,NiT 
E(J)= 0. 

32 CONTINUE 

F(l)= 1. 

ESIGmA= 1. 


THI(1)= 0. 

Ir( NTH.LT. 1 ) GO TO bO 
■kk-kk compute CUT-ON MOOES, ANGLFS, 
00 40 J=1,NTH 
FJ = J 
F1=FF/FJ 

F2 = SQPT(F1**4-: i*x2*( l.-XMl*^t2)) 
F3 - XM 1 **2 -<-Fl 2 

F4 = (F2-XMD/F3 
J1 = 2*J 


ANO ENERGY 


Figure 35. Program Listing - Matrix Inversion Program (Continued). 
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142U J2 = 2*J+1 

143U i'HI(l)= U. 

144U lHI(Jl)= T0UHG*Ai^CaS(F4) 

14bO THI(J2)= -THICJl) 

1460 FC = (FJ*AS1*12 . )/(2.*FI*W4HAN) 

Mob FC = FC*S0HT(1.-XM1**2) 

1470 FWSO = l-C/FFbQ(L) 

14«0 IF( FIlSO.Gi'. 1 .02b ) GO TO 40 

16/0 F(J1) = 1. 

I6b0 F(J2)= H(Ji) 

IbVO t:SIGMA= FS1GMA+2.^H( J1 ) 

1600 40 CUNTINUF 

1610 bO COiJlINUE 

1620 

1630C CUMFUTF HNHHGY OlSi^IBUiION 

1640 UL oO K=1,J2 

16b0 ciK)- H(K)/FSIGMA 

1660 60 CONiINUH 

lo70 
16 bO 

16900 ■*■*•*■*** Iiji\lhh! LOOP TU BUILD MATbICcS 

17 00 

i/lo SUMT= 0. 

1/20 DO 18b K=l,iFlT 

l/3o 62 fHFN = THICK )*fOAAD 
1/40 </viYGAV- 0. 

17bO n,..iXSAV= 0. 


. 760 

M 

= 

1 


1 / /u 

6b M 

- 

M+l 


1 /GO 

IF( 

M • G i 

’.NSlA ) Gu iu loo 


1 /90 





IbUOC 

CALCuLATb 

1 ANGLES AmJ h'Ali.^3 


1810 





1820 



M-1 


1830 

V(H) 

- 

U ( N ) *1 AiJ ( iUPAD^Ph I ( N ) ) 

-VN(N) 

1840 

\/(M) 

= 

U ( M ) *TAN ( iUAAD*Ph KM)) 

-VN(N) 

IbbO 

MYCN 

) == 

V(fj)/AS(N) 


I860 

HYCM 

) = 

V ( M ) /AS ( M ) 


18 /O 

MACH 

(N) = 

SQPl (MX (i. ) *m-2 + MY(N) **2) 

1880 

MACl I 

(M) = 

SQin (MX (M ) **2+/.'.Y (M) *-*2 ) 

1890 

N 

- 

M-1 


1900 

GA = 

GAM 

( I ) 


1910 

GB = 

(GA 

+ 1. )/(2.*(GA-l. )) 


1920 

AAST 

AR = 

( (2.+CGA-1. )*MACH(M)**2)/(GA+1. ) ) **GB/MACH ( M ) 

1930 

XMN 

= 

l.-MXCN )Vr*2 


1940 

XPiJ 

= 

l.+MX(f])**2 


19bO 

X MM 

= 

1 . -MX (M ) *>^2 


1960 

SINN 

z; 

SIN Cl HFN) 


19/0 

COSN 

= 

CDS ( THEN ) 



19bOC **** CHECK FOk UfSTEEAM FHOPAGATION 
1 990 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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V 


2000 
2010 
2020 
2030 
2040 
2050 
2060 
20 70C 
20b0 
2090 
2100 
2110 
2120 
2 1 30 
2140 
2150 
2 160 
2170 
2160 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
22ou 
22/0 
2280 
22 VO 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
23VO 
r-a.ju 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 


HHSFU= U (N ) +AS (N ) *CuSN 
IF( HHSPD.Lh'.O. ) GO TO 175 

OMN = KNQKM ( U ) *S I NN/( 1 . +MX ( N ) *COSN+MY ( W ) *S I NN ) 
IN = XMN*SINi. 

ru - X PN *cas n + 2 . *mx ( n ) 

rH8N = ATAH2( l'N,To ) 

I EhM = -GMi J*MX ( M ) * ( 1 . -GMN’^tM Y ( M ) ) 

CHECK FUR TUxAL REF'LECTIUN 


RUCL- ( 1 . -GiV1N*m Y ( M ) ) ■**2-XMM*UMN **2 
IF( RDCL.LH.O. ) GO TO 1/5 
RADI CL= GMN*SQRT ( ROCL ) 

FN = -TERM+RAUICL 
/ C = ( 1 . -GMI,' Y ( M ) ) **2 -GMH **2 
THEM = ATAH2( TH,TU ) 

IF( iJ.NE.l ) GO TO 70 
THR(K)= T0L)EG*THBN 
/O THEM = ATAN2( TN-2.*TERM , iD ) 
MACHM = MACH(M) 
mACHW = MACH(N) 

ALE AN - ATAN2 ( M Y (N ) , MX (iJ ) ) 

BEiAM = AfAN2( HY(M) , MX(M) ) 

IF( MACM(M) .LT. 1. ) GO TO 71 
001' = COS(bE'i'AM) 

8M = BETAM 

3ETAM = ARCOSIAAgTAR^DOT ) 

BElAM = 6 IGJI ( 1 , 8 m ) i^bbTAM 
71 CONTINUE 
KMYKM = GMN 
UYKN = KMYSAV 

KMXKM = (1.-GMN*MY(M))/MX(M) 

KNXKN = KMX6AV 
KM/6AV= KMYKM 
•<;'USAV= KMXKM 
OKNR = 0. 


/5 


I E ( N . EQ , 1 ) GO ro 75 

OKNN = 1 . /SORT ( KN XKN**2+KN YKN**2 ) 

QKMM - 1 . /SQRi ( KMXKM**2+I04YKM**2 ) 

A ( 1 , 1 , N ) = M X ( N ) +COS ( THEN ) 

A ( 1 , 2 , N ) = MX ( N ) -COS ( THBN ) 

A ( 1 , 3 , N ) = KN Y KN *QK1FN 


1 . +MACHN*COS { ALFAN-THETJ ) 

1 . -MACHN*CUS ( ALFAN+THBN ) 

OKNN* ( MX ( N ) *KN YKN-MY ( N ) *KNXKN ) 

( MX ( M ) +COS ( THEM ) ) /AMQ AN ( M ) 

( MX (M ) -COS (THEM ) ) /AMQ AN ( M ) 

KM YKM*QKMM/AMQAN ( M ) 

RHGRAT ( M ) * ( 1 . +MACHM*COS ( BETAM-THFM ) ) 

RHURAl ( )* ( 1 . -MACHM*COS (BETAM +THBM ) ) 

RHORA f (M) * (QKMM*(MX (M) * kmYKM-MY (M) *KMXKM ) ) 


Figure 35. 
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251U A(3,2,N)= 0. 

2d2U A(3,J,i'j)= 0. 

2530 U(3,1,N)= SIN(BEi'AM-THFM) 

2540 0(3,2,N)= -S I U ( BETAM+THBM) 

2550 0 ( 3 , 3 , in = QK M;. ( KMY KM*S I iM ( BET AM ) +KMXKM*CUS ( BtiTAM ) ) 

2560 00 TO 80 

25/0 78 FUhMAT( 5X. ^*^*DUOiiOT[?EAM RELATIVE FLO/J AT RC//',I3, 

2680cx 1X,^IS f^UFEwSUN'IC***'///) 

2'.>yOC ***** COMPUTE liiVEliSE OF MATrilX AMD STORE 

2600 

2610 80 CALL OliJVER( N ) 

2620 I'HFN i'HFM 

2630 00 TO 65 

2640 

2650C ***** COMPUTE MATRIX PRODUCT 

2 660 

2670 100 CALL MAPROD( NSTA-1 ) 

2680 

2690C ***** STORE AMPLITUDES 

2 /OO 

2/10 110 1HT(K)= THFM*'fUDEG 

2/20 B(K) = -PR0L)(2, 1)/PR0D(2,2) 

2730 i'W(K) = PRODCl, 1)+B(K)*PR0D(1,2) 

2/40 0(K) = PR0D(3,1)+B(K)*PR0D(3,2) 

2 /50 

2/60C ***** COMPUTE Ti^ANSM I SSI ON LOSS 

2 //O 

’/80 120 ASl = AS(1) 

2790 aSN = AS(NSTA) 

2800 ^’X(l) = U(1)/AS1 

2810 Ma(NSTA)= U(NSTA)/ASN 

2 820 VVl = U (1 ) *TAN ( TOrAD*PH 1(D) 

2 830 VVN = U(NS'iA)*iAN (TORAD*PHI (NSTA)) 

2840 mY( 1) = VVl/ASl 

2850 1Y(NSTA)= VVN/ASN 

2860 RHORA = P(NSTA)*i( 1 )*AS U/(P( 1 )*T(NSTA)*AS1 ) 

2870 :ERM2 = (l.+MX(NSTA)*COS(THFM)+ 

2 8803. ( .,1 Y ( NST A ) + VR ( I' STA- 1 ) /ASN ) *S I N ( THFM ) ) * ( COS ( THFM ) +MX ( NSTA ) ) 

2890 iriIN = TH I ( K ) *'l ORAo 

2900 I'ERMl = (1.+MX(1 )*CljE(TH1 - J+MYXD* 

2910& SIN (THIU))*(COS( THIN )+MX(l )) 

2 920 FLOSS ( K ) = 10. *AL JOl ; ' ( iLMOHA*ABS ( TERMl / 1 ERM2 ) /TN ( K ) **2 ) 

2925 IF( TW(K).LT.O. ) GO TO 176 

2930 IF( TW(K).GE.l.) GO TO 1/6 

2935 GO TO 180 

2940 l/'j rL05S(K)= 20. 

2945 GO TO 177 

2947 17o TLOSS(K) = 25. 

2950 17/ B(K) = 1. 

2900 TN(K)= 0. 

29/^ 180 CONTINUE 


Figure 35. Program Listing - Matrix Inversion Program (Continued). 


95 


2975 IF( TLOSS(K) .LT.O. ) TLOSSOO = 1. 

2980 SUMI= 1. 

2 990 TL = TLQSS( 10/10. 

3000 SUMT= SUMT+E(10/10.**TL 

3010 185 CONTINUE 


3020 

3030C 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120& 

3130 

3140 


200 


235 

240 

245 

300 


I'LS 1 0'MA= 10 . *AL0G1 0 ( SUM I /SUM! ) 

***** PRINT OUTPUT 

JCT = K 
UG 240 1=1, KT 

^ »'THR( I ),THT( I) ,TW(I ) ,B( I) ,E( I ),TLOSS( I) 
PGRMAT(P9.3,2F10. 3,F9.4,3F10.4) 

CONTINUE 

PRINT 245,FREQ(L),TLSIGMA 
F0RMAT(/14X, 'FREQUENCY=' , 14, IX, 'HZ', 5X, 

^ TRANS MI SSI ON LUSS=', F6. 2////) 

CUNITNUE 


TG 15 


3150 


IF( JENO.NE.-l ) GO 

3100 



3170 

40u 

STOP 

3180 


END 

3190CGAMX 

FUNCTION GAM. 

3200 


F UN C LION GAMX(T) 

3210 


IF( T.LE.800. ) GO 

3220 


1F( 1. GE.3600. ) GO 

3230 


vjAMX. = 2.2 3/OB/T**, 

;240 


GO TO 15 

3250 

10 

GAMX = 1.4 

3260 


GO 10 15 

3270 

12 

GAMX = 1.254 

3280 

15 

RETURN 

3290 


END 

3300 




TO 10 
TO 12 
.070271 


33 lOCU INVER 
3320 


3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460C 

3470 


CALCULATE INVERSE OF MATRTX J 
SUBROUTINE DINVERC N ) 

COMMON /CMAfRX/ 0 ( 3 , 3, 21 ) , 01 (3, 3, 21 ) , A( 3 3 21) PROOf^ > 
JIMENSION DU(9,21),DDI(9,21),LABEL(3) ’ ’ ’ 

^^D(l,l),0(l,i,l )),(DDI(1,1),DI(1,1,1) ) 
L^IiMENSION PPR0LK9) ,TEMP(3, 3) , TEMPI (3,3) ,TEMP2(9) 
EQUIVALENCE ( PPUOD ( 1 ), PR0O( 1 , 1 )), (TEMP2 (1 ), TEMPI ( 1 , 1 ) ) 


NN = N 
00 10 1=1,9 
DDI (I,N)= DD( I,N) 

10 CONTINUE 

CALL MTINV(DDI(1,N) ,3,3,3,LABEL) 

20 RETURN 

**** entry MAPROO ** COMPUTE PRODUCT OF 01 


AMD A 


Figure 35. Program Listing - Matrix Inversion Program (Continued). 
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3480 


ENTRY MAPRGU( N ) 

3490 


NN = N 

35 00 


DO 30 1=1,3 

3310 


DO 30 J=l,3 

3520 


PROD(I,J)= 0. 

3530 


IF( l.EQ.J ) PRQD(I,J) = 1. >' 

3540 

30 

CONTINUE 

3530 


3560 


DO 100 L=1,NN 

35/0 


DO 60 J=l,3 

3580 


)0 50 1 = 1,3 

3590 


lEMP(I,J)= 0. 

3600 


DO 40 K=l,3 

3610 


TEMP(I,J)= TEMP(I,J)+DI ( I , K , L) *A (K, J, L ) 

3620 

40 

COuTINUE 

3630 

50 

CONTINUE 

3640 

3630 

60 

CONTINUE 

3660 


DO 90 J=l,3 

36 /O 


DO 80 1=1,3 

3680 


fEMPl(I,J)= 0. 

3690 


DO 70 K=l,3 

3700 


TEMP 1 ( I , J ) = TtHP 1 ( I , J ) +TEMP ( I , K ) *PROD (K , J ) 

3710 

70 

CONTINUE 

3720 

80 

CONTINUE 

3 730 

90 

CON'lTMUE 

3/40 



3/50 


DO 95 1=1,9 

3760 


PPROD(I)= TEMP2(I) 

3 770 

93 

CONTINUE 

3780 

100 

CONTINUE 

3/90 

200 

RETURiJ 

3800 

3810 


END 

3820CCUTUFF DETEIRvlIN LUivER/UPPER CUTOFF LIMITS 


3b3U SUdROUTINH CUTOFF 

Jb40 COMMON /CANGP / P I , i UOEC/fORAD 

38b0 CO.viMON /CUTJFA/ THCL,THCU 

3860 COMMON /CAEPO / V ( 2 1 ) , MX ( 2 1 , »MY ( 2 1 ; , AmQAN (21), 

38 /0& PM0PAT(21) , HACH( 21 > , AS (21 ) 

3880 REAL MX , MY , MACH , KNQKM( 2 1 ) 

3890 EQUIVALENCE (KNQKM , AMQAR ) 

3900 COMMON /CINPUT/ NSTAGE , lOPT, I PRINT, PTO, TLO, 

3910& STACEX(5, 13) ,NSTA 

3920 COMMON /CAEROI/ IJ ( 2 1 ) , PHI (21 ) , Vk ( 21 ) , H « ' . ) , T( 2 1 ) , GAm ( 21 ) 

3930 DIMENSION THCl ( 43 ) , TfiC2 ( 2 1 ) , THC3 ( 2 1 ),A.-.GF 02) 

3940 EQUIVALENCE (THC2 , THCl ( 22 ) ) , (THC3, THCl ( 43 / ; 

3930 DAI A BITS/0377777 //777Y/ 

3960 

3970 rAN(X)= SIN(X)/COS(X) 

3980 
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39900 
4 000 


***** CALCULATE I'HE 3 CUTFF INCIDENCE ANGLES 


4010 

4020 

4030 

4040 

40b0 

4060 

4070 

40«0 

4090 

4100 

4110 

4120 

4 130 

4140 

41b0 

4160 

4170 

4180 

4190 

4200 

4210& 

4220 

4230 

4246 

42b0 

i j 

42/0 

4280 

4290 

4300 

4310 

4 o2 OC 

4330 

4340 

43b0 

4360 

4 370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

44b0 

4460 

4470 

4480 

44 90 


00 b 1=1, 4b 
THC1(I)= BITS 
b CONTINUE 

00 100 I=1,NSTA 
II = I+l 

THCl (I )= . 90.+T0DEG*ARSlN(MX( I ) ) 

THC2(I)= -THCKI ) 

IF( I.EO.NSTA ) GO TO 100 

VVl = U(I )*TAN(T0HAD*PHI(I ))-VN(I) 

VV2 = U( II)*TAN(TORAD*FHI ( ID )-VH(I) 

MY( D = VVl /AS ( I ) 

MY( II)= VV2/AS( II) 

MACH(I)= SQi/r(MX( I) **2+MY( I )**2) 

MACH ( II )= SORT(MX (II) **2+MY( II ) ** 2 ) 

TERM = SQRT(1.-MX( II )**2) 

ANGFM ( 1 ) = ATAN2 ( TERM, -MX (II)) 

ANGFM(2)= ATAN2 (-TE1R/1,-MX( II ) ) 

KGO = 1 

10 ANG = ANGFM (KGO) 

GNM = CIN ( ANG) /(KNOKM( 1 I ) * ( 1 ,+MX( II )*COS ( ANG)+ 
MY( 1 1 )*SIN(ANG) ) ) 

X = (i*-GNM*MY(I))**2-(l.-MX(I)**2)*GNM**2 

IF( X.GE.O. ) GO TO 20 
IF( KG0.E0.2 ) GO TO 100 
KGO = 2 
GO TO 10 

20 XNU = GWM*MX( I )*( 1 .-GNM*MY( I ) )+GNM*SORT( X ) 

XDEN = ( 1 . -GNM*MY ( I ) ) **2-GNM**2 
iHC3(I)= TOOEG*A TAN2 ( XNU, XDEN ) 
lOO CONTINUE 


***FIND LARGEST - ANGLE AND SMALLEST + ANGLE 


THCL = -boo. 

DO no 1 = 1,44 

IF( THCKI ) .EQ.bliS .OR. THCl ( I ) . GT. 0. 
LHCL = AMAXK 1HCL,THC1(I) ) 
no CONTINUE 

THCU = 500 . 

DO lib 1=1,44 

IF( THCKD.EQ.brrs .OR. THCl ( I ) . LF.O. 
THCU = AM INK THCU,THCKI) ) 
m CONTINUE 


) GO TO no 


) GO TO lib 


120 RETURN 
END 

SUBROUTINE VIKI 

RETURN 

£Nl) 
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******* filename: mtin.v ******* 

IOC MTliNV 

30* *************MATHI X INVEiio lOi'l************************* 

40 bUbKJUilNE MTINV(A,NFARG,NCAhfG, DIM, LABEL) 
bO DIMENSION A ( IDIM,NCAAO ) , LALEL(nnAAO') 

60 1 iJrV=4HAt^G 
/U OC=)'«CA/LJ 
HO JJ ^1 Jl=l,Nk 
90 21 LAbEL( J1 )=J1 
loo JU 291 Jl=l,Nk 

IK,* **************F HEMAIiKiJO kOW CGN'iAllJINo LARGE > i *** 

l,:iu* **************AbSULUTE VALUE IH FIVdTAL CnLU-iU****’^*** 

13(j 101 1L-1P=G.0 
140 L.v, 121 J2=Jl,Nk 

IbO IF (AbS(A(J2,Jl )). L'i .VE ir ) G(J lu 121 
IbO iEMF=^ABS(A( J2, J1 ) ) 

1/0 IblG=J2 

180 1,^1 CUiITIhUE 

190 IFdbiG.EQ.JDGU id 201 

?00* ******** *****-/x-RE A kill AN OF -idWb IG PLACE LAOOcSi ABSGLIJIE 

210* **************VaLUl In PlVdF PGS III UN «■** <************* 

220 Jd 141 J2-1,NC 
230 FEMP-A(dl, J2) 

240 A(Jl,J2)^A(IblG,J2) 

2b0 141 A( I BIG, J2 )=TEi/iP 
260 I=LABEL(J1) 

2 70 LAbELlJl ) = LABE'L( IbiG) 

280 LABEL! IBlo)=I 

300* : : : :CdMP'JI'E CUEFFICIENTS IN PIVUTAL HUN:::: 

3U0 201 lE./lF-^A ( J1 , J1 ) 

310 A(Jl,Ji .0 

320 UG 221 J2-1,NC 

330 221 A( J1,J2)-A(J1, J2 )/iEhK 

340* **************C,d b-u i E CdEFFICIF:NTS IN GTHFil RUNS****** 

360 DU 281 J2=l,Nk 

360 Ih(J2.EQ.Jl) Gd To 281 

370 iEmH=A(J2, J1 ) 

380 A(J2,J1 )-0.0 
390 Ud 241 J3=1,NC 

400 241 A( J2. J3)=A( J2, J3)-iEMr^*A( Jl, J3) 

410 2b 1 CUNilNUE 
420 291 CGNlINUE 

430* **************INfEkCHAMGE CULUMNS ACCUkDING ru******** 

^40* **************INTEkCHANGE5 GF RUNS GF GklGlNAL MATkIX* 

460 30l iNl^Nk— 1 , 

460 DG 391 Jdl,Nl 
470 DU 321 J,.-=J1,NU' 

480 IF (LABEL (J2) .NE.Jl ) GG Id 321 
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'-ivu IK J2 J1 ) UJ ro 39 1 
t'’)U 3U id 341 
■diU 321 CUiJTlNUt: 

520 341 OU 361 J3-l,i'itI 
5 30 ih(.i2=^A( J3, Jll 
54U A(J3, J1 )=A(J3, J2) 
ddO 3ol A( J3, J 2 )='aiMi- 
560 LA3hL( J2 ) = LABh'L( Jl ) 

5 70 3vl CdiiiluUi: 
dbO 5001 HdTUBI'J 
5 90 dNO 


Figure 35. Program Listing - Matrix 

Inversion Program (Concluded). 



M0DMLT Flow Chart 


0-2 KXK ^ V 



Figure 36. 


Flow Chart - Multistage, Multimode Computer Program 
Using Matrix Inversion. 
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Data File 


Tip Radius (R®) - 
Hub/Tip* Ratio (a) - 


Turbine/Power - 
No. of Stages (N) - 

(use 0 for isolated row) 
Input (2N + 1) times: 


Axial Abs. Flow Wheel Static Static 

Stn Vel. (U) Angle (cj)) Speed (Vr) Pressure (Pg) Tc.mp (Tg) 
(m/sec) (deg) (m/sec) (kN/m^ or kPa) (K) 
























01 /\ 2/11 11.317 

INPUT FILE NAME = DFLPII f T' •, 

a? Mk v,;'.: ' ■ 

3-STG LPTs 2.0 PR 100% N 

3 STAGES 

* AERO-THERMO PARAMETERS * 


STAGE 

STATION U- FPS 

PHI- DEG 

VR- FPS 

PS- PSIA 

TS- DEG R 

1 

1 

607.000 

0. 

0. 

34.630 

730.000 

1 

2 

42 1 .000 

62.100 

409.000 

28.470 

693.000 

2 

3 

327.000 

-33.600 

0. 

26.170 

679.000 

2 

4 

300.000 

62.100 

428.000 

23.400 

660 . 000 

3 

5 

255.000 

-2.300 

0. 

21 .390 

648.000 

3 

6 

219.000 

55.100 

437.000 

20.280 

64 1 . 000 


7 

21 3.000 

45.630 

0. 

1 9.560 

639.000 


'k'k'k'k'k 

FIRST CUT-ON 

OCCURS AT 

219. HZ 

***** 



THETA- I 

THETA-R 

THETA -T 

T 

B 

E 

T-LOSS 

0. 

0. 

0. 

0.5267 

1 .5109 

1 .0000 

5.3529 


FREQUENCY= 50 HZ TRANSMISSION L0SS= 5.35 

0. 0. 0. 0.5267 1.5109 1.0000 5.3529 

FREOUENCY= 63 HZ TRANSMISSION LOSS= 5.35 

0. 0. 0. 0.5267 1.5109 1.0000 5.3529 

FREQUENCY= 80 HZ TRANSMISSION LOSS= 5.35 

0. 0. 0. 0.5267 1.5109 1.0000 5.3529 

FREOUEHCY= 100 HZ TRANSMISSION LOSS= 5.35 

0. 0. 0. 0.5267 1.5109 1.0000 5.3529 

FREQUENCY= 1 25 HZ TRANSMISSION LOSS= 5.35 

0. 0. 0. 0.5267 1.5109 1.0000 5.3529 

FREQUENCY= 160 HZ TRANSMISSION LOSS= 5.35 


Figure 38 


Typical Output 


0. 


0. 

0. 

0.5267 

1.5109 

0. 3333 

5.3529 

105. 

152 

0. 

0. 

0. 

1 . 0000 

0. 3333 

20.0000 

-105. 

152 

-51 .205 

0. 

0. 

1 .0000 

0.3333 

20.0000 



FREQUENCY= 200 

HZ TRANSMISSION 

L0SS= 9.84 



0. 

0. 

0. 

0.5267 

1.5109 

0.3333 

5.3529 

76.61 7 

0. 

0. 

0. 

I .0000 

0.3333 

20.0000 

-76.617 

-32.3J0 

53.120 

0.2976 

0.4328 

0.3333 

10.0400 


FREQUENCY 

= 250 HZ 

TRANSMISSION 

L0S5= 8.74 



0. 

59.477 
-59.4 77 

0. 0. 

0. 0. 

-23.662 -41.815 

0.5267 

0. 

0.3052 

I .5109 
1 . 0000 
0.7204 

0.3333 

0.3333 

0.3333 

5.3529 

20.0000 

10.4325 


FREQUENCY= 315 HZ 

TRANSMISSION 

LOSS= 8.84 



0. 

46.429 

-46.429 

105.152 

-105.152 

0. 

17.875 
-1 7.875 

0. 

-51.205 

0. 

40.278 

-32.994 

0. 

0. 

0.5267 

0.2935 

0.3463 

0. 

0. 

1.5109 
0.9117 
0.9371 
1 . 0000 
I .0000 

0.2000 

0.2000 

0.2000 

0.2000 

0.2000 

5.3529 

8.9339 

9.5133 

20.0000 

20.0000 


FREQUENCY= 400 HZ 

TRANSMISSION 

L0SS= 

9.58 


0. 

37.004 

-37.004 

76.617 

-76.617 

0. 

1 3.991 
-1 3.991 

0. 

-32.310 

0. 

30.951 

-26.537 

0. 

-53. 120 

0.5267 

0.3139 

0.3872 

0. 

0.2976 

1.5109 
1 .0457 
1 .0872 
1 . 0000 
0.4328 

0.2000 

0.2000 

0.2000 

0.2000 

0.2000 

5.3529 
8. 6506 
8.5684 
20.0000 
1 0.0400 


FREQUENCY= 500 HZ 

TRANSMISSION 

LOSS= 

8.69 


0. 

29.310 
-29.310 
59.477 
-59.477 
96.118 
-96. 1 18 

0. 

10.955 

-10.955 

0. 

-23.662 

0. 

-44.3.99 

0. 

23.895 
-2.1 .202 

0. 

-41.815 

0. 

0. 

0.5267 

0.3529 

0.4261 

0. 

0.3052 

0. 

0. 

1.5109 
1.1762 
1 .2054 
1 .0000 
0.7204 
1 . 0000 
1 .0000 

0.1429 
0.1429 
0.1429 
0. 1429 
0.1429 
0.1429 
0.1429 

5.3529 
7.9020 
7. 7004 
20.0000 
10.4325 
20.0000 
20.0000 


FREQUENCY= 630 HZ 

TRANSMISSION 

L0SS= 

9.74 



0. 

0. 

0. 

0,5267 

1.51 09 

0. 1 

1 1 1 

5. 3529 

23.056 

8.554 

1 8.450 

0.3954 

1 .2860 

0. 1 

1 1 1 

7. 1 333 

-23.056 

-8. 554 

-16.81 1 

0.4594 

1 .2969 

0. 1 

1 1 1 

6.9859 

46.429 

1 7.875 

40.278 

0.2935 

0.91 17 

0. 1 


3.9339 

-46.429 

-17.875 

-32.994 

0.3463 

0.9371 

0.1 

1 1 1 

9.5133 

71 .205 

0. 

0. 

0. 

1 . 0000 

0. 1 

1 1 1 

20.0000 

-71 .205 

-29.421 

-49.603 

0.2878 

0.5184 

0. 1 

1 1 1 

1 0.5856 


Figure 38. Typical Output (Continued). 



105.152 

-105.152 


0 . 

1 8.434 
-I 8.434 
37.004 
-J7.n04 
56.048 
-56.048 
76.617 
-76.617 
105.152 
-105. 152 


0 . 

1 4.742 
-1 4.742 
29. 546 
-29.546 
44.532 
-44.532 
59.979 
-59.979 
76.617 
-76.617 
97.357 
-97.357 


0 . 

11.515 
-I 1 .515 
23.056 
-23.056 
34.668 
-34.668 
46.429 
-46.429 
58.498 
-58.490 
71 .205 
-71 .205 
85.390 
-85.390 
105.152 
105.152 


0 . 

9.21 I 
-9.2 1 I 


K)(h 


0 . 0 . 0 . 

-51 .205 0. 0. 


1.0000 0. I 1 1 I 20.0000 

1.0000 0.1111 20.0000 


FREQUENCY= 800 HZ TRANSMISSION L()SS= 9.36 


0. 

0. 

6.810 

14,566 

-6.810 

-13.529 

1 3.991 

30.951 

-13.991 

-26.537 

0. 

0. 

-22.086 

-39.510 

0. 

0. 

-32.310 

-53. 120 

0. 

0. 

-51.205 

0. 


FREQUENCY=tOOO HZ 


0.5267 

1.5109 

0.4300 

1 .3624 

0.4830 

1.3598 

0.3139 

1 .0457 

0.3872 

1 .0872 

0. 

1 .0000 

0.3142 • 

0.7784 

0. 

1 .0000 

0.2976 

0.4328 

0. 

1 .0000 

0. 

1 .0000 


TRANSMISSION 


0.0909 

5.3529 

0.0909 

6.5610 

0.0909 

6.4870 

0.0909 

8.6506 

0.0909 

8.5684 

0.0909 

20.0000 

0.0909 

10.2468 

0.0909 

20.0000 

0.0909 

1 0.0400 

0.0909 

20.0000 

0.0909 

20.0000 


= 9.46 


0. 

0. 

5.432 

1 1.539 

-5.432 

-1 0.879 

1 1 .046 

24. 105 

-1 1 .046 

-21 .366 

1 7.076 

38.329 

-17.076 

-3 1 . 700 

0. 

0, 

-23.897 

-42. 152 

0. 

0. 

-32.310 

-53. 120 

0. 

0. 

-45.279 

0. 


0.5267 

1.5109 

0.4574 

1 .41 60 

0.5000 

1 .4055 

0.3515 

1.1721 

0. 4249 

1 .2019 

0.2946 

0.9351 

0.3538 

0.9677 

0. 

1 .0000 

0.3040 

0.71 18 

0. 

1 .0000 

0.2976 

• 0.4328 

0. 

1 . nooo 

0. 

1 .0000 


0.0769 5.3529 
0.0769 6.1467 
0.0769 6.1268 
0.0769 7.9292 
0.0769 7.7276 
0.0769 8.9583 
0.0769 9.3381 
0.0769 20.0000 
0.0769 10.4557 
0.0769 20.0000 
0.0769 10.0400 
0,0769 20.0000 
0.0769 20.0000 


»-REQUEMCY=1250 HZ TRANSMISSION L05S= 


0. 

0. 

4.235 

8.942 

-4.235 

-8.542 

8.554 

1 8.450 

-8.554 

-16.811 

13.058 

28.762 

-13.053 

-24.923 

1 7.875 

40.273 

-1 7.875 

-32.994 

0. 

0. 

-23.208 

-41.158 

0. 

0. 

-29.421 

-49.603 

0. 

0. 

-37.382 

0. 

0. 

0. 

-51.205 

0. 


0.5267 

1.5109 

0.4795 

1 .4548 

0.5125 

1 .4409 

0.3954 

1 .2860 

0.4594 

1 .2969 

0.3239 

1 .0841 

0.3986 

1 . 1236 

0.2935 

0.911 7 

0.3463 

0.9371 

0. 

1 .0000 

0.3076 

0.7370 

0. 

1 .0000 

0.2878 

0.5184 

0. 

1 .0000 

0. 

1 . 0000 

0. 

1 .0000 

0. 

1 . 0000 


0.0588 

5.3529 

0.0588 

5. 8393 

0.0583 

5.8530 

0.0588 

7. 1333 

0.0538 

6.9859 

0.0583 

8. 4587 

0.0538 

8.3105 

0.0583 

8.9339 

0.0583 

9.5133 

0.0583 

20. noon 

0.0538 

10.3843 

0.0583 

20.0000 

0.0588 

10.5856 

0.0588 

20.0000 

0.0588 

20.0000 

0.0538 

20.0000 

0.0588 

20.0000 


F‘RE0UENCY=1600 HZ TRANSMISSION LOSS= 9.35 


0 . 

3.384 

-3.384 


0 . 

7. I 15 
-6.859 


0.5267 

0.4934 

0.5196 


1 .5109 
1 .4769 
1.4627 


0.0476 

0.0476 

0.0476 


5.3529 
5.661 3 
5.6873 


Figure 38. Typical Output (Continued). 


105 


18.434 

6.810 

14.566 

0.4300 

1 .3624 

0.0476 

6 

-1 3.434 

-6. 810 

-13.529 

0.4830 

1.3593 

0.0476 

6 

27.689 

10.327 

22.461 

0.3631 

1 .2048 

0.0476 

7 

-27.689 

-10.327 

-20.069 

0.4347 

1.2296 

0.0476 

7 

37.004 

13.991 

30.951 

0.3139 

1 .0457 

0.0476 

8 

-37.004 

-1 3.991 

-26.537 

0.3872 

1 .0872 

0.0476 

8 

46.429 

1 7.875 

40.278 

0.2935 

0.911 7 

0.0476 

8 

-46.429 

-1 7.875 

-32.994 

0.3463 

0.9371 

0.0476 

9 

56.048 

0. 

0. 

0. 

I .0000 

0.0476 

20 

-56.043 

-22.086 

-39.510 

0.3142 

0.7784 

0.0476 

10 

66.01 1 

0. 

0. 

0. 

1 .0000 

0.0476 

20 

-66.01 1 

-26.793 

-46. 176 

0.2923 

0.6081 

0.0476 

10 

76.617 

0. 

0. 

0. 

1 .0000 

0.0476 

20 

-76.61 7 

-32.310 

-53. 120 

0.2976 

0.4328 

0.0476 

10 

88.591 

0. 

0. 

0. 

1 .0000 

0.0476 

20 

-88.591 

-39.371 

0. 

0. 

1 .0000 

0.0476 

20 

105.152 

0. 

0. 

0. 

1 . 0000 

0.0476 

20 

-105.152 

-51.205 

0. 

0. 

1 .0000 

0.0476 

20 


FREQUENCY=2000 HZ 

TRANSMISSION 

LOSS= 9.26 


0. 

0. 

0. 

0.5267 

1.5109 

0.0400 

5 

7.368 

2.705 

5.668 

0.5032 

1 .4908 

0.0400 

5 

-7.368 

-2.705 

-5.504 

0.5239 

1 .4773 

0.0400 

5 

1 4.742 

5.432 

11.539 

0.4574 

1 .4160 

0.0400 

6 

-14.742 

-5.432 

-10.879 

0.5000 

1 .4055 

0.0400 , 

6 

22. 131 

8.203 

1 7. 664 

0.4022 

1.3019 

0.0400 

7 

-22. 131 

-8.203 

-16. 157 

0.4642 

1 .3099 

0.0400 

6 

29.546 

1 1 .046 

24. 105 

0.3515 

1.1721 

0.0400 

7 

-29.546 

-1 1 .046 

-21 .366 

0.4249 

1 .2019 

0.0400 

7 

37.004 

1 3.991 

30.951 

0.3139 

1 .0457 

0.0400 

8 

-37.004 

-13.991 

-26.537 

0.3872 

1 .0872 

0.0400 

3, 

44.532 

I 7.076 

38.329 

0.2946 

0.9351 

0.0400 

8, 

-44.532 

-1 7.076 

-3 1 . 700 

0.3538 

0.9677 

0.0400 

9. 

52.169 

20.354 

46.443 

0.3079 

0.8687 

0.0400 

8- 

-52.169 

-20.354 

-36.891 

0.3259 

0.8430 

0.0400 

9, 

59.979 

0. 

0. 

0. 

I . 0000 

0.0400 

2 ^. 

-59.979 

-23.897 

-42. 152 

0.3040 

0.7118 

0.0400 

10, 

68.067 

0. 

0. 

0. 

1 . 0000 

0.0400 

20, 

-63.067 

-27.318 

-47.537 

0.2897 

0.5725 

0.0400 

10, 

76.617 

0. 

0. 

0. 

i .0000 

0.0400 

20, 

-76.61 7 

-32.310 

-53. 1 20 

0.2976 

0.4328 

0.0400 

10, 

86.01 6 

0. 

0. 

0. 

1 . 0000 

0.0400 

20, 

-86.016 

-37.765 

0. 

0. 

1 .0000 

0.0400 

20, 

97.357 

0. 

0. 

0. 

1 .0000 

0.0400 

20, 

-97.357 

-45.279 

0. 

0. 

1 .0000 

0.0400 

20, 


EREQUEMCY=2500 HZ 

TRANSMISSION 

LOSS= 8.96 


0. 

0. 

0. 

0.5267 

I .5109 

0.0303 

5. 

5.847 

2.146 

4.483 

0.5101 

1 .4997 

0.0303 

5. 

-5.847 

-2.146 

-4.380 

0.5265 

1 .4084 

0.0303 

5. 

1 1 .698 

4.302 

9.088 

0.4783 

1 .4529 

0.0303 

5. 

-1 1 .698 

-4.302 

-8.675 

0.5119 

1 .4390 

0.0303 

5. 

1 7.555 

6.481 

1 3.839 

0.4367 

1 .3759 

0.0303 

6. 

—1 7.555 

-6.481 

-12.900 

0.4873 

1.3711 

0.0303 

6. 

23.424 

8.694 

1 8.763 

0.3927 

1 .2797 

0.0303 

7. 

-23.424 

-3.694 

-1 7.071 

0.4574 

1 .2917 

0.0303 

7. 

29.310 

10.955 

23.895 

0.3529 

1.1762 

0.0303 

7, 

-29.310 

-10.955 

-2 1 . 202 

0.4261 

1 .2054 

0.0303 

7. 

35.223 

1 3.279 

29.279 

0.3214 

1 .0748 

0.0303 

8. 

-35.223 

-13.279 

-25.308 

0.3958 

1.1150 

0.O3O3 

8. 


Figure 38. Typical Output (Continued). 


.5610 
.4870 
. 7103 
.5136 
. 65^6 
.5684 
.9339 
. 5 1 33 
.0000 
.2468 
.0000 
.6358 
.^000 
.0400 
.0000 
. 0000 
.0000 
.0000 


.3529 
.5462 
'.5756 
. 1467 
. 1268 
.0165 
. 8830 
.9292 
. 7276 
6506 
5684 
9583 
3331 
3731 
9852 
0000 
4557 
0000 
6455 
0000 
0400 
0000 
0000 
0000 
0000 


3529 
4703 
4990 
8555 
8673 
4575 
3975 
1 797 
0270 
9020 
7004 
5076 
3725 
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41.176 

15.685 

34.974 

0.3005 

0.9815 

0.0303 

8. 8893 

-41 . 176 

-15.685 

-29.404 

0.3680 

1 .0214 

0.0303 

9.0083 

47. 1 84 

18.196 

41.065 

0.2937 

0.9032 

0.0303 

8.9084 

-47.184 

-1 8. 196 

-33.508 

0.3435 

0.9248 

0.0303 

9.5805 

53.213 

20.842 

47.680 

0.3173 

0.8729 

0.0303 

8.0982 

-53.213 

-20.842 

-37.637 

0.3224 

0.8247 

0.0303 

10.0647 

59. All 

0. 

0. 

0. 

1 . 0000 

0.0303 

20.0000 

-59. All 

-23.662 

-41.815 

0.3052 

0.7204 

0.0303 

10.4325 

65.849 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

-65.849 

-26.714 

-46.068 

0.2926 

0.6109 

0.0303 

10.6337 

72.470 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

-72.470 

-30.082 

-50.431 

0.2880 

0.4970 

0.0303 

10.5258 

79.481 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

-79.481 

-33.909 

-5A.950 

0.3527 

0.4328 

0.0303 

8.3979 

87.146 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

-87.146 

-38.463 

0. 

0. 

1 .0000 

0.0303 

20.0000 

96.118 

0. 

0. 

0. 

I .0000- 

0.0303 

20.0000 

-96. 118 

-44.399 

0. 

0. 

I .0000 

0.0303 

20.0000 

109.677 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

-109.6-77 

-54.998 

0. 

•0. 

1 .0000 

0.0303 

20.0000 


FREQUENCY=3150 HZ 

TRANSMISSION LOSS= 9.06 


0. 

0. 

0. 

0.5267 

1 .5109 

0.0244 

5.3529 

4.605 

1 .689 

3.521 

0.5151 

1 .5051 

0.0244 

5.4228 

-4.605 

-1.689 

-3.457 

0.5279 

1 .4957 

0.0244 

5.4475 

9.2J 1 

3.384 

7. 1 15 

0.4934 

1 .4769 

0.0244 

5.661 3 

-9.21 1 

-3.384 

-6.859 

0.5196 

1 .4627 

0.0244 

5.6873 

13.820 

5.089 

10.792 

0.4639 

1 .4279 

0.0244 

6.0529 

-1 3.820 

-5.089 

-10.21 3 

0.5038 

1 .4161 

0.0244 

6.0441 

18.434 

6.810 

1 4.566 

0.4300 

1 .3624 

0.0244 

6.5610 

-1 8.434 

-6.810 

-13.529 

0.4830 

1.3598 

0.0244 

6.4870 

23.056 

8.554 

1 8.450 

0.3954 

1 .2860 

0.0244 

7. 1 333 

-23.056 

-8.554 

-1 6.81 1 

0.4594 

1 .2969 

0.0244 

6.9859 

27.689 

10.327 

22.461 

0.3631 

1 .2048 

0.0244 

7.7103 

-27.689 

-10.327 

-20.069 

0.4347 

1 . 2296 

0.0244 

7.5136 

32.337 

12. 137 

26.619 

0. 3355 

1 .1236 

0.0244 

8.2338 

-32.337 

-12. 137 

-23.308 

0.4103 

1.1595 

0.0244 

8.0474 

37.004 

13.991 

30.951 

0.3139 

1 .0457 

0.0244 

8.6506 

-37.004 

-13.991 

-26.537 

0.3872 

1 .0872 

0.0244 

8.5684 

41 .699 

15.900 

35.489 

0.2993 

0.9739 

0.0244 

S.9082 

-41 .699 

-15.900 

-29.763 

0.3657 

1.0131 

0.0244 

9.061 3 

46.429 

1 7.875 

40.278 

0.2935 

0.91 1 7 

0.0244 

8.9339 

-46.429 

-1 7.875 

-32.994 

0.3463 

0.9371 

0.0244 

9.5133 

51.207 

1 9.931 

45.378 

0.3026 

0.8706 

0.0244 

8.5509 

-51.207 

-19.93J 

-36.239 

0.3291 

0.8590 

0.0244 

9.9127 

56.048 

0. 

0. 

0. 

1 .0000 

0.0244 

20.0000 

-56.048 

-22.086 

-39.510 

0.3142 

0.7784 

0.0244 

1 0.2468 

60.972 

0. 

0. 

0. 

1 .0000 

0.0244 

20.0000 

-60.972 

-24. J63 

-42.817 

0.301 / 

0.6949 

0.0244 

10.4982 

66.01 1 

0. 

0. 

0. 

1 .0000 

0.0244 

20.0000 

-66.01 1 

-26.193 

-46.176 

0.2923 

0.6081 

0.0244 

10.6358 

71.205 

0. 

0. 

0. 

1 . 0000 

0.0244 

20.0000 

-71.205 

-29.421 

-49.603 

0.2878 

0.5184 

0.0244 

10.5855 

76.617 

0. 

0. 

0. 

1 .0000 

0.0244 

20.0000 

-76.617 

-32.310 

-53. 120 

0.2976 

0.4328 

0.0244 

10.0400 

82.351 

0. 

0. 

0. 

1 . 0000 

0.0244 

20.0000 

-82.351 

-35.565 

0. 

0. 

1 .0000 

0.0244 

20.0000 

88.591 

0. 

0. 

0. 

1 . 0000 

0.0244 

20.0000 

-88.591 

-39.371 

0. 

0. 

1 .0000 

0.0244 

20.0000 

95.741 

0. 

0. 

0. 

1 .0000 

0.0244 

20.0000 

-95.741 

-44. 134 

0. 

0. 

1 .0000 

0.0244 

20.0000 

105.152 

0. 

0. 

0. 

1 .0000 

0.0244 

20.0000 

-105.152 

-51 .205 

0. 

0. 

1 . 0000 

0.0244 

20.O0O0 


EREQUENCY= 

4000 HZ 

TRANSMISSION L()SS= 

9.15 



Figure 38. Typical Output (Continued). 
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3-STG LPT* 3.n PR 


100% N 


3 STAGES 

* AERO-THERMO PARAMETERS * 



STATION 

U- FPS 

PHI- DEG 

VR- FPS 

PS- PS I A 

TS- 

1 

1 

673.000 

0. 

0. 

33.480 

723 

1 

2 

523.000 

62.100 

409.000 

23.410 

657 

2 

3 

413.000 

-42.900 

0. 

21 .060 

641 

2 

4 

407.000 

62.000 

423.000 

16.910 

608 

3 

5 

349.000 

-26.500 

0 . 

15.080 

591 

3 

6 

31 7.000 

55.100 

437.000 

1 3.-770 

579 


7 

322.000 

1 8.260 

0. 

12.830 

569 



***** FIRST CUT-ON OCCURS AT 193. HZ 

***** 


THETA- I 

THETA -R THETA-T T B 

E 

T 

0. 

G. 0. 0.4913 1.6899 

1 . 0000 

4 


FREQUENCY= 50 HZ TRANSMISSION 

L()SS= 4.08 



0. 

0 . 0 . 

0.4913 1 

.6899 

1 . OOOO 4 


FREQUENCY= 63 

HZ TRANSMI 

SSI ON 

LOSS= 4.03 

0. 

C 

• 

o 

0.4913 1 

.6899 

1 . 0000 4 


FREQUEMCY= 80 HZ TRANSMISSION LOSS= 4.08 

G- 0. 0. 0.4913 1.6899 1.0000 4 

EREQUENCY= 100 HZ TRANSMISSION LOSS= 4.08 

G. 0. 0.4913 1.6899 1.0000 4 

EREQUENCY= 125 HZ TRANSMISSION L0SS= 4,08 

G. 0. 0.4913 1.6899 1.0000 4 

FREQUENCY= 160 HZ TRANSMISSION L0SS= 4.08 

G- G. 0. 0.4913 1.6899 0.3333 4 

Figure 38. Typical Output (Continued). 


DEG R 
000 
000 
000 
000 
000 
000 
ono 


-LOSS 

.0764 


.0764 


. 0764 


0764 


0764 


0764 


0764 
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OtiKiCNAL 

i,v yi.m ' 


101 .591 
-101.591 

1 . 689 
-42.721 

3.521 

-3.457 

0. 

0. 

1 .0000 
1 .0000 

0. 3333 
0.3333 

20.0000 

20.0000 


FREQUENCY^ 200 HZ 

TRANSMISSION 

L0S5= 

8.63 


0. 

77.135 

-77.135 

0. 

1.689 

-28.544 

0. 

3.521 

-3.457 

0.491 3 

0. 

0. 

1 .6899 
1 .0000 
1 .0000 

0.3333 

0.3333 

0.3333 

4.0764 

20.0000 

20.0000 


FREUUENCY= 250 HZ 

TRANSMISSION 

L0S5= 

8.63 


0. 

60.501 

-60.501 

0. 

1 .609 
-2 1 . 080 

0. 

3.521 
-44 . 766 

0.491 3 

0. 

0.2306 

1.6899 
1 .0000 
0.8926 

0. 3333 
0.3333 
0.3333 

4.0764 
20. OOOO 
10.8433 


FREQUENCY= 315 HZ 

TRANSMISSION 

LOSS= 

7.93 


0. 

47.488 
-47.4B8 
101 .5V1 
-101 .591 

0. 

15.977 
-15.977 
3.384 
-42. 721 

0. 

38.565 

-35.095 

7.115 

-6.859 

0.4913 

0.2034 

0.2575 

0. 

0. 

1 .6899 
1.1035 
1.1111 
1 .0000 
1 .0000 

0.2000 

0.2000 

0.2000 

0.2000 

0.2000 

4.0764 
10.2793 
10.044 / 
20.0000 
20. OOOO 


FREQUENCY= 400 HZ 

TRANS 

MISSION 

L0SS= 

9. I B 


0. 

37.962 

-37.962 

77.135 

-77.135 

0. 

12.524 
-12.524 
3.384 
-2 8. 544 

0. 

30.220 

-28.073 

7.115 

-6.859 

0.4913 

0.2007 

0.2932 

0. 

0. 

1 .6899 
1 .2347 
1 .2554 
1 .0000 
1 .oooo 

0.20O0 

■0.2000 

0.2000 

0.2000 

0.2000 

4.0764 
1 0. 4037 
8. 9554 

20. 0000 

20.0000 


FREQUENCY= 500 HZ 

TRANSMISSION 

LOSS= 

9.00 


0. 

30. 127 
-30. 127 
60.501 
-60.501 
94.733 
-94.733 

0. 

9.815 
-9.815 
3.384 
-21 .080 
5.089 
-38.226 

0. 

23.642 
-22.316 
7.115 
-44.766 
1 0.792 
-1 0.21 3 

0.4913 

0.2371 

0.3344 

0. 

0.2306 

0. 

0. 

1 .6899 • 
1 .3525 
1 .3659 
1 .0000 
0.8926 
1 .0000 
1 . 0000 

0. 1 429 
0. 1 429 
0. 1429 
0. 1 429 
0. 1 429 
0. 1 429 
0. 1 429 

4.0764 
9.5699 
7.8009 
20.0000 
1 0.8433 
20.0000 
20.0000 


FRHQUEMCY= 630 riZ TF?ANSMISSI()N L()SS= 9.53 


0. 

0. 

0. 

0.4913 

1 .6899 

0. 

1 1 1 

4.0764 

23.728 

7.668 

1 8.428 

0.2796 

1 .4558 

0. 

1 1 1 

8. 345 7 

-23. 72 8 

-7.668 

1 7.616 

0.3759 

1.4517 

0. 

1 1 1 

6. 7459 

47.488 

15.977 

38.565 

0.2034 

1 .1035 

0. 

1 1 1 

1 O. 2793 

-47.488 

-15.977 

35.095 

0.2575 

1 . 1 1 1 1 

0. 


1 0. 0447 

71 .956 

5.089 

10.792 

0. 

1 . OOOO 

0. 

! 1 } 

20.0000 

-71 .956 

-26.080 

53.357 

0.2378 

0.6830 

0. 

111 

10.3005 

101 .591 

6.310 

14.566 

0. 

1 .OOOO 

0. 

1 1 1 

20. OOOO 

-1 01 .591 

-42.721 

13.529 

o. 

1 .0000 

0. 

1 1 1 

20.0000 


frequency 

= 800 HZ 

TRANSMISSION 

L0SS= 

^.?7 



Figure 38. Typical Output (Continued). 
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0. 

0. 

0. 

0.4913 

1 .6899 

0.0909 

4.076 4 

1 8.984 

6.107 

1 4.642 

0.3239 

1 .5338 

0.0909 

7.2 1 29 

-18.984 

-6.107 

-1 4. 1 26 

0.4098 

1 .5127 

0.0909 

5, 9530 

37.962 

12.524 

30.220 

0.2087 

1 .2347 

0.0909 

1 0.4037 

-37.962 

-12.524 

-28.073 

0.2932 

1 .2554 

0.0909 

8 . 955 4 

57. 105 

5.089 

10.792 

0. 

1 .0000 

0.0909 

20.0000 

-57.105 

-19.696 

-42.231 

0.2353 

0.9521 

0.0909 

10.7237 

77.135 

6.810 

14.566 

0. 

1 .0000 

0.0909 

20.0000 

-77.135 

-28.544 

-13.529 

0. 

1 . 0000 

0.0909 

20,0000 

101 .591 

8. 554 

1 8.450 

0. 

1 .0000 

0.0909 

20, nooo 

-101 .591 

-42.721 

-16.81 1 

0. 

1 . 0000 

0.0909 

20,0000 


FREQUENCY=1000 HZ 

TRANSMISSION 

L()SS= 9.67 


0. 

0. 

0. 

0.491 3 

1 .6899 

0.0769 

4.0764 

15.189 

4.871 

11.655 

0.3655 

1 .5923 

0.0769 

6.2727 

-15.189 

-4.871 

-1 I .326 

0.4369 

1 .5593 

0.0769 

5 . 35 1 5 

30.368 

9.897 

23.841 

0.2358 

1 .3488 

0.0769 

9.6069 

-30.368 

-9.897 

-22.493 

0.3330 

1 .3626 

0.0769 

7.8392 

45.578 

15.268 

36.857 

0.201 6 

1 . 1283 

0.0769 

10.4281 

-45 . 578 

-15.268 

-33.633 . 

0.2635 

1.1410 

0.0769 

9.8574 

60.998 

6.810 

14.566 

0. 

1 . 0000 

0.0769 

c 

o 

c 

• 

c 

-60.998 

-21.285 

-45. 1 37 

0.2300 

0.8838 

0.0769 

10.8544 

77.135 

8.554 

1 8.450 

0. 

1 .0000 

0.0769 

20. OOOO 

-77.135 

-20.544 

-16.811 

0. 

1 . 0000 

0.0769 

20.0000 

95.753 

1 0.327 

22.461 

0. 

1 . ooon 

0.0769 

20. 0000 

-95.753 

-38.804 

-20.069 

0. 

1 .0000 

0.0769 

20. 0000 


ERE0UEMCY=I250 HZ 

i'RAi' 

IS. MI SSI ON 

L0SS= 9.32 



0, 

0. 

0 . 

0.491 3 

1 .6899 

0.0588 

4.0764 

11.867 

3. 798 

9.068 

0.4033 

1.6361 

0.0588 

5.5097 

-1 1 .867 

-3.798 

-8.867 

0.4589 

1 .5976 

0.0588 

4.8796 

23.728 

7.663 

1 8.428 

0.2796 

1 .4553 

0.0588 

8.345 7 

-23.728 

-7.668 

-17.616 

0.3759 

I .4517 

0.0588 

6. 7459 

35.588 

1 1 . 692 

20.202 

0.2149 

1 .2695 

8.0588 

1 0.2328 

-35.588 

-11.692 

-26.327 

0.3045 

1 .2896 

0.0583 

8.6271 

4 / . 4 80 

15.977 

38, 565 

0.2034 

1.1035 

0.0588 

10.2793 

-47.488 

-15.977 

-35.096 

0,2575 

1.1111 

0.0588 

10.0447 

59.534 

8. 554 

1 8.450 

0. 

1 .OOOO 

0.0588 

20.0000 

-59.534 

-20.681 

-44.043 

0.231 7 

0.9098 

0.0588 

10.3167 

71.956 

10,327 

22.461 

0. 

1 .0000 

0.0583 

20.0000 

-71 .956 

-26.080 

-53.357 

0.2378 

0.6830 

0.0588 

10.3005 

85.318 

12.137 

26.619 

0. 

1 .0000 

0.0588 

20.0000 

-85.31 8 

-32.763 

-23.308 

■T. 

1 .0000 

0.0588 

20. OOOO 

101.591 

1 3 . 99 1 

30.951 

0. 

1 .0000 

0.0583 

20.0000 

-101 .591 

-42.721 

-26,537 

0. 

1 .0000 

0.0588 

20.0000 


FRE0UEMCY=1600 HZ 

TRANSMISSION 

LOSS= 9.26 


0. 

0. 

0. 

0.4913 

1 .6899 

0.04/6 

4 . 0764 

9.494 

3.035 

7.234 

0.4289 

1 .6610 

0.0476 

5.0376 

-9.494 

-3.035 

-7.106 

0.4725 

1 .6229 

0.0476 

4.5900 

1 8.984 

6. 107 

1 4.642 

0. 3239 

1 .5338 

0.0476 

7.2129 

-1 8.984 

-6. 107 

-14. 126 

0.4098 

1 .5127 

0.0476 

5. 9530 

28.470 

9.254 

22.280 

0.2462 

1 .3787 

0.0476 

9.2965 

-28.470 

-9.254 

-21.100 

0.3445 

I .3885 

0.0476 

7.5343 


Figure 38. Typical Output (Continued). 
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I 



37.962 
37.962 
47.488 
-47.488 
■yi.\ob 
-57.105 
66.91 8. 
- 6() . 9 I 8 
77.135 
-77. 1 35 
88.214 
-88.21 4 
101 .591 
101 .591 


0 . 

7.596 
-7.d96 
15. 1 89 
-15.189 
22.779 
-22. 779 
30.368 
-30.363 
37.962 
-37.962 
15.578 
-45.573 
53.242 
-53.242 
60.998 
-60.998 
68.919 
-68.91 9 
77. 1 35 
-77.135 
85.888 
-85.808 
95.753 
-95.753 
109.031 
109.031 


0 . 

6.029 
-6.029 
12.056 
-12.056 
18.081 
-1 8.081 
24.104 
-24. 104 
30.127 
-30.127 
36.153 
-36. 153 
42.190 
-42.190 


UHAUilVAi. 8 

POOR QU/j 


12.524 

30. 220 

0.2087 

1 .2347 

0.0476 

1 0. ^'^37 

-12.524 

-28.073 

0.2932 

1 .2554 

0.0476 

8.9554 

1 5.977 

38.565 

0.2034 

1 .1035 

0.04/6 

10.2793 

-15.977 

-35.095 

0.2575 

1.1111 

0.0476 

10.0447 

10.327 

22.461 

0. 

1 .0000 

0.0476 

20. OOOO 

- 1 9. 696 

-42.231 

0.2353 

0.9521 

0.0476 

1 0. 7237 

12. 137 

26.619 

0. 

1 .0000 

0.0476 

20. oooo 

-23. 81 I 

-49.571 

0.2276 

0. 7760 

0.0476 

10.81 9-) 

1 3 . 9v 1 

30.951 

0. 

1 .0000 

0.0476 

20. OOOO 

-28.544 

-26.537 

0. 

1 . 0000 

0.0476 

2^.0(100 

1 5 . 900 

35.489 

0. 

1 .0000 

0.0476 

20. OOO') 

-34.367 

-29.763 

0. 

1 . noOO 

0.0476 

20.0000 

1 7.875 

40.278 

0. 

1 .0000 

0.0476 

20. OOOO 

-42. 721 

-32.994 

0. 

1 .rtoon 

0.0476 

20.0000 

1-REQU2 

NCY=2000 HZ 

'HUMSMISSION 

LOSS= 9.33 


0. 

0. 

0.4913 

I .689V 

0.0370 

4.0764 

2.426 

5. 775 

0.4473 

1 .6761 ■ 

.0.037^ 

4.7190 

-2.426 

-5.693 

0.4815 

1 .6413 

0.037') 

4.39:>5 

4.371 

11.655 

0 . 3655 

1.5923 

0.0370 

6.2121 

-4.371 

-1 1 .326 

0.4369 

1 .5593 

0.0370 

5.3515 

7.354 

1 7. 666 

0.2876 

1 .4715 

0.0370 

8. 1 300 

-7.354 

-16.913 

0.3826 

1 .4641 

0.03/0 

6.5862 

9.897 

23.341 

0.2358 

1 .3488 

0.0370 

9.6069 

-9.897 

-22.493 

0.333'’' 

1 . 3626 

0.03 JO 

7.8392 

12.524 

30.220 

0.2087 

1 .2j47 

0.0370 

10.408/ 

-12.524 

-20.073 

0.2932 

1 .2554 

0.0370 

8. 9554 

1 5.268 

36.857 

0.2016 

1 . 1 283 

0.0370 

10.4281 

-1 5.26'3 

-33.683 

0.2635 

1- . 1 4 1 0 

0.0370 

9.8571 

18. 1 70 

43.829 

0.2364 

1 .1^594 

0.0370 

3. 7630 

-18.1 70 

-39.358 

0.2427 

1 .01 77 

. 0.0370 

10.506? 

1 3 . 99 1 

30.951 

0, 

1 .0000 

0.037'^ 

20.0000 

-21 .285 

-45. 1 37 

O.23on 

n. 883.3 

0.0370 

10.8544 

15.900 

35.489 

9. 

1 . 0000 

0.0370 

20.O00O 

-24.698 

-5 1 .075 

0.229.3 

0.7388 

0.03 /o 

10.7057 


17.8/5 40.2 78 0. 

28.544 -32.994 0. 

19.931 45.378 0. 

33.074 -36.239 0. 

0 . 0 . 0 . 

•38.864 -39.510 0. 

0 . 0 . 0 . 

■48.219 -42.017 0. 

8R80UHMCY=2500 HZ 


1.0000 0.0370 
1 .^'000 0.0370 

1.0000 0.0370 

1.0000 0.03/0 

1.0000 O.0370 

1.0000 0.0370 
I .OQOO 0.0370 

1.0000 0.0370 

IRANS.'HSSIOM LOSS= 9.31 


0 . 

1 .925 
-1.925 
3.859 
-3.859 
5.812 
-5.81 2 
7.793 
-7.793 
9.81 5 
-9.315 
1 I . R89 
-1 I .889 
14.032 
-14.032 


0 . 

4.575 
-4.524 
9.214 
-9.007 
I 3.928 
-1 3.460 
1 8.731 
-I 7.893 
23.642 
-22.316 
28.680 
-26.742 
33.871 
-31.1 34 


0.4913 
0.4607 
0.4873 
0.4012 
0.4573 
0.3334 
0.4163 
0.2765 
0.3733 
0.2371 
0. 3344 
0.21 33 
0.301 7 
0.2022 
0.2755 


1 .6899 
1 .6051 
I .6550 
1 .6339 
I .5955 
I .5403 
1 .5240 
1 .4496 
I .4467 
1 .3525 
1 .3659 
1.2611 
1.2815 
1.1746 
1 . 1 929 


0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0303 
0.0 J03 
0.0303 


Z 0.0000 
20.^’000 
20. OOO'.l 
20.^000 
20 . noOO 
20.0000 
20.0000 
20. OOOO 


4.0764 

4.5012 

4.2643 

5. 5502 
4. 9045 

6. 987 5 
5. 8056 
8.4293 
6.0093 
9.5697 
7.8009 

10.2800 
8.70 71 
10.5237 
9. 4059 


Figure 38. Typical Output (Continued). 


48.247 

1 6.261 

39.249 

0.2047 

1 .0940 

0.0303 

-48.247 

-16.261 

-35.656 

0.2553 

1 .0991 

0.0303 

54.343 

15.900 

35.489 

0. 

1 . 0000 

0.0303 

— 54.343 

-18.600 

-40. 1 76 

0.2403 

0.9992 

0.0303 

60.501 

1 7.875 

40.278 

0. 

1 .0000 

0.0303 

-60.501 

-2 1 . 080 

-44.766 

0.2306 

0.8926 

0.0303 

66.760 

1 9.931 

45.378 

0. 

1 .0000 

0.0303 

—66. 760 

-23.741 

-49.453 

0.2276 

0.7789 

0.0303 

73.175 

0. 

0. 

0. 

1 .0000 

0.0303 

— 13 . 1 75 

-26.647 

-54.274 

0.2470 

0.6624 

0.0303 

79.838 

0 . 

0. 

0. 

1 . 0000 

0.0303 

-79.838 

-29.891 

-42.817 

0. 

1 .0000 

0.0303 

86.913 

0. 

0. 

0. 

1 .0000 

0.0303 

-86.913 

-33.639 

-46. 1 76 

0. 

I .0000 

0.0303 

94. 733 

0. 

0. 

0. 

1 .0000 

0.0303 

-94.733 

-38.226 

-49.603 

0. 

1 .0000 

0.0303 

104.250 

0. 

0. 

0. 

1 .0000 

0.0303 

-104.250 

-44.604 

-53. 120 

0. 

1 .0000 

0.0303 


0. 

0. 

4.748 

1.515 

-4.743 

-1.515 

9.494 

3.035 

-9.494 

-3.035 

1 4.240 

4.564 

-I 4.240 

-4.564 

1 8.984 

6. 107 

-18.984 

-6. 107 

23.728 

7. 668 

-23.728 

-7. 668 

28.470 

9.254 

-28.470 

-9.254 

33.214 

10.871 

-33.214 

-10.871 

37.962 

12.524 

-37.962 

-12.524 

42.719 

14.223 

-42.719 

-14.223 

47.488 

15.977 

-47.488 

-15.977 

52.280 

17.797 

-52.280 

-17.797 

57.105 

0. 

-57.105 

-19.696 

61 .977 

0. 

-61 .977 

-21 .693 

66.918 

0. 

-.66.91 8 

-23.81 1 

71 .956 

0. 

-71 .956 

-26.080 

77.135 

0. 

-77.135 

-28.544 

82.518 

0. 

-82.518 

-31 .271 

88.214 

0. 

-88.214 

-34.367 

94.419 

0. 

-94.419 

-38.031 

101 .591 

0. 

101 .591 

-42.721 

1 1 1 .521 

0. 

M 1 .521 

-50.228 


HREQUENCY=3I50 HZ 

0 . 0 . 

1.515 3.598 

1.515 -3.566 

3.035 7.234 

3.035 -7.106 

4.564 10.913 

4.564 -10.624 

6.107 14.642 

6.107 -14.126 

7.668 18.428 

7.668 -17.616 

9.254 22.280 

9.254 -21.100 

7.871 26.207 

3.871 -24.583 

^.524 30.220 

?.524 -28.073 

U223 34.334 

1.223 -31.574 

j.977 38.565 

■>.911 -35.095 

^797 42.936 

^797 -38.644 

). 0 . 

>.696 -42.231 

>. 0 . 

.693 -45.869 

'. 0 . 

i.811 -49.571 

I. 0. 

.080 -53.357 i 

• 0. j 

.544 -53.120 ( 

. 0 . ( 

.271 0. { 

. 0 . ( 

.367 0. ( 

. 0 . ( 

.031 0. ( 

. 0 . ( 

.721 0. f 

. 0. c 

.228 0. C 


TRANSMISSION 
0.4913 1.6899 


0.4701 
0.4907 
0.4289 
0.4725 
0.3764 
0.4435 
0.3239 
0.4098 
0.2796 
0.3759 
0.2462 
0.3445 
0.2231 
0.3168 
0.2087 
0.2932 
0.2018 
0.2735 
0.2034 
0.2575 
0.222 7 
0.2448 
0 . 

0.2353 

0 . 

0.2291 

0 . 

0.2276 

0 . 

0.2378 

0 . 


I .6901 
I .6650 
I .6610 
I .6229 
I .6057 
I .5705 
1 .5338 
1 .5127 
1 .4558 
1.4517 
1 .3787 
1 .3885 
1 .3050 
1 .3232 
I .2347 
1 .2554 
1.1673 
1 . 1 849 . 
1.1035 
I. I 1 1 1 
1 .0553 
1 .0336 
I .0000 
0.9521 
I . 0000 
0.8663 
1 .0000 
0.7760 

1 .0000 
0.6830 

1 . 0000 
1 .0000 
1 .0000 

1 .0000 

1 .0000 
I .0000 

1 . 0000 

1 .0000 

1 . 0000 

1 .0000 
1 .0000 
1 .0000 


L()SS= 9.26 

0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.0233 
0.02 33 
0.0233 
0.0233 
0.0233 
0.0233 ; 

0.0233 I 
0.0233 ; 

0.0233 1 

0.0233 ; 

0.0233 1 

0.0233 ; 

0.0233 I 
0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 

0.0233 2 


FREQUENCY=4000 HZ TRANSMISSION L0SS= 9.33 
Figure 38, Typical Output (Continued). 
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3-STG LPTs 4.0 PR 100% N 


3 STAGES 

* AERO-THERMO PARAMETERS * 


STAGE 

STATION 

U- FPS 

PHI- DEG 

VR- FPS 

PS- PSIA 

TS- DEG R 

1 

1 

681 .000 

0. 

0. 

33.340 

722.000 

1 

2 

554.000 

62.100 

409.000 

21.820 

645.000 

2 

3 

445.000 

-45.100 

0. 

1 9.260 

626.000 

2 

4 

462.000 

62.000 

428.000 

14.320 

581 .000 

3 

5 

403.000 

-33.700 

0. 

1 2.430 

562.000 

3 

6 

377.000 

55.100 

437.000 

1 1 .040 

546.000 


7 

386.000 

5.900 

0. 

10.030 

534.000 



***** FIRST 

CUT-ON OCCURS AT 

1 87 . HZ ***** 


THETA- I 

THETA-R 

THETA-T 

T 

B E 

T-LOSS 

0. 

0. 

0. 

0.4803 

1.6387 1.0000 

2.9605 


FREQUENCY= SO HZ TRANSMISSIOH L()SS= 2.96 

0. 0. 0. 0.4003 1.6387 1.0000 2.9605 

. FREQUENCY= 63 HZ TRANSMISSION L()SS= 2.96 

0. 0. 0. 0.4803 1.6387 1.0000 2.9605 

FREOUENCY= 80 HZ TRANSMISSION LOSS= 2.96 

0. 0. 0. 0.4803 1.6387 1.0000 2.9605 

FREQUENCY= 100 HZ TRANSMISSION LOSS= 2.96 

0. 0. 0. 0.4803 1.6387 1.0000 2.9605 

FREQUEMCY= 125 HZ TRANSMISSION LOSS= 2.96 

0. 0. 0. 0.4803 1.6387 1.0000 2.9605 

FREQUENCY= 160 HZ TRANSMISSION LOSS= 2.96 

0. 0. 0. 0.4803 1.6387 0.3333 2.9605 

Figure 38. Typical Output (Continued). 
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ni .257 
-lOI .257 

1 . i 5 
-41 .HO 7 

3.598 

-3.566 

0. 

0. 

1 .0000 
1 . 0000 

0..3333 

0.3333 

20.0000 

20.0000 


FRHOUE 

ilCY= 200 ilZ 

TRANSMISSIOf! 

LOSS= 

7.56 


0. 

77. 1 85 
-77. 105 

0. 

1.515 

-28.085 

0. 

3.598 

-3.566 

0.4803 

0. 

0. 

1.6637 
1 . onoo 
1 . 0000 

0.3333 
0.3333 
0. 3333 

2.9605 

20.0000 

20,0000 


t\lUuUS 

MCY= 250 HZ 

THAMS-MISSIorJ 

LOGS= 

7.56 


n. 

60.616 

-60.616 

0. 

1 .515 
-20.162 

0. 

3.598 

-46.740 

0.4803 

0. 

0. 1 900 

1 .6387 
1 .0000 
0.9404 

0. 3333 
0.3333 
0.3333 

2.9609 
20. 0000 
1 1 . 1763 


r'REOUE 

MCY= 315 HZ 

TfMNSMISSlOiJ 

L055= 

7.05 


- • 

47.6 1 1 
-47.61 1 
101.257 
-101 .257 

0. 

15.744 
-15. 7^4 
3,035 
-4 1 .807 

0. 

37.799 
-36.505 
7.234 
-7. 106 

0.4303 
0. 1593 
0.21 01 

0. 

0. 

1 .6337 
1 . 1 509 
1 . 1 540 

1 . 0000 

1 .0000 

0.2000 

0.2000 

0.2000 

0.2000 

0.2O'’>0 

2.960'_) 
1 1 . 1 5n 7 
1 0. 440o 
20.0000 
20 . onoo 


hREOUHNCY= 400 ilZ 

TRAM^i 

i.lISSIOM 

LOSS= 

6 . 5 8 


0. 

60.074 
-68.074 
77.1 85 
-77. 105 

0. 

12.344 
-12.344 
3.035 
-2 8.035 

0. 

29.914 
-29. 106 
7.234 
-7. 106 

0.4803 
0. 1599 
0.2^1 5 

0. 

0. 

1 .6387 
! .2793 
1 .2381 
1 . 0000 
1 .0000 

o. 2000 
0. 2000 
0.20'V) 
0.2000 
0.200 ) 

2. 9605 
11.4833 

9. 2664 
20. 0000 
20.nOj ) 


FRhOUL'RCY= 5oo !IZ 


3 1 35 ION 

LOSS= 

8.45 



0. 

30.223 

-30.223 

60.616 

-60.616 

94.579 

-94.578 

0. 

9.675 
-9.675 
3.035 
-20.762 
4,564 
-j7. 504 

0. 

23.570 
-23.067 
7.234 
-46. 740 
10.913 
-10.624 

0. 4;103 

0. 1 81 3 

0.2816 

0. 1 900 

0. 

0. 

1 .6387 
1.3795 
1 .3847 
1 .0000 
0.9404 

1 . 0000 

1 . 0000 

0.1429 
0.1 429 
0. 14^9 
n. 1 429 
0. 1429 
O. 1 429 
0. 1 429 

2.9609 
1 0.6591 
7.9244 
20.00'^0 
1 1 . i 733 
20.0000 
20 .no')o 


FREOUENCY= 630 HZ 

TRAKSHISSroN 

L0S5= 9.11 


0. 

23.807 
-23.807 
47.61 1 
-47.61 1 
72.036 
-72.036 
101 .257 
-101 .257 

0. 

7.559 
-7.559 
15.744 
-15.744 
4 , 564 
-25.67! 

6. 107 
-41.807 

0. 

1 8.468 
-18. 1 59 
37.799 
-36.505 
10.913 
-55.857 
1 4.642 
-14.126 

0.4003 
0.2194 
0.3250 
0. 1593 
0.2101 
0. 

0.2009 

0. 

0. 

I .6387 
I .4602 
1 .4551 
1 . 1 509 
1 . 1 540 
1 . OOOO 
0.7270 
1 .0000 
1 .0000 

0.1111 
0.1111 
0.1111 
0.1111 
0.1111 
0.1111 
0.1111 
0.1111 
0. 11 1 1 

2.9605 
9. 1 96 7 
6.6254 
11.1567 
10 . 4406 
20.0000 
10. 1 223 
20.00OO 
20.000') 


FREQUiZ 

MCY= 800 HZ 

TRANSMISSION 

L0SS= 8. 98 



Figure 38* Typical Output (Continued) 



0. 

0. 

0. 

0.4803 

1.6387 

0.09 ) 9 

2 . 96' ) ' 

1 9.0b0 

6.020 

14.727 

0.2656 

1.519' 

0.0909 

7. 671 ■ ) 

-1 9.050 

-6.020 

-14.530 

0. ^649 

1 .5027 

0.0909 

5 . 605 1 

38.074 

1 2.344 

29.914 

0. 1599 

1 .2793 

0.0909 

11.483^ 

-38.074 

-12.344 

-29. 106 

0.2415 

1 .2881 

0.0909 

9.2661 

■d1.2?A 

4.564 

10.913 

0. 

1 . 0000 

0.0909 

20.0000 

-bl.22A 

-19. 403 

-44.054 

0. 1929 

0.999J 

0.0909 

1 I .'^91 9 

77. I8'j 

6. 107 

1 4.642 

0. 

1 .0000 

0.0909 

r 

C 

c 

• 

c 

-77.18b 

-20.005 

-14. 1 26 

0. 

1 .0000 

0.0909 

20.000'', 

101 .?57 

7.660 

1 3. 428 

0. 

1 .0000 

0.0909 

20.0000 

-101.25/ 

-41 .307 

-17.616 


1 .0000 

0.0909 

20.0003 


FHFOUn 

MCY=IOOO HZ 

TRANS. 

USS I ON 

L05S= 9.38 


o. 

0. 

0. 

0.4803 

1 .6387 

0.0769 

2 . 9605 

1 b.242 

4.802 

1 1 . 754 

0.3145 

i .5647 

0.0769 

6. 304 1 

-I 5.242 

-4.802 

-1 1 .629 

0.3985 

1 .5331 

0.0769 

4. 8039 

30.464 

9. 755 

23.763 

0. 1002 

1.3765 

0.0769 

1 0.6997 

-30.464 

-9. 755 

-23.252 

0.2801 

1 .3819 

0.0769 

7.970 3 

45.699 

1 5.046 

36.201 

0. 1570 

1 . 1 771 

0.0769 

1 1 . 35 '>3 

-45.699 

-15.046 

-35.016 

0.2151 

1 . 1 824 

0.0769 

10.246/ 

61.112 

6. 107 

1 4.642 

0. 

1 . 0000 

0.0769 

20.0000 

-61.112 

-20.964 

-47. 1 33 

^.1896 

0.9316 

0.0769 

11.1 855 

77.105 

7. 668 

1 8.428 

1 ^ 

1 . 0000 

0.0769 

20.0000 

-77.185 

-28.085 

-17.616 

0. 

1 .0000 

0.0769 

20.3000 

95.577 

9.254 

>2.280 

0. 

1 .0000 

0.0769 

20.0000 

-95.577 

-30. 1 1 9 

-21.100 

0. 

1 .oooo 

0.0769 

20.0003 


FREOUli 

RCY=1250 HZ 

TRANSMISSION 

LOSS= 9.07 


0. 

0. 

0. 

0.4803 

1 .6387 

0.0588 

2.9605 

1 1 .909 

3.745 

9. 1 66 

0.3628 

1 .5989 

0.0583 

5. 1 452 

-1 1 .909 

-3.745 

-9.090 

0.4275 

1 . 566'9 

0.058B 

4. I 552 

23.807 

7. 559 

1 8.468 

0.2194 

1 .4602 

0.0588 

9. 1 967 

-23.807 

-7.559 

-1 8. 1 59 

0.3258 

1 . 455 I 

0.0588 

6.6254 

35.695 

1 1 .524 

27.979 

0.1642 

1-3100 

0.058.8 

1 1 . 3372 

-35.695 

-11.524 

-27.272 

0.2521 

1.3187 

0.0588 

8.8933 

47.61 1 

15.744 

37.799 

0. 1593 

1 . 1 509 

0.0508 

11.1 567 

-47.61 1 

-1 5. 744 

—3 6 . jO‘ j 

'".2101 

1 . I 540 

0.0588 

1 0. 4406 

59.650 

7. 663 

1 3. 428 


1 .0000 

0.0583 

20. 3000 

-59.650 

-20.371 

-40. 973 

0. 1906 

0.9575 

0.0538 

11.1 63/ 

72.036 

9.254 

22.200 

n 

/ m 

1 .0000 

0.0583 

20.0000 

-72.036 

-25.671 

-55.857 

0.2009 

0.7270 

0.0583 

10. 1 223 

85.300 

10.871 

26.207 

0. 

1 .'"^000 

0.0508 

20.000f) 

-85.300 

-32.206 

-24.583 

0. 

1 . 0000 

0.0508 

20. 0000 

101 .257 

12.524 

30. 220 

0. 

1 .0000 

0.05 38 

20. 0000 

-lot .257 

-4 1 . 807 

-20.073 

0. 

1 . 0000 

0.0588 

20. 0003 


FREQU1ZNCY=1600 WZ 

TRANS.MISSIf;N 

L(V 3= '-^.01 



0. 

0. 

r> , 

0.4803 

1 .6387 

0.0476 

2. 

9.52 8 

2.992 

l'.324 

0. 3972 

1.6184 

0.0476 

4 . 4 1 * 

-9.528 

-2.992 

-7.275 

0.4465 

I .5853 

0.0476 

3 . lAbo 

19.050 

6.020 

14.727 

0.2656 

1 .5197 

0.0476 

7. 6710 

-1 9.050 

-6.020 

-1 4.530 

0.3649 

1.5027 

0.0476 

5. 60 3 1 

28.562 

9. 122 

22.243 

3. 1 .33 9 

1 . 4004 

0.0476 

1 0.3020 

-28.562 

-9. 122 

-21.79 

-'.2A2.n 

1 .4036 

0.0476 

7.602 t 


Figure 38. Typical Output (Continued). 


115 



38.074 

12.344 

29.914 

0. 1599 

1 .2793 

0.0476 

1 1 .4833 

-38.074 

-12.344 

-29. 106 

0.241 5 

) .2881 

0.04/6 

9.2664 

47.61 1 

1 5. 744 

37.799 

0.1593 

1 . 1 509 

0.0476 

I I . 1567 

-47.61 I 

-15.744 

-36.505 

0.2101 

1 .1540 

0." 

•0.4406 

7.224 

9.254 

22.280 

0. 

1 . 0000 

0. , 4 76 

20.0000 

-:j/.224 

-1 9.403 

-44.054 

0.1929 

0.9.995 

0.0476 

1 1 .091 9 

67.017 

10.871 

26.207 

0. 

1 .0000 

0.0476 

20. 0000 

-67.017 

-23.445 

-51.837 

0.1903 

0.8228 

0.0476 

1 1 .0489 

77. 185 

12.524 

30.220 

0. 

1 .0000 

0.0476 

20.0000 

-77.185 

-28.085 

-28.073 

0. 

1 .0000 

0.0476 

20.0000 

8 . 162 

14.223 

34.334 

0. 

1 .0000 

0.0476 

20.0000 

-8b. »62 

-33.766 

-31 .574 

0. 

1 . 0000 

0.0476 

20.0000 

f 01 .257 

15.977 

38.565 

0. 

1 . 0000 

0.0476 

20.0000 

-101 .257 

-41 . 807 

-35.095 


1 . 0000 

0.'^476 

20.0000 


hREQUENCY=2000 HZ 

TRAMS 'll SSI ON 

LOSS= 9.08 


0. 

0. 

0. 

0.4803 

1 .6387 

0.0370 

2.9605 

7.623 

2.392 

5.854 

0.4225 

1 .6300 

0.0370 

3.9217 

-7.623 

-2.392 

-5.323 

0.4598 

1 .5998 

0.0370 

3.4636 

1 5.242 

4.802 

I 1 .754 

0.3145 

1 .5647 

0.0370 

6,3041 

-15.242 

-4.802 

-.1 1.629 

0.3935 

1 .5381 

0.0370 

4.8039 

22.856 

7.250 

17.718 

0.2273 

1 .4722 

0.0370 

8.9178 

-22.856 

-7.250 

-I 7.433 

0.3333 

1 .4649 

0.0370 

6.4229 

30.464 

9.755 

23.763 

0. 1 802 

1 .3765 

0.0370 

10.6997 

-30.464 

-9.755 

-23.252 

0.2001 

1 .3819 

0.0370 

7.9703 

38.074 

12.344 

29.914 

0. 1599 

1 .2793 

0.0370 

1 1 .4 333 

-38.074 

-12.344 

-29. 106 

0.2^15 

1 .2881 

0.0370 

9.2664 

45.6-99 

15.046 

36.201 

0. 1570 

I. 1 771 

0.0370 

1 1 . 3588 

-45.699 

-1 5.046 

-35.016 

0.2151 

1.1824 

0.0370 

10.2467 

53.364 

10.871 

26.207 

0. 

i .0000 

0.0370 

20.0000 

-53 . 364 

-17.901 

-41.012 

0. 1 9«3 

1 .0639 

0.0370 

10.8972 

61.112 

12.524 

30.220 

0. 

1 .0000 

0.0370 

20.0000 

-61.112 

-20.964 

-47. 1 33 

0. I 896 

0.9316 

0.0370 

1 1 . 1865 

69.01 1 

1 4.223 

34. 334 

0. 

I . 0000 

0.0370 

20.0000 

—69.01 1 

-24.316 

-53.433 

0. 1934 

0.7847 

0.0370 

10.8684 

77.185 

1 5.977 

38.565 

0. 

1 .0000 

0.0370 

20.0000 

-77. 1 85 

-28.085 

-35.095 

0. 

1 .0000 

0.0370 

20.0000 

85. 864 

17.797 

42.936 

0. 

I .0000 

0.0370 

20.0000 

-85.864 

-32.509 

-38.644 

0. 

1 .noon 

0.0370 

20.0000 

95.577 

0. 

0. 

0. 

1 .0000 

0.0370 

20.0000 

-95.577 

-38. 1 19 

-42.231 

0. 

1 .0000 

0.0370 

20.0000 

108.308 

0. 

0. 

0. 

1 .0000 

0.0370 

20.0000 

-108.308 

-46.903 

-45.869 

0. 

1 .0000 

0.0370 

20.0000 


FREQUENCY=2500 HZ 

TRANSMISSION 

LOSS= 9.22 


n. 

0. 

0. 

0.4803 

1 .6387 

0.0303 

2.9605 

6.050 

1 .898 

4.643 

0. 4408 

1 .6368 

0.0303 

3. 5867 

-6.050 

-1.898 

-4.623 

0.4690 

1.6104 

0.0303 

3.2692 

12.098 

3.804 

9.31 ^ 

0.3601 

1 .5972 

0.0303 

5.207-j 

-12.098 

-3.804 

-9. 2*^4 

0.4259 

1 .5653 

0.0303 

4 . 1 900 

18.143 

5.729 

14.01 ; 

0.2764 

1 .5308 

0.0303 

7.3508 

-I 8. 143 

-5.729 

-13.839 

0.3728 

1.5114 

0.0303 

5.4) iP 

24.184 

7.682 

1 8.767 

0.2165 

1 .4555 

0.0303 

9. 303 1 

-24.184 

-7.682 

-1 8.448 

0.3229 

1.4511 

0.0303 

6, 7053 

30.223 

9.675 

23.570 

0. 1 31 3 

1 .3795 

0.0303 

1 0.6591 

-30.223 

-9.675 

-23.067 

0.2816 

1 .3847 

0.0303 

7.924-4 

36.262 

11.719 

28.438 

0. 1630 

1 .3027 

0.0303 

1 1 .y!n 

-36.262 

-1 1.719 

-27.708 

0.2495 

1 .31 15 

0.0303 

8.9854 

42.307 

13.829 

33.3 88 

0.1562 

1 .2232 

0.0303 

1 1.5306 

-42.307 

-13.829 

, -32.380 

0.2255 

1 .2310 

0.0303 

9.1 51 3 


Figure 38. Typical Output (Continued). 



48.370 

16.023 

38.436 

-43.370 

-1 6.023 

-37.098 

54.465 

IA.?23 

34.334 

-54.465 

-18.325 

-41 .078 

60.616 

1 5.977 

38.565 

-60.616 

-20.762 

-46.740 

66.859 

17.797 

42.936 

-66.859 

-23.377 

-51.711 

73.248 

0. 

0. 

-73.248 

-26.227 

-42.231 

79.870 

0. 

0. 

-79.870 

-29.402 

-45.869 

86.877 

0. 

0. 

-06.877 

-33.058 

-49.571 

94.578 

0. 

0. 

-94,578 

-37.504 

-53.357 

103.812 

0, 

0. 

-103.812 

-43.581 

-53. 1 20 


8REQUENCY=3I50 HZ 

0. 

0. 

0. 

4.765 

1.494 

3.654 

-4.765 

-1 .494 

-3.642 

9.528 

2.992 

7.324 

-9.528 

-2.992 

-7.275 

1 4.290 

4.500 

11.014 

-14.290 

-4.500 

-1 0.903 

19.050 

6.020 

1 4.727 

-1 9.050 

-6.020 

-14.530 

23.007 

7.559 

1 8.468 

-23.807 

-7.559 

-1 8. 1 59 

28.562 ■ 

* 9. 122 

22.243 

-28.562 

-9. 122 

-2 1 . 795 

33.31 7 

10. 71 5 

26.056 

-33.31 7 

-10. 715 

-25.442 

38.074 

)2,34A 

29.914 

-38.074 

-1 2.344 

-29. 106 

42.837 

14.017 

33.826 

-42.837 

-14.01 7 

-32.791 

47.61 I 

1 5.744 

37. 799 

-47.61 1 

-15.744 

-36.505 

52.403 

1 7.534 

41.846 

-52.403 

-f 7.534 

-40.256 

57.224 

0. 

0. 

-57.224 

-19.403 

-44.05^ 

62.089 

0. 

0. 

-62.089 

-21.366 

-47.909 

67.01 7 

0. 

0. 

-67.017 

-23.445 

-51.837 

72.036 

). 

0. 

-72.036 

-25.671 

-55.857 1 

//. ' 85 

0. 

0. 1 

- / ' . 1 u5 

-28.085 

-53. 120 

32.520 

0. 

0. ( 

-32.523 

-30.750 

0. < 

30. 1 62 

0. 

0. ( 

-38. 162 

-33 .766 

0. ( 

94.270 

0. 

n. r 

-94.270 

-37.3! 7 

0. r 

101 .257 

0. 

0. c 

101 .257 

-41 .807 

0. 

1 10.556 

0. 

0. r 

1 I0.‘ • 6 

-48.668 

0. f 


0. 1608 

1 .1406 

0.0303 

0.2083 

1.1425 

0.0303 

0. 

1 .0000 

0.0303 

0. 1965 

1 .0458 

0.0303 

0. 

1 .0000 

0.0303 

0. 1900 

0.9404 

0.0303 

0. 

1 . OOOO 

0.0303 

0. 1902 

0.8258 

0.0303 

0. 

1 . OOOO 

0.0303 

0. 

1 .0000 

0.0303 

0. 

1 . OOOO 

0.0303 

0. 

1 .0000 

0.0303 

0. 

1 . OOOO 

0.0303 

0. 

1 .0000 

0.0303 

0. 

1 . OOOO 

0.0303 

0. 

1 .0000 

0.0303 

0. 

1 . OOOO 

0. 0303 

0. 

1 . OOOO 

0.0303 

Z TRANSMISSION 

LOSS= 9. 1 

0.4803 

1.6387 

0.0233 

0.4536 

1 .6403 

0.0233 

0.4749 

1 . 6 1 82 

0.0233 

^.3912 

1.6184 

0.0233 

0.4465 

1 .5858 

0.0233 

0.3200 

1 .5752 

0.0233 

0.4069 

1 .5466 

0.0233 

0.2656 

1 .5197 

0.0233 

0.3649 

1 .5027 

0.0233 

0.2194 

1 .4602 

0.0233 

0.3253 

1 .4551 

0.0233 

0.1889 

• 1.4004 

0.0233 

0.2920 

1 .4036 

0.0233 

0. 1 702 

1 .3404 

0.0233 

^.2640 

1 .3481 

0.0233 

0. 1599 

1 .2793 

0.0233 

0.241 5 

1 .2881 

0.0233 

0. 1561 

1.2161 

0.0233 

0. 2238 

1 .2235 

0.0233 

0. 1593 

1 .1509 

0.0233 

0,2101 

1 . 1 540 

0.0233 

0, 1800 

1 .0908 

0.0233 

0 . 1 999 

1 .0794 

0.0233 

0. 

1 . OOOO 

0.0233 

1929 

0.9995 

0.0233 

0. 

1 . OOOO 

0.0233 

0. 1 892 

0.9I4I 

0.0233 

0, 

1 .0000 

0.0233 

0. 1903 

0.8223 

0.0233 

0. 

1 .0000 

0.0233 

0.2089 

0.7270 

0.0233 

0, 

1 . OOOO 

0.02 33 

• 

1 .0000 

0.0233 

0. 

I . OOOO 

0.0233 

0. 

i . OOOO 

0.0233 

0. 

1 .0000 

0.0233 

0. 

1 . OOOO 

0.0233 

0. 

1 . OOOO 

0.0233 

0. 

1 . OOOO 

0.02 33 

0. 

1 . OOOO 

0.0233 

0. 

1 .0000 

0.0233 

0. 

1 . OOOO 

0.0233 

c. 

1 .0000 

0.0233 


^Ra)UC'JCy=4000 HZ transmission L0SS= 9.0'-> 


Figure 38. Typical Output (Continued). 


I I .04dl 
10.51 1 y 

20.0000 

1 0. 9622 

20.0000 

I I . I 783 

20.0000 

I I .0552 

20.0000 

20.0000 

20.0000 

20.0000 

20.0000 

20. 0000 

20.0000 

20. 0000 

20.0000 

20.0000 

I 

2,9605 
3.3651 
3 . I 405 
4.4152 
3.7456 
5 . 96 jO 
4.6117 
7. 67 1 0 
5.6051 
9. 1967 
6.6254 
10. 3520 
7.6024 
M. (031 
3.491 I 
1 I .4833 
9.2664 
1 1.5167 
9.91 73 
1 1 . 1567 
10. 4406 
9.9079 

I 0.8344 

20.0000 
1 1.0919 
?':.0000 
I I . i 893 
20 . "^ooo 

1 1 .0489 

20.0000 

10. 1223 

20.0000 

20.0006 

20.0000 
2^.0000 

20.0000 

20.0000 

20,0000 

20.0000 
20.0006 
20. OOOO 
2''.0000 
20.0006 
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3-STG LPTs 5.2 PR 100% 'A 


3 STAGES 

* AERO-THERMO PARAMETERS * 


STAGE 

STATION 

U- EPS 

PHI- DEG 

VR- EPS 

PS- PSIA 

TS- DEG 

1 

1 

682.000 

0. 

0. 

33.310 

722.000 

1 

2 

567.000 

62.100 

409.000 

21.180 

639.000 

2 

3 

460.000 

-46.000 

0. 

1 8.440 

619.000 

2 

4 

508.000 

62.000 

428.000 

12.550 

561 .ono 

3 

5 

461 .000 

-38.800 

0. 

10.340 

535.000 

3 

6 

444.000 

55.000 

437.000 

8.840 

516.000 


7 

464.000 

-4.700 

0. 

7.750 

490.000 



***** EIRST CUT-ON 

OCCURS AT 

1 86. HZ ***** 

THETA-I 

0. 

THETA-R THETA-T 

0. 0. 

'T 

1 

0.4585 

8 E T-L0S3 

1.5652 1.0000 1.8801 


FREQUENCY^ 50 

HZ TRANSHISSIO!! LOSS= 1.38 


0. 

0. 

n. 

0.4585 1.5652 

1 . 0000 

1 . 8801 


FREQUENCY^ 

63 

HZ TRAUSMTSSIOH 

L0SS= 1 . 88 


0. 

0. 

0. 

0.4585 1.5652 

1 .0000 

1 . 880 1 


FREQUEMCY= 80 dZ TRAHSTISSION LOSS= I . 88 


O. 0. 0. 0.4585 1.5652 1.0000 1.8801 

EREQUEHCY= 100 HZ TRANS .IISSION LOSS= 1.88 

O. n. 0. 0.45-Qr, 1.5652 1.0000 1.8801 

rREQUENCY= 125 rZ i hANSMISSION L0S5= 1.88 

0. 0. 0.4585 1.5652 1.0000 1.8801 

r ■•.'.JUENCY= '60 HZ TRANSMISSION LOSS= 1.88 


O. 


0. 0.4o85 1.565' O.J 1.8801 

Figure 38. Typical Output (Continued). 
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ORIGINAl. PAGii.; 
OF P()OR QUALITx’’ 


iOl .232 

1.494 3.654 0. 

1 .0000 

0. 3333 

— - 

20. OOOO 

-101 .232 

-41.708 -3.642 0. 

1 . 0000 

0.3333 

20.0000 


FRCOUENCY= 200 HZ 

TRANSMISSION 

L0SS= 6.52 



0. 

0. 

0. 

0.4585 

1.5652 

0.3333 

1 . 8801 

77.198 

1 .494 

3.654 

0. 

1 . 0000 

0. 3333 

20.0000 

-77.198 

“28.034 

-3.642 

0. 

1 .0000 

0.3333 

20. OOOO 


FREQUHNCY= 250 HZ 

rRANSMISSIOM 

LOSS= 6.52 


0. 

0. 

0. 

0.4585 

1 .5652 

0. 3333 

1 . 880 1 

60.635 

1 .494 

3.654 

0. 

1 . 0000 

0. 3333 

20.0000 

-60.635 

-20.727 

-49.338 

0. 1522 

0.9530 

0.3333 

1 1.6097 


FREQUENCY^ 315 HZ 

TRANSMISSION 

LOSS= 6.15 


0. 

0. 

0. 

0.4585 

1 .5652 

0.2000 

1 . 880 1 

47.629 

15.717 

37. 148 

0.1192 

1.1619 

0 . 20^0 

12.2963 

-47.629 

-15.717 

-38.382 ■ 

0.1674 

1 .1630 

0.2000 

10. 880 7 

101 .232 

2.992 

7.324 

0 . 

1 .0000 

0.2000 

20. 0000 

-101 .232 

-41 .703 

-7.275 

0. 

I . 0000 

0.2000 

20.0000 


FREQUENCY= 400 HZ 

TRANSMISSION 

LOSS= 7.91 


0. 

0. 

0. 

0. 4885 

1 .5652 

0.2000 

1 . 880 1 

38.091 

12.324 

29.710 

0. 1 1 88 

1.2882 

0.2000 

12.6840 

-38.091 

-12.324 

-30.490 

0. 1947 

1.2910 

0. 2000 

9. 594 -3 

77.198 

2.992 

7.324 

0. 

1 .0000 

O.2OO0 

20. ooOO 

-77.198 

-28.034 

-7.275 

0. 

1 .0000 

0.2000 

20.000) 


FREQUENCY= 500 HZ 

TRANSMISSION 

LOSS= 7.79 



0. 

0. 

0. 

0.4585 

1 .5652 

0.1429 

1.8801 

30.237 

9.659 

23.592 

0. 1351 

1.3775 

0.1429 

1 1 . 81 56 

-30.237 

-9.659 

-24.080 

0.2324 

1 .3791 

0. 1 429 

8.04y > 

60.635 

2.992 

7.324 

0. 

1 .0000 

0. 1429 

ZO.noO' 

-60.635 

-20.727 

-49.338 

0.1522 

0.9530 

0. 1 429 

1 1 .6097 

94.570 

4.500 

11.014 

0. 

1 .0000 

0. 1429 

? 0.0000 

—94 . 570 

-37.426 

-1 0.903 

0. 

1 .0000 

0.1429 

20. 0000 


FREQUENCY= 630 HZ 

TRANSMISSION 

CO 

11 

00 

• 

CD 

CO 



0. 

0. 

0. 

0.4585 

1 .5652 

0. 1 1 II 

1 . 880 1 

23.81 8 

7.547 

1 8.595 

O. 1686 

1.4415 

0.1111 

10.0769 

-23.818 

-/.547 

1 8.896 

0.2769 

1 .4392 

0. 1 ' M 

6.5020 

47.629 

15.71 7 

37.148 

0.1192 

1.1619 

0. 1 1 ■ 1 

12.2968 

-47.629 

-15.717 

38.382 

0.1674 

1 . 1 630 

0. 1 Ml 

10.8807 

72.051 

4.500 

11.014 

0. 

1 . OOOO 

0. 1 1 1 1 

20.0000 

-72.051 

-25.626 

1 0.903 

0. 

1 . OOOO 

0. 1 1 I 1 

20.0000 

101 .232 

6.020 

1 4.727 

0, 

1 .0000 

0. I 1 1 1 

20.0000 

101 .232 

-41.708 

1 4.530 

0 . 

1 .0000 

0. 11 11 

20.0000 


FREQUENCY 

= 800 HZ 

TRANSMISSION 

LOSS= 8.93 



Figure 38. Typical Output (Continued) . 
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^ A*" 


0. 

0. 

0. 

0.4585 

1 .5652 

0.0909 

1 . 880 1 

1 9.059 

6.010 

14.889 

0.2137 

1 .484''' 

0.0909 

8. 1364 

-19.059 

-6.010 

-15.081 

0.3183 

1 .4770 

0.0909 

5.2606 

38.091 

12.324 

29.710 

0.1188 

I .2882 

0,0909 

12.6840 

-38.091 

-12.324 

-30.490 

0.1947 

1 .2910 

0.0909 

9.5943 

57.243 

4.500 

11.014 

0. 

1 .0000 

0.0909 

20.0000 

-57.243 

-19.370 

-46.460 

0.1538 

1 .OJ 15 

0.0909 

1 1.5525 

77.198 

6.020 

14.727 

0. 

1 .0000 

0.0909 

20. 0000 

-77.198 

-28.034 

-14.530 

0. 

1 .0000 

0.0909 

20.0000 

101 .232 

7.559 

18.468 

0. 

1 .0000 

0.0909 

20,0000 

-101 .232 

-41.708 

-1 8.159 

0. 

1 .0000 

0.0909 

20.0000 


FREQUENCY=I000 HZ 

TRANSMISSION 

L()SS= 9.00 


0. 

0. 

0. 

0.4585 

1 .5652 

0.0769 

1.8801 

1 5.250 

4.795 

1 1 .921 

0.2658 

1.5143 

0.0769 

6.3283 

-15.250 

-4.795 

-12.044 

0.3555 

I .5031 

0.0769 

4.2709 

30.478 

9.739 

23.780 

0. 1343 

1 .3749 

0.0769 

11.8613 

-30.478 

-9.739 

-24.276 

0.2310 

1 .3766 

0.0769 

8. 1028 

45.717 

15.021 

35.657 

0.1176 

1 . 1 387 

0.0769 

12.4894 

-45.717 

-15.021 

-36.791 

O. 1716 

1 . 1 905 

0.0769 

10.6750 

61.130 

6.020 

14.727 

0. 

1 .0000 

0.0769 

20.0000 

-61 . 130 

-20.929 

-49.759 

0. 1514 

0.9442 

0.0769 

1 1 .6474 

77. 198 

7.559 

1 8.468 

0. 

1 .0000 

0.0769 

20.0000 

-77. 198 

-28.034 

-18.159 

0. 

1 ,0000 

0.0769 

20.0000 

95.567 

9. 122 

22.243 

0. 

1 .0000 

0.0769 

20.0000 

-95.567 

-38.038 

-2 1 . 795 

0. 

1 . 0000 

0.0769 

20.0000 


• FREQUENCY=1250 HZ 

TRANSMISSKH'I 

L()SS= 8.74 



0. 

0. 

0. 

0.4585 

1 .5652 

0.0588 

1 . 880 1 

11.915 

3.739 

9.321 

0.3209 

1.5367 

0.0588 

4.766'^ 

-11 .915 

-3.739 

-9.396 

0.3890 

1 .5228 

0.0583 

3.4586 

23.81 8 

7.547 

18.595 

0. 168o 

1.4415 

0.0588 

10.0769 

-23.81 8 

-7.547 

-1 8.896 

0.2769 

1 .4392 

0.05B8 

6.5020 

35.71 1 

1 1.505 

27.856 

0. 1218 

1 .3166 

0.0588 

12.5478 

-35.711 

-1 1.505 

-28.539 

0.2044 

1 .3194 

0.0588 

9. 1 726 

47.629 

15.717 

37. 1 48 

0.1192 

1.1619 

0.0583 

12.2968 

-47.629 

-15.717 

-38.382 

0. 1674 

1 . 1 630 

0.0588 

10. 8807 

59.669 

7.559 

I 8.468 

0. 

1 . 0000 

0.0588 

20.0000 

-59.669 

-20.336 

-48.517 

0.1526 

0.9699 

0.05 88 

( 1 .5977 

72.051 

9. 122 

22.243 

0. 

I .0000 

0.0588 

20.0000 

-72.051 

-25.626 

-21.795 

0. 

I .0000 

0.0538 

20.0000 

85.306 

10. 715 

26.056 

0. 

1 .0000 

0.0583 

20.0000 

-85.306 

-32. 1 44 

-25.442 

0. 

1 .0000 

0.0588 

20.0000 

101 .232 

12.344 

29.914 

0. 

1 . 0000 

0.0583 

20.0000 

101 .232 

-41.708 

-29. 106 

0. 

1 .0000 

0.0588 

20.0000 


FREQUENCY=1600 HZ 

TRANSMISSION 

L05S= 8.36 



0. 

0. 

0. 

0.4585 

1 .5652 

0.0476 

1 . 880 1 

9,533 

2.988 

7.461 

0.3615 

1 .5495 

0.0476 

3.7791 

-9.533 

-2.988 

-7.509 

0.411 8 

1 .5350 

0.0476 

2.9379 

1 9.059 

6.010 

1 4.889 

0.2137 

1 .4840 

0.0476 

8. 1364 

-19.059 

-6.010 

-15.081 

0.3183 

1 , 4770 

0.0476 

5.2606 

28.576 

9. 107 

22.299 

0,1415 

1 .3948 

0.0475 

I 1 .4644 

-23.576 

-9. 107 

-22.734 

0.2426 

1 .3957 

0.0476 

7. 6696 


Figure 38, Typical Output (Continued) . 
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38.091 

12.324 

29.710 

0.1188 

1 .2882 

0.0476 

1 2 . 6840 

-38.091 

-12.324 

-30.490 

0. 1 947 

1 .2910 

0.0476 

9.5943 

47.629 

15.717 

37. 148 

0.1195 

1.1619 

0.0476 

12.2724 

-47.629 

-15.717 

-38.382 

0. 1674 

1 . 1 630 

0.0476 

10.8822 

57.243 

9.122 

22.243 

0. 

1 . 0000 

0.0476 

20.0000 

-57.243 

-19.370 

-46.460 

0. 1538 

1.0115 

0.0476 

1 1 . 5526 

67.034 

10.71 5 

26.056 

0. 

1 .nooo 

0.0476 

20.0000 

-67.034 

-23.404 

-54.801 

0.1547 

0.8355 

0.0476 

1 1 .3803 

77. 198 

12.344 

29.914 

0. 

1 . 0000 

0.047d 

20.0000 

-77.198 

-28.034 

-29.106 

0. 

1 .0000 

0.0476 

20.0000 

88. 164 

14.01 7 

33.825 

0. 

1 . 0000 

0.0476 

20.0000 

-88.164 

-33.701 

-32.791 

0. 

1 .0000 

0.0476 

20.0000 

101.232 

15.744 

37.799 

0. 

1 .0000 

0.0476 

20.0000 

101 .232 

-41 .708 

-36.505 

0. 

1 .0000 

0.0476 

20.0000 


0 . 

7.627 
-7.627 
15.250 
-1 5.250 
22.867 
- 22.861 
30.478 
-30.478 
38.091 
-38.091 
45.717 
-45.717 
53.383 
-53.383 
61.130 
-61 .130 
69.028 
-69.023 
77. 198 
-77. IvB 
85.870 
-85.870 
95 . 567 
-95.567 
108.245 
108.245 


8REQUENCY=2000 HZ TRANSMISSION LOSS= 3.76 


0. 

0. 

0.4585 

2.388 

5.972 

0.391 8 

-2.388 

-6.003 

0.4283 

4.795 

11.921 

0.2658 

-4.795 

-1 2.044 

0.3555 

7.238 

1 7.854 

0. 1 759 

-7.238 

-18.131 

0. 2846 

9.739 

23.780 

0. 1343 

-9.739 

-24.276 

0.2310 

12.324 

29.710 

0.1187 

-12.324 

-30.490 

0. 1947 

1 5.021 

35.657 

0. I 176 

-15.021 

-36.791 

0.1716 

10.715 

26.056 

0. 

-1 7.871 

-43.203 

0. 1578 

12.344 

29.914 

0. 

-20.929 

-49. 759 

0.1514 

14.01 7 

33.825 

0. 

-24.273 

-56.51 1 

0. 1588 

15.744 

37.799 

0. 

-28.034 

-36.505 

0. 

1 7.534 

41 .846 

0. 

-32.446 

-40.256 

0. 

0. 

0. 

0. 

-38.038 

-44.054 

0. 

n. 

0. 

0. 

-46.766 

-47.909 

0. 


1 .5652 

0.0370 

1 . 880 1 

1 .5572 

0.0370 

3.1164 

1 .5436 

0.0370 

2.574! 

1 .5143 

0.0370 

6.3283 

1 .5031 

0.0370 

4.2709 

1 .4503 

0.0370 

9.7300 

1 .4472 

0.0370 

6.2573 

1 .3749 

0.0370 

11.8613 

1 .3766 

0.0370 

8. 102 8 

1 .2882 

0.0370 

12.6880 

1 .2910 

0.0370 

9.5949 

1 . 1 887 

0.0370 

12.4894 

1 . 1 905 

0.0370 

10.6750 

1 .0000 

0.0370 

20.0000 

1 .0750 

O'. 05 70 

1 I .3592 

1 .0000 

0.0370 

20. 0000 

0.9442 

0.0370 

1 1.6474 

1 . 0000 

0.0370 

20.0000 

0.7971 

0.0370 

1 1 . 1200 

1 . 0000 

0.0370 

20.0000 

1 .0000 

. 0.0370 

20.0000 

1 .0000 

0.0370 

20.0000 

1 .0000 

0.0370 

20.0000 

1 .0000 

0.0370 

20.0000 

1 . 0000 

0.0370 

20.0000 

1 .0000 

0.0370 

20.0000 

1 .0000 

0.0370 

20. 0000 


TRE0UENCY=2500 HZ TRANSMISSION L()SS= 8.88 


0. 

• 

0. 

0. 4585 

1 .5652 

0.0303 

1 . 8801 

6.053 

1.895 

4.742 

0.41 37 

1 .5619 

0.0303 

2.6714 

-6.053 

-! .895 

-4.761 

0.4401 

1 .5499 

0.0303 

2.3192 

12.104 

3.798 

9.468 

0.3176 

1 .5356 

0.0303 

4.8503 

-12. 104 

-3.798 

-9.546 

0.3871 

1.5217 

0.0303 

3.5023 

18.152 

J.720 

14. 1 83 

0.2248 

1.4916 

0.0303 

7.7179 

-18.152 

-5.720 

-14.357 

0.3269 

1 .4835 

0.0303 

5.0223 

24.196 

7.670 

1 0. 889 

0. 1658 

1 .4379 

0.0303 

10.2081 

-24.196 

-7.670 

-1 9.200 

0.2739 

1 .4360 

0.0303 

6.5984 

30.237 

9.659 

23.592 

0.1351 

1 .3775 

0.0303 

11.8156 

-30.237 

-9.659 

-24.080 

0. 2324 

1.3791 

0.0303 

8.0492 

36.277 

1 1 . 699 

28.297 

0. 1209 

1.3100 

0.0303 

12.5894 

-36.277 

-1 1 .699 

-29.003 

0.201 9 

1 .3128 

0.0303 

9.2768 

42.325 

1 3.806 

33.01 1 

0. 1 165 

1 .2345 

0.0303 

12.6987 

-42.325 

-13.806 

-33.978 

0. 1 805 

1 .2370 

0.0303 

10.2453 


Figure 38. Typical Output (Continued). 
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48.3 88 

15.997 

37.741 

0.1223 

I . 151 I 

0.03C3 

12.0433 

-48.388 

-15.997 

-39.016 

0.1659 

1 . 1519 

0.0303 

10.9557 

54.484 

14.01 7 

33.825 

0. 

I . 0000 

0.0303 

20.0000 

-54.484 

-18.294 

-44. 129 

0. 1566 

1 .0572 

o.o:’o3 

1 1 .4078 

60.635 

15.744 

37.799 

0. 

1 . 0000 

0.0303 

20.0000 

-60.635 

-20.727 

-49.338 

0.1522 

0.9530 

0.0303 

1 1 .6083 

66.877 

1 7.534 

41 .846 

0. 

1 .0000 

0.0303 

20.0000 

-66.877 

-23.336 

-54.666 

0.1542 

0.8385 

0.0303 

11.4140 

73.263 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

-73.263 

-26. 180 

-44.':)54 

0. 

1 .0000 

0.0303 

20.0000 

79.S81 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

-79.081 

-2 9.34 8 

-47.909 

0. 

1 .0000 

0.0303 

20.0000 

86.881 

0. 

0. 

0. 

I .0000 

0.0303 

20.0000 

-86.881 

-32.995 

-51.837 

0. 

1 .0000 

0.0303 

20.0000 

94.570 

0. 

0. 

0. 

1 . 0000 

0.0303 

20.0000 

-94.570 

-37.426 

-55.857 

0. 

1 .0000 

0.0303 

20.0000 

J03.776 

0. 

0. 

0. 

1 .0000 

0.0303 

20.0000 

103.776 

-43.472 

-53. 120 

0. 

1 .0000 

0.0303 

20.0000 


h'REQUENCY=31 50 HZ 

fR A'VSVlISSlnN 

L()SS= 8.77 



0. 

0. 

0. 

'^.4585 

1 .5652 

0.0233 

1 . 880 1 

4.767 

1 .492 

3.736 

0.4289 

1 .5645 

0.0233 

2.381 5 

-4.767 

-1 . 492 

-3.748 

0. 4481 

1 . 5544 

0.0233 

2. 1 472 

9.533 

2.983 

7.461 

0.3615 

1 .5495 

0.0233 

3.7791 

-9.533 

-2.988 

-7.509 

0.4118 

1 .5350 

0.0233 

2.9379 

I -.297 

4.492 

11.178 

0.2809 

1.5212 

0.0233 

5.8703 

-1-.297 

-4.492 

-1 1 .286 

0.3651 

1 .5090 

0.0233 

4.031 6 

1 . -5y 

6.010 

1 4. 889 

0.21 37 

1 .4840 

0.0233 

8.1364 

— .'59 

-6.010 

-1 5.081 

0.3183 

1 .4770 

0.0233 

5.2606 

0.318 

7.547 

1 8.595 

0.1685 

1.4415 

0.0233 

10.0769 

-23.81 8 

-7.547 

-1 0.896 

0.2769 

1 .4392 

0.0233 

6.5020 

28.576 

9. 107 

22.299 

0. 141 o 

1 .3948 

0.0233 

1 1 .4644 

-28.576 

-9. 107 

-22.734 

0.2426 

1 .3957 

0.0233 

7.6696 

33.332 

10.697 

26.003 

0.1263 

1 .3438 

0.0233 

12.3053 

-33. ^32 

-10.697 

-26.597 

0.2! 

1 .3463 

0.0233 

8.7093 


Figure 38. Typical Output (Continued) . 
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J 8 . 09 I 

12.324 

29.710 

0 . 

- 38.091 

- 12.324 

- 30.490 

0 . 

42.854 

1 3. 994 

33.424 

0 . 

- 42.854 

-1 3.994 

- 34.416 

0 . 

47.629 

1 -J .71 7 

37.140 

0 . 

- 47.629 

- 15.717 

- 38.382 

0 . 

52.422 

1 7.505 

40. 888 

0 . 

- 52 . 42 ? 

- 17.505 

- 42.394 

0 . 

57.243 

0 . 

0 . 

0 . 

- 57.243 

- 1 9. 370 

- 46.460 

0 . 

62.107 

0 . 

n. 

O. 

- 62 . 1 07 

- 21 .329 

- 50.592 

0 . 

67.034 

0 . 

0 . 

0 . 

- 67.034 

- 23.404 

- 54.801 

0 . 

72.051 

0 . 

0 . 

0 . 

- 7. '.051 

- 25.626 

- 55.857 

0 . 

77. 1 98 

0 . 

0 . 

0 . 

- 77. 1 98 

- 23.034 

- 53.120 

0 . 

82.536 

0 . 

0 . 

0 . 

- 32.536 

- 30.693 

0 . 

0 . 

38. 1 64 

0 . 

0 . 

0 . 

- 88.164 

- 33.701 

0 . 

0 . 

94.263 

0 . 

0 . 

0 . 

- 94.263 

- 37.239 

0 . 

0 . 

m .232 

0 . 

0 . 

0 . 

- 101 .232 

- 41 .703 

0 . 

0 . 

1 10.472 

0 . 

0 . 

0 . 

- 1 1 0.472 

- 48.510 

0 . 

0 . 


FREQUENCY =4000 HZ TfJAH^ 


1 88 

1 .2802 

0.0233 

12.6840 

947 

1 .2910 

0.02 '3 

9.5948 

168 

1 .2275 

0.0233 

12.6581 

790 

1 .2299 

0.0233 

10.3182 

1 95 

1.1619 

0.0233 

12.9724 

674 

1 . 1 630 

0.0233 

10. 8822 

380 

1 .0939 

0.0233 

10.8213 

592 

1 .0902 

0.0233 

1 1 .2929 


1 .0000 

0.0233 

20.0000 

538 

I . 01 1 5 

0.0233 

I 1 .5526 


1 .0000 

0.0233 

20.0000 

521 

0.9267 

0.0233 

1 1 .5943 


1 . 0000 

0.0233 

20.0000 

547 

0.8355 

0.0233 

1 1 .3808 


1 .0000 

0.0233 

20. 0000 


1 . 1 ’ 

0.0233 

20.0000 


1 . 0000 

0.0233 

20.0000 


1 .0000 

0.0233 

20.0000 


1 .0000 

0.0233 

20.0000 


1 . '^000 

0.0233 

20.0000 


1 . 0000 

r ^.0233 

20. 0000 


1 .0000 

0.0233 

20.0000 


1 .0000 

0.0233 

20. 0000 


1 .0000 

0.0233 

20.0000 


1 .0000 

0.0233 

20.0000 


1 .0000 

0.0233 

20. 0000 


1 .0000 

0.0233 

20. 0000 


1 .0000 

0.02 33 

20.0000 

MISSION L 0 SS= 

8.84 



Figure 38, Typical Output (Concluded). 
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APPENDIX B 


GENERALIZED ITERATIVE PROCEDURE COMPUTER PROGRAM 


The matrix inversion procedure cannot be used for turbines containing 
choked blade rows because the matrices decouple at these rows. The problem 
is really that it becomes impossible to determine the three downstream 
(Pm> 2jjj, Qjn) because one of the three available equations is 
independent of the downstream amplitudes for a choked row. However, it 
remains possible to calculate F^, Bn, and Qm in terms of Fn, Bm, and’p . 

The solution procedure utilized here exploits this fact. The value of’^Bm is 
guessed at each choked blade row. The amplitudes Fn and Qn from the computa- 
tions for the preceding row are used in the choked flow equation to calculate 
^ ®n** If Bn* is not the same as the Bn from the immediately preceding 
computations, the program returns to the last blade row where B was guessed, 
and a new guess is made. If Bn* agrees with Bn, then F^ and 0^ are calculated 
using the other two equations available at that blade row, and the program 
continues to the next row. The final verification of a correct guess is made 
at the last blade row where the assumption of an anechoic termination requires 
Bn, - 0. 

This procedure is summarized in the flow chart of Figure 39. Each 
choked blade row entails verification of the last guess in the form of 
®n 7 Bn*? Once this is achieved, a new guess (Bn,) follows. The final 
verification is (Bm)last stage = 0* The one deviation is at the first blade 
row; if this is unchoked, the program guesses at Bj. Since Fi = 1 Qi = 0 
P2> ^ 2 , and Q 2 can be calculated. ’ ’ 

A listing of the program can be found in Figure 40, and a sample output 
is provided in Figure 41. Both input and output frequency distribution and 
energy assignment are the same as in the matrix inversion program. 
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P^GBJS 


Figure 39. Flow Cha^t - Mutistage, Multimode Computer Program 
Using Iterative Solution. 
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********** MLTGS3 ********** 


0010*#»<UNH *; 


OO 2 0*# L I b A R Y /MT I N V , H= { UL I B ) U S E RL I B/TDS , R 
0030C ******* FILENAME MLIGS3 ******* 

0040 CUMMUN /CANGP / PI , TGDEG,Ti)RAU 

OObO COMMON /CINPUT/ NSTAGE, lOPT, IPRINT,PT0,TTU,STAGEX(5 , 15) . 

OU60^ NSTA, TITLE 

0070 CHARACTER TITLE*40 


0080 COMMON /CIJTOFA/ THCL,THCU 

009l common /CAERO / V(lb),MX(lb) ,MY(15),AMQAN(lb),RHDRAT(15) , 

OlOU^i MACH(lb),AS(15) 

110 WEAL KiIQKM(lb) ,MX,MY,MACH,MABS 
0120 INTEGEh FREU(20) 

01 3u EQUIVALENCE (KNQKM, A.MQAN ) 


0140 
Ol'jO 
0100& 
01 /O 
OlbO 
01 9U 


COMMON 

COMMON 

COMMON 

COMMON 


/CAEROI/ 
/COUT / 

/CMATRX/ 

/CATTCH/ 


U(lb) ,PHKlb),VR(lb) ,P( lb),TClb),GAM( lb) 
TLOSS(IOO), THI(100),THR(100) , THT(IOO), 
Q(100),B(100),B1(100) ,TW(100),F(100) 
Ll(3,3,15),DI(J,3,15),A(3,3»lb),PR0D(3,3,lb) 
CF1,CF2,CF3,BUF(380) 


02 UO 

02ub 

0210 

0220 

0230 

0240 

02b0 

02O0& 

0€70 

0280 

0290 

0300 

0310 

031bd, 

0320 

0330 

0340 

03b0 

0360 

03/0 

0380C 

0390 

0400 

0410 

0420 

0430 

0440 

04:d0 


DIMENSION STA(iEP(/b) 

DIMENSION SMACH(Id) 

DIMENSION E(IOO) 

EQU I VALENCE ( STAGER ( 1 ) , S I AGEX (1,1)) 

CHARACTER CF1*1 CF2*8, CF3*l/“ ; •'/ 

CHARACTER f ITLE*40, bLANK*40 

REAL MACHN,MACHM,KMYKM,KNV-:n,KMXKM,KNXKiJ,KMYSAV,KMXSAV, 
KNN.KMM 

EQU I VALEN CE ( I B I IS , b I TS ) 

DATA BITS/D3777 / /7 r/777/, 0KA/04 (j 3 /OOOOOOOO/, Jp0/004007b040007/ 
DATA R0,SIGMA/1 I .6,0.882/ 

DAPA FREQ/b0,63,Bu,l{)0, 12b, 160,200,250, 31b, 400, bOO, 

630,800,1000, 1250, 1600,2000,2500,3150,4000/ 

DATA PI, T0DEG,T0RAU/3. 1415927, 57. 29578, .0174532925/ 

DATA BLANK/" i. / 

NAMELIST /TNOISE/ I0PT,PT0, T'l J, STAGER, TI TLE, GAM, lAERO 
TAN(X)== SIN(X)/COS(X) 

SET UP NAMELIST INHUT FILE 

lAERO = 0 

CALL FPARAM(3,JP0) 

PRINT," INPUT FILE NAME " 

READ,CF2 

CALL ATTACH ( 1 , CFl , 1 , 0, 5TAT, blJF) 

IF( STAT.EQ.O. .OR. STAT.EQ.OKA ) GO TO d 


Figure 40. Program Listing - Generalized Iterative Procedure. 
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I 



0460 

04 to 
0480 
0490 
OdOOC 
0510 
0520 
0530 
0540 
0550 
0560 

05 /O 
0580 
05 VO 
0600 
061 OC 
0620 
0630 

0650 
0060 
0670 
OotsO 
06V0 
0/00 
0710 
0/20 
0/30 
') /40 
0750 
0760 
0 / / O 
0/80 
0 7 VO 
0800 
0810 
0820 & 
0«30ei 
0o40 
0850 
0860 
08/0 
088o 
08V0 
OVOO 
( Viu 
OV20 
OVJO 
0V40& 
0V50 
0V60 


Figure 


PRINT 1,STAT 

1 RORMATC INPUT FILt STATUS=*' , 012 ) 

STOP 

INI ri ALI zati qn ********* ********** 

5 IGPl = 1 
DO 10 1=1, l3 
STAGEX(1,I)= BITS 
uam(I)= bits 
10 CONTINUE 

PTO = 14.6V6 
TTJ == 518.7 
niLE = BLANK 

i/EAU INPUT FILE ****■ COUNT NH. OF STATIONS 
15 ( I , TNOISE, END = 400) 

DO 1/ 1=1,15 

IFCSTAGEXd, I) .EQ.bl T5 ) GO TO 18 
U(I) = STAGEXd,!) * 3.048 

PHI (I) = STAGEX(2,I) 

VRd) = STAGEX(3,I) * 3.048 

PCI) = STAGEX(4,n / 6.8V5 

I'd) = STAGEX(5.I) * 1.8 

I F ( GAM d ) . N E . B I rs ) GO 1 0 17 
GAM(1)= GAMXd'd)) 

17 CONTINUE 

18 NS I A = I-l 

IF( NSTA.EG.14 ) NSiA=15 
NSTAGE= (NSTA-D/2 
PRINT 21, TITLE,,/ STAGE 

21 F0nMAT(//16a,A40//32X. 12," STAGES"//) 

PRINT 22 

22 F0RMA1'(28X, "* AERO-iHEkMO PARAMETERS *"// 

2X, "STAGE",3X, "SiAlT0N",3X, "U- FPS",3X, "PHI- DEG", 
3X,"VR- FR:v',:^X,“i-S- PS I A", 2X, "TS- DEG ri"/) 

NSTG = 0 
00 24 l=l,NSlA 

1F(( (I/2)*2)/I.EQ.0 ) NSTG=NSTG+1 
IF( I.EQ.N iA ) NSi..= IBITS 

PRINT 23,NSxu,l,'J' ') ,PHId), VR(I),Pd),T(I) 

23 F0RmAT(4X, 1 1 , 7a, d.rl2. 3,4F10.3) 

24 CONTINUE 

IF( IAEhO.EO.O. ) vjU to 27 
PRINT 25 

25 FQRMAT(/2X, "STAGE", 3X, "STATION", 5X, "MX", 8X, "MY", 

/X, "MACH",SX, "KNOKM", 7X, "V"/) 

NSTG = 0 
00 26 I=1,NSTA 


40. Program Listing - Generalized Iterative Procedure 
(Continued) . 


Oy/u 

OVbU 

oyyo 
1000 
1010 
1020 
1030& 
1040 
lObOC 
1000 . 
10 70 ' 

lut . ■ 

1 Oyy 
1100 
1110 
1120 
1130 
1140 
1130 
Iloo 
Ilob 
11/0 
llbO 
1190 
1200 
1210 
1220 
12302. 
1240 
1200 
1200 
12 10 
12bO 
1290 
1300 
1310 
1320 
13 JO 
1 J40 
13bO 
13o0 
1370 
IJoO 
1390C 
1400 
1410 
1420 
1430 
1440 
14bO 
1400 


Figure 


IF( ( (I/2)*2)/I.t'Q.O ) NSTG=ilSlJ+l 
IF( I.EQ.l ) GO TO 26 

FRINf 23, N3TG, I ,MX( I ) ,MY( I ) ,MACH( I ) ,K.JOKM( I ) , V( I ) 

2o CONTINUE 
2/ CONTINUE 

28 FORMAT ( //2X , ^THETA- 1 3X , 'THhTA-R-', 3X, ^THETA-T' , 6X, 

, 'T^,9X,'b', .a,'E', 6X, 'T-LUSS^/l 

CALC^tE AEwO-THERMO PARAMETERS 

‘ 'C '-' v .- - 

, 'gfcf^'-2 9’ I=1,NSTA 

*‘AS(I) = 41 .42*SQrT(GAM( I )*T( I) ) 

MX ( I ) = U ( I ) /AS ( I ) 

IF( I.EQ.l ) GO 0 29 ' 

AMQAN ( I )= SQRT(oAM( I )*T( I )/(GAM( I-l )*T( I-l ) ) ) 

RHOH ATI I )= 1 ( I) »F( I-l )/ri ( I-l )*P( I ) ) 

29 CONlIHUE 

ASl = AS(i) 

fMEAN= R0*SQRT( ( 1 .+SIGMA**2 )/2. ) 

RMEAN = RMEAN/2.b4 
fFI = 1. 

FREOCO^ ( (FF1*(AS1*12. ) )/(2.*Pl*RMEAN 1 )*SQRT( 1 .-XM1**2) 
FREQC1= AINT(FREQCO) 

PRINT 30,FREQC1 
PRINT 28 

30 FORMATI/ZIOX,'**** FIRST C'Ji'-UN OCCURS Af ' , 1 X , F4 .'), IX , 

'HZ »***'///) 

JO 300 L=l,20 

FP = (2.*PI>-pi^EQ(L)*RMEAiJ)/( AS1*12. ) 

XMl = MX(l) 

NTH=FP/SORT( 1 . -XMl ** 2 ) 

IF( NTH.GT.50 ) N rn=bO 
NTT = 2*NTH+1 
FRSQ= 0. 
rHKl)= 0. 

JO 32 J=1,NTT 
E(J)= 0, 

32 CONTINUE 
E(l)= 1. 

ESlGMA= 1. 

IHI (1)= 0. 

IF( NTH.LT.l ) GO 10 bO 
**** COMPUTE CU'I-ON MODES, ANGLES, AND ENERGY 
DO 40 J=l,Nri-i 
FJ = J 
F1=FP/FJ 

F2 = SORT! F1**4-F1**2*( 1 .-XM1**2 ) ) 

F3 = XMl **2+F 1**2 
F4 ^ (F2-XMD/F3 
J1 = 2*J 


40. Program Listing - Generalized Iterative Procedure 
(Continued) . 
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I 



14/0 
I4b0 
1490 
1 -jOO 
IblO 
lb20 
l!330 
1040 
1640 
lObO 
1 660 
16 70 
16bO 

lovo 
1 /OOC 

1 no 

1720 
1730 
1740 
1 /bO 
1 /60C 
1 no 
1 /80 
1820 
1830 
1840 
18b0 

lovo 

189b 
1 8v6 
189/ 
1898 
1900 
IvlO 
1920 
1930 
1940C 
I9b0 
19o0 
19/0 
1980 
1 990 
2000 
2010 
2020 
2030 
2040 
20b0 
2060 
2070 
2080 


J2 = 2*J+1 
1HI(1)= 0. 

rHI(Jl)= TGUE0*At<CUS(F4) 
iHl(J2)= -THKJl) 

FC = (FJ*AS1*12. )/(2.*PI*HMEAN) 
FC = FC*5Q8T( l.-XMl**2) 

FkSO = FC/FHtU(L) 

1F( Fl/SQ.GT. 1.02b ) GG lu 40 

t(Jl) = 1. 

t(J2)= E(J1) 

h5l0MA= bS I GMA+2 . *b l J 1 ) 

4u GGlFl INUb 
bO CGNlINUb 



**** COMPUTE b'Nb'KGY UISTHIBUlTUN 
00 60 K=1,J2 
b'(K>= b(K)/ESlGMA 
60 COiniNJb 


****** INNbR LOOP TO BUILD MATRICES ***** 

SUmT= 0. 

00 240 K=1,NTT 
62 I’HFN = THI(K)*TOr/AO 
KMYSAV= 0. 
iCMXSAV= 0. 

OGS = 1 
iMISIAK = 1 
F'(l) = 1.0 
8(1) = -30. 

0 ( 1 ) = 0 . 

M = 1 

6b M = M+1 

1F( M.GT.NSfA ) GO TO 9b 

CALC' 'LATE ANGLES AND RATIOS 

6o N = M — 1 

V(N) = U(N)*TAN(TORAO*PHI(N))-VR(N) 

V(M) = U(M)*TAN(T0HAD*l-:iI(.4) )-VR(N) 

MY(N) = V(n)/AS(N) 

MY(M) = V(M)/AS(M) 

MACH ( N ) = SO RT ( M X ( N ) **2+M Y ( N ) **2 ) 

MACH(M)= SQRT(MX(M) **2+MY (M ) **2 ) 

N — M — 1 

PHIN= PHI(N)*TDRAO 
GA = GAM ( I ) 

GB = (GA+1. )/(2.*(GA-l. )) 

AASlAR = ((2. + (GA-l. )*MACH(M)**2)/(GA+1. ) )**GB/MACH( M ) 
XMN = l.-MX(N)**2 


Figure 40, Program Listing — Generalized Iterative Procedure 
(Continued) , 
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2uvo 
2100 
2110 
2120 
2130C 
2140 
21bU 
2160 
2170 
21b0 
21V0 
2200 
2210 
2220C 
2230 
2240 
22b0 
2260 
2270 
2260 
2290 
2300 
2310 
2320 
2330 
2340 
23b0 
2360 
23/0 
23b0 
2390 
2400 
2410 
2^^0 
2430 
2440 
2450 
2460 
24 10 
2480 
24V0 
?o00 
2510 
2520 
2530 
2o40 
2550 
2 d60 
2d 70 
2580 
2590 


Figure 


XFN = l.+MX(N)**2 
XMM = l.-MX(M)**2 
SINW = SIN(THFN) 

CG5N = casriHFW) 

**** CHFCK FOR UPSTREAM PRUPAGATION 


PHSPD= U(N)+AS(N)*CUSW 
IF( PHSPO.LE.O. ) GU TO 175 

GMN = KNQKM (M )*SINN/( 1 .+MX (N )*CC)SN+MY (N ) *S INN ) 

IN = XMN*SINN 

iD = XPN*CUSN+2.*MX(N) 

i'HBN = ATAN2( TN,TU ) 

TERM = -GMN *MX ( M ) * ( 1 . -(5MN*M Y ( M ) ) 

**** CliECK FOR TOTAL' REFLECTION 

RUCL= ( 1 . -GMN*MY { M ) )**2-XMM*GMNA*2 
IF( RDCL.LE.O, ) GO TO 175 
RAUICL= GMN *SQRT( ROOD 
TN = -TERM+RAUrCL 
TO = (l.-GMN*MY(M))**2-GMN**2 
I'HFM = ATAN2( TN,TU ) 

IF( N.NE.l ) GO TO 70 
fHR ( K ) = T0UEG*THBN 

70 THBM = ATAN2( TN-2.*T’ERM , lO ) 

MACUM = MACiHM) 

MACHN = MACH(N) 

ALFAN = ATAN2( MY(N),MX(N) ) 

BETAM = ATAN2( MY(M , MX(M) ) 

IF( MACH(M) .LT. 1. ) iO TO /I 

00 r = COSCBETA-VO 

BM = BETAM 

BEFAM = ARCOSIAAS'iArt* -JT) 

BE i'AM = S I GfJ ( 1 , BM ) *b ETAH 

71 CONTINUE 
KMYKM = GMN 
KNYKN = KMYSAV 

KMXKM = (1.-GMN*MY(m))/MX(M) 

KNXKN = KMXSAV 
KMYSAV= KMYKM 
KMXSAV= KMXKM 
OKNN = 0. 

IF( N.EQ.l ) GO 10 75 
QKNN = l./SQRT(K^iXKN**2+KNYKN**2) 
lb QKMM = 1./SQRTIKMXKM**2+KMYKM**2 
A(1,1,N)= MX(N)+COS(THFH) 

A C 1 , 2 , N )= MX (N) -COS (THBN ) 

A(1,3,N)= KNYKN*QKNN 
A ( 2 , 1 , N ) = 1 . +' ACUN*COS ( ALF AN-THFN ) 

A ( 2 , 2 , N 3 = 1 . -MACHN*COS ( ALFAN +THBN ) 

A (2 , 3, N ) = QKNN* ( MX ( N ) *KNYKN-MY (N ) *KNXK • ) 

0 ( 1 , 1 , N ) = ( MX (M ) +COS ( I'HFM ) ) / AMQAiI (M ) 


40, Program Listing - Generalized Iterative Procedure 
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2oU0 
2olO 
2620 
2630 
2640 
26b0 
2660 
26/0 
2 6 BO 
2090 
2700 
2 700 
2/10 
2 72 0 
2 130 
2 I40& 
2750C 
2 /60 
2 770 
2 /oOC 
2/90 
2 800 
28 lo 
2840 
2880 
2600 
2870 
2 880 
2890 
2900 
2910 
2920 


D ( 1 , 2 , N ) = ( MX ( M ) -COS ( TH8M ) ) /AMOAN ( M ) 

0( 1 , 3,8)= KMYKM*OKMM/AMOAN(M) 

0 ( 2 , 1 , N ) = HHORAT ( M ) * ( 1 . +MACHM*CUS ( FiETAM- fUFM ) ) 
c) ( 2 , 2 , N ) = RHORAT ( M ) * ( 1 . -MACHM^fCOS ( BETAM+ I HBW ) ) 
o (2 , 3 , N ) = RHOl/AT ( M ) * ( QK MM* ( mX ( M ) *KMYKM-M Y ( M ) *KMXKM ) ) 
A(3,1,N)= 0. 

A(3,2,N)= 0. 

A(3,3,N)= 0. 

J(3,1,N)= SIN(BETAM-THFM) 

0(3,2,N )= -SIN(BETAM+THBM) 

0 ( 3 , 3 , N ) = OKMM* ( KMYK M*S I N ( BETAM ) +KMXKM*COS ( BETAM ) ) 
SMACH(M)=0. 

IF( MACH(M).UT.l. ) GO TO 90 
GO lU 80 

78 FORMAKbX, '***UOWNSTREAM RELATIVE FLU8 AT ROW'', 13, 

IX, 'IS SURERSOMIC***'///) 

***** COMMUTE INVERSE OF MATR-IX AND STORE 

80 CALL DINVER( N ) 

***** COMPUTE MATRIX PRODUCI ***** 

GO lO 15b 

Oij dMACH(M)=1. 

92 MAt)S= SQl7r(MX(N)**2+(MY(N;+VR(W)/AS(N))**2) 

. )= (GAM(N)+1.)*MABS 

t-6= 2. *( 1 . +mABS*^2* ( GAM ( N ) -1 . ) /2. )/MX (N ) 

F7= (GAM(N)-l. )*(MAbS**2-l . ) 

' 1,1, N)= E /+F6*C0S ( I'HFn > -Fb*C(JS ( PHIN-THFN ) 

A C 3 , 2 , ) = F 7-F8*C0S ( THBN ) +F-. *COS ( I’H I N+THBN ) 

A ( 3, 3, N ) = F6*KN YKN-f-b* ( Ko YfoJ*CJS ( PlIIN )-KNXKN*SI 0 ( PH I N ) ) 
J(J,1,N)=0. ; U(3,2,N) = 0. ; L. 3,3,N)=0. 

00 10 6b 


292b 

2927 

2930 

2940 

294'j 

2960 

2968 

2966 

2967 
2970 
2980 
2 990 
3000 

3030 

3031 

3032 
3038 
3036 
303/ 


98 vi = 1 

3AO=u . 

1'. /. = M + 1 

IFi ^'.GT.NSfA ) GO i.J IbO 
Ni ^3-1 

( o'viACH(M) .EO. 1. ) GO 10 110 
1F( MISTAK.EO.O ) B(l'l) = B(N) + 0.01 

80. ii 1 - . ( 1 ) 

B1(K)= BGNl 
GO lu 160 
110 CONTINUE 


BGS = b(N) 

BOM) = -( A(3, 1,N)*F(N) + A( 8, 3, J )*Q( N ) ) /A ( 3, 2 , N ) 
BKK) = B(l) 

IF( N.Gi.l ) GO To 123 
B(;^i + l) = -30. 

IF( EREQ(L) ,EQ. 1600 . BOv+l)= -8b. 

JO iJ 130 


Figure 40 


Program Listing - Generalized Iterative Procedure 
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3038 

3039 

3040 

3041 

3042 
304 8 
3044 
304-j 
304u 
30-4 ' 
304.. 
3049 
3000 
3060& 
30 /O 
30bU(3. 
3090 
3100& 
3110 
312 u 
3180 
3140 
31b0 
31o0 
31b0 
3190 
3200 

3202 

3203 
3206 
3208 
321u 
3220 
32 JJ 
3240 
3280 
3260 
32/0 
3280 
3290 
3300 
3310 
3320C 
3380 


l<i8 JNl = ABS(BGS-B(N)) 

y(N-t-l) = -40. 

1F( M.BQ.NSTA ) IJ(N■^1)= 0. 
IF( ONB.LT.0.05 ) G(J TO 130 
IF( NGS.GT.l ) GO TO 125 
fJlISlAK = 0 


-M = NGS 
GO TO 100 

l2b B(NGS) = BGS-^.01 
M = NGS 
N = NGS - 1 
130 HGS= N + 1 


13b 

140 

150 


A(1,2,N)*B(N) 

A(2,2,N)*I3(N) 


+ A(1,3,N)*Q(N) 
+ A(2»3,N)*Q(N) 


Cl= A(l, 1,N )*F(N) + 

0( 1 , 2 , N ) *B(N+1 ) 

02= A(2, 1,N)*F(N) + 

D(2, 2, N)*B(N+1) 
h'(N + U= (I)(2,3,N)*G1 - U(1,3,N)*C2)/ 

Fc. ruT " <U{2,1,N)*1>(1,3.N))) 

IF( iHI(K).Gi.O. ) GO 10 13b 

IF( thick ).LT. 0. ) 00 TO 13b 

0(iJ + l)= (C2-[)(2, l,N)*F(iM+l ))/U(2,3,N) 

GO TO 140 

0(Nh-1)= (C1-u( 1, l,N)*F(N-t-l))/D(l,3,N) 

CONI INUt 
GO TO 100 
CONTINUE 


IF(B(N+1 ).EQ.O. ) GO TO lb2 
IFCBAG. EQ.O. ) BAG=B(N+1) 
8AT= B(N+1 )/dAG 
IFCdAT.LT.O. ) GO TO 152 
bag = B(N+1) 

GO TO 16b 


152 F,V(K)= F(N + 1) 
^(K) = Q(N+1 ) 
B(K)= Bl(K) 


iHT(K)= THFM*TUUEG 
GO TO 1/0 


’5b CALL MAH 200 ( N ) 
i'HFN = ihi-M 
GO 10 6b 


***** STONE AMFLITUOES 


3400 

3410 

3420 

3430 

3440 

3450 

3460 


160 


16b 


FTN)*PNUD(1,1,N)+B(N)*PW0D(1,2,N)+Q(N)*FR0L 
Bt N+1 ) F(N ) *PP0D(2, 1 ,N) +B(N )*PROD (2, 2,N )-^Q(N )*PROL 

rn'i \ ^ ''' ^ ( 3 , 1 , N ) +B ( N ) *PROD ( 3, 2 , N ) +Q (N ) *PROE 

ou 10 100 

CONTINUE 

IF( NGS.GT.l ) BGS=B(NGS) 


1.3, N) 

2.3, N) 

3.3, N) 


Figure 40. Program Listing - Generalized Iterative Procedure 
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3461 

3462 

3463 
3404 
3466 
346B 
34 IOC 
34B0 
3490 
3b00 
3bl0 
3b2U 
4d3U 
3b4U 
3bb0 
3bo0 
3b VO 
3boO 
3d 9m 
3600 
3610 
3620.1 
3630 
3632 
3634 
3636 
3o40 
36bO 
36b 1 
36b2 
36b4 
3680 
3690 
3/00 
3/10 
3/20 
37bOC 
3 /60 
3790 
3800 
3810 
381b 
3820 
3830 
3840S. 
384b 
38bO 
3860 
3870 
3880 
3890 


IF( NGS.EQ.l ) GO TO 167 
GO iG 12b 

167 B(l) = B(l) + .0001 
M = 2 
N = 1 
GO TO 160 

COMPUTE TRANSMISSION LOSS 

1 /O ASl = AS(1 ) 

A5N = AS(NSTA) 

mX(1) = U(1)/AS1 

MX(NSTA)= U(NSTA)/ASN 

VVl = U(1)*TAN(T0 kAL)*PHI(1)) 

VVN = U(NSTA)*TAN(T0HAD*PHI (NSTA)) 

MY(1) = VVl /ASl 
MY (NSTA )= VVN/ASN 

HHOkA = p(NSTA)*T( 1)*ASN/(P(1)*1(NSTA>*AS1) 
ihPM2 = ( l.+MX( NSTA) *CGS (THEM) + 

(MY (NSTA )+VR(NSTA-l )/ASN)*S IN (THEM) )* (COS (THEM )+MX( NSTA ) ) 
i'HIN = THI(K)*TORAu 
lERMl = (1 .+mX(1)*C0S(THIN)+MY(1)* 

S IN ( THI N ) ) * ( COS ( THIN ) +MX ( 1 ) ) 

: . . -.S(K)= 10.*AL0G10(RH0KA*ABS(TERM1/TERM2)/TW(K)**2) 

1 • iW(K) .LT.O. ) GO TO 177 
IE/ iN(K).GE.l. ) GO TO 177 
IE. ^•JJSS(K) .LT.O. ) TLOSS(K) = 1. 

• >.1 iO 180 
1 /b TLOSS(K)= 20. 

BKK) = 1. 

GO TO 1/8 
177 fLOSS(K) = 2b. 

1/8 TW(K) = 0. 

180 CONTINUE 
SUMI= 1. 

TL = TL0SS(K)/10. 

S'lMT= SUMT+E(K)/10.**TL 

print output 

PRIN. 23b, THI (K),THR(K) ,THT( K ) ,TW.( K) , B1 ( K ) , E( K J , TLOSS (K ) 
23b E0RMAT(F9.3,2E10.3,E9.4,3ElO.4) 

240 CONTINUE 

LSIGMA = 10.*AL0G10(SUMI/SUMT) 

PRINT 24b,FRE0(L) .TLSIUMA 
2 4b FORMAT ( / 1 4X , ' FREQU ENCY= , 1 4 , 1 X , ' HZ' , bX , 

' TR AN SM I SS I Oiv LO SS= ^ , E6 . 2 //// ) 

300 CONTINUE 

GO TO 15 

400 STOP 
END 


Figure 40. Program Listing — Generalized Iterative Procedure 
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3900CU 

AMX 

FUNCTION GAMX(T) 

3910 


FUNCTION GAMX(T) 

3920 


1F( T.LE.800. ) GO i'O 10 

3930 


1F( T.GE.3600. ) GO TO 12 

3940 


GAMX = 2.23708/1*^.070271 

3950 


GO Tu 15 

396>j 

10 

GAMX =1.4 

39/u 


GO TO 15 

3980 

12 

GAMX = 1.254 

3 99u 

15 

RETURN 

4000 


END 


4010 

4020CtJINVE^i CALCULATE INVERSE OF MATRIX 0 

4030 SUBROUTINE OINVERI N ) 

4040 common /CMAfRX/ 0( 3 , 3, IS ) ,0 1 ( 3, 3, 15 ) , A( 3, 3» . RROD( 3, 3, ) 

4050 DIMENSION DU ( V, 15 ) , DUI ( 9, 15 ) , LABEL(3 ) 

4060 EQUIVALENCE ( DU( 1 , 1 ) , D( 1 , 1 , 1 ) ) , ( DDI ( 1., 1 ) , DI ( 1 , 1 , 1 ) ) 

40/0 DIMENSION i'EMR(3,3) 

4090 DIMENSION AA(9,1S) 

4100 EQUIVALENCE ( AA ( 1 , 1 > , A ( 1 , 1 , U ) 

4110 

4120 NN = N 

4130 DO 10 1=1,9 

4140 DDI(I,N)= DD(I,N) 

4150 10 CONTINUE 

4160 CALL MTINV(DDI(1,N),3,3,3,LABEL) 

^170 20 RETURN 

4180 

4190C **** ENTRY MAPRO!) ** COMHJTl PRODUCT OF 01 AND A 

4200 

4210 ENiRY MAPROIK N ) 

4220 NN = il 

4230 

4240 L = NN 

4250 DU 60 J=l,3 

4260 DO 50 1=1,3 

4270 TEMP(I,J)= .. 

4280 DO 40 K=l,3 

4290 TEMP(I,J)= TEMPd ,J)+DI ( I,K,L)*A(K, J,L) 

4300 40 JONIINUH 

4310 50 v^.hTINUE 

4320 60 COniINUE 

4324 DO 100 J=l,3 

4326 DO 100 1=1,3 

4328 100 PROD(I,J,L) = TEMP (I, 3) 

4330 

4370 200 RETURN 

4380 END 

4390 

5060 S r 3P 

5f)70 - , 

Figure 40. Program Listing - Generalized Iterative Procedure 
(Concluded) . 


I 


134 



07/13/77 


09.031 


LOADER DIAGNOSTICS 

<rt> LOADED PREVIOUSLY 

INPUr FILE NAME = DFLP39 


NASA CORE HOT HPTs 3.0 PR 100% N 

1 STAGES 

* AERO-THERMO PARAMETERS * 


STAGE 

STATION 

U- EPS 

PHI- DEG 

VR- EPS 

PS- PS I A 

TS- DEG R 

1 

1 

388.000 

0. 

0. 

54.730 

1375.000 

1 

2 

512.000 

67.300 

819.000 

36.490 

1 261 .000 


3 

857.000 

-38.300 

0. 

1 8.230 

1050.000 



**** FIRST 

CUT -ON 

OCCURS AT 369. HZ **** 



THETA- I 

0. 

THETA-R 

0 . 

THETA- 

0. 

■T T B 

0.2808 0.9600 

E 

1 .0000 

T-LOSS 

4.7659 


FREQUENCY^ SO HZ TRANSMISSION LOSS= ^.77 

0. O. 0.P808 0.9600 . 1.0000 4,7659 

FREOUENCY= 63 HZ TRANSMISSION LO.-35= a.H 

0. 0. 0.2808 0.9600 1.0000 4.7659 

t-REOUENCY= 80 HZ TRANSMISSION LOSS= 4.77 

0. 0. 0.2808 0.9600 1.0000 4.7659 

FREOUENCY= 100 HZ TRANSMISSION L0SS= 4.77 

0. 0. 0.2808 0.9600 1.0000 4.7659 

FREOUENCY= 125 HZ TRANSMISSION L0SS= 4.77 

0. 0. 0.2808 0.9600 1.0000 4.7659 

FRE0UENCY= 160 HZ TRANSMISSION LOSS= 4.77 

0. 0. 0.2808 0.9600 1.0000 4.7659 

Figure 41. Sample Output. 
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0 . 


0 . 


0 . 

76.385 
—76. 385 


0 . 

56.016 
-56.01 6 


0 . 

42.856 

-42.856 


0 . 

33.nyrj 

-33.098 

76.385 

-76.385 


0 . 

26.206 

-26.206 

56.016 

-56.016 


0 . 

20.835 
-20.835 
43.237 
-43.237 
71 .545 


rREQUENCY= 200 HZ TRAMS IISSIOM L0SS= 4.77 


0 . 0 . 

FREQUENCY= 250 

0 . 0 . 

FREQUENCY= 315 

0 . 0 . 

0 . 0 . 

-53.807 0. 

FREQUENCY= 400 


0 . 0 . 

0 . 0 . 

-37.869 -79.743 

HREOUhMCY= 500 


0.2808 0.9600 

HZ TRANSMISSION 

0.2808 0.9600 

HZ TRANS.MISSION 

0.2808 0.9600 

0 . 1.0000 
0 . 1.0000 

HZ TRANSMISSION 


0.2808 0.9600 
O. 1.0000 
0.3443 0.1700 

HZ TRANS.MISSION 


1.0000 4.7659 

LOSS= 4.77 

1.0000 4.7659 

L05S= 4.77 

0.3333 4.7659 

0.33.33 20.0000 

0.3333 20.0000 

L()SS= 9.28 

O. 3333 4. 7659 

0.3333 20.0000 

^.3333 5.8353 

LOSS= 6.96 


28.409 37.114 

-28.409 -58.120 


0 . 

0.4 12^: 


0.7300 

0.1200 


JO j 

0 .. 33.33 
0.3333 


FREQUENCY= 630 HZ TRANSMISSION LOSS= 5.41 


0 . 0 . 

21.698 30.422 

-21.698 -43.183 

0 . 0 . 

-53.807 0. 

rREOJENCY= 800 HZ 

0 . 0 . 

17.076 25.098 

-17.076 -33.188 

0 . 0 . 

-37.869 -79.743 

HREQUEMCY=1 000 HZ 


O. 

I 3.525 
-I 3.525 
23.676 
- 23,616 

0 . 


0 . 

20.593 

-25,1A2 

37.355 

-58.721 

0 . 


0.2808 0.9600 

0.2810 0.7300 

0.5050 -0.3500 

1.0000 

I .0000 


0.2000 

0.2000 

0.2000 

0.2000 

0.2000 


0.2808 0.9600 

n.28ll 0.7800 

-8.8000 
1.0000 
0.3443 0.1700 


0.20O0 

0.2000 

0.2000 

0.2000 

0.2000 


4.7659 

25.0000 

2.7835 


4. 7659 
4.4967 
1 . OOOO 

20.r\mn 
20. 0000 


transmission L0SS= 5.22 


4. 7659 
4.6035 
< ■. tOOO 
20.0000 
5.8353 


O.PBOo 

0.2729 

0 . 

0 . 

0.4092 

0 . 


nssioN 

LOSS= 7.19 


0.9600 

0. 1429 

4.7659 

0.8200 

0.1429 

4.9128 

2 . 0000 

0. 1429 

25.^^000 

0.7200 

0. 1 429 

25.^000 

0.1300 

0. 142'-^ 

2.8780 

1 .0000 

0. 1429 

20.0000 


Figure 41, Sample Output (Continued), 
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- / I 

-49.848 

-1 1 2.2.,i4 

0.3437 

0.0300 

0. 1429 

1 1 . 930 7 


FREQUI-: 

NCY=!250 HZ 

T'v 

NSMISSION 

LOSS= 7.47 


0. 

0. 

0. 

0.2808 

0.9600 

0.1111 

4 . 7659 

16.209 

10.495 

16.455 

0.2766 

0.8500 

0. 1 1 1 1 

4.8281 

-1 6.209 

-10.495 

-I 9.588 

0. 1495 

1 .3600 

0 . 1 n 1 

10.4055 

33.098 

21.698 

30.422 

0.28I-; 

0.7300 

0.1111 

4. 4967 

-33.098 

-21 .698 

-43. 1 83 

0.5850 

-0.3500 

0.1111 

1 .0000 

61 .832 

0. 

0. 

0. 

1 .0000 

0. I I 1 1 

20.0000 

-6 1 . 832 

-34.802 

-72.654 

0.3594 

0 . 1 800 

0. 1 1 1 1 

4.8793 

76.385 

o. 

0. 

0. 

1 . 0000 

n. 1 1 1 1 

20.0000 

-76.385 

-53.807 

0. 


1 .0000 

n. 1 1 ! 1 

20.0000 


h^EQiiENCY=l 600 HZ 

TRANSMISSION 

f)5S= 6.00 


0. 

0. 

n. 

O.2808 

0.9600 

0.0909 

4. 76 3 9 

12.938 

R.365 

1 3.382 

0.2809 

0.8700 

0.0909 

4.7115 

-12.938 

-8.365 

-1 5.383 

0.2073 

1 . 1 900 

0.0909 

7.5142 

26.206 

1 7.076 

25.098 

r'.2H : 1 

n . 7800 

n.0909 

4.6035 

-26.206 

-17.076 

-33.1 8,8 

0. 

-8. 8000 

' .0909 

25.0000 

40.266 

26.606 

35.433 

o. 

1 .0200 

0.0909 

25.0000 

-40.266 

-26.606 

-54 . 069 

' . .374 

0.0700 

0.0909 

2.0704 

56.016 

0. 

0. 

0. 

1 .0000 

0.0909 

20.0000 

— 8' . -)1 6 

-37.869 

-79. 743 

0.3443 

0. 1 700 

0.0909 

5.8353 

76.385 

0. 

0. 

0. 

1 . oooo 

0.0909 

20.0000 

-76.385 

-53.307 

0 . 

0. 

1 .0000 

0.0909 

20.0000 


FHEQUENCY=200o 

IRANSMISSIcm 

LOSS= 7.16 


0. 

0. 

0. 

0.2808 

0.9600 

0.0769 

4.7659 

1 0. 335 

6.676 

1 0.849 

^.2152 

0.8900 

0.0769 

4.8978 

—10. 335 

-6.676 

-12. i.;8 

0. ^271 

1.1200 

0.0769 

6.6908 

20.835 

’ 3.526 

20. 59'^ 

0.2/29 

0.8200 

0.0769 

4.91 28 

-20.835 

- i 3.525 

-25.742 

0. 

2.0000 

0.0769 

25.0000 

'oo 

20.754 

29.384 

0.2822 

0.7400 

0.0769 

4.4865 

-ji . /on 

-20. 754 

-41.118 

0.6498 

-0.5500 

0.0769 

1 .0000 

4 3.28 7 

23.676 

37.355 

0. 

0.7200 

0.0769 

25.0000 

-4 3.23 / 

-23.676 

-58.721 

0.4092 

0. 1 300 

0.0769 

2.8786 

56.01 


0. 

0. 

1 .0000 

0.0769 

20.0000 

-56.01 6 

-J / .B69 

-79.743 

0.3443 

0. 1 700 

0.0769 

5.8353 

71 .548 

0. 

0. 

0. 

1 .0000 

0.0769 

20.0000 

-71 .545 

-49.343 

-112.284 

0. 3437 

0.0300 

0.0769 

1 1.9307 


t-RE0UENCY=25O0 3/ 

TRANSMISSION 

L0SS= 6.07 


0. 

0. 

0. 

0.2308 

0.9600 

0.0588 

4.7659 

8.195 

‘ .291 

8.706 

0.2840 

0.9000 

0.0588 

4.6343 

-8. 1 95 

-t>.29l 

-9.51 I 

0.2454 

1 .0700 

0.0588 

5.9928 

16.470 

1 0. 666 

16.695 

0.2721 

0.8500 

0.0588 

4.9678 

-16.470 

- 1 0. 666 

-19.928 

0.1427 

1 .3800 

0.0588 

10,8168 

24.91 3 

16.220 

24.046 

0.2780 

0.7900 

0.0588 

4.7129 

-24.918 

-16.220 

-31 .373 

0. 

12.5800 

0.0588 

25.0000 

33.655 

22.076 

30.331 

0.3049 

0.7200 

0.0588 

3.7790 

-33.655 

-22.076 

-44.012 

0.5640 

-0.2800 

0.0588 

1 .0000 

42.856 

28.409 

37. 1 1 4 

0. 

0.7300 

0.0588 

25.0000 

-42.856 

-23.409 

-58. 1 2n 

0.4122 

0. 1 200 

0.0588 

2.7835 

52.809 

0. 

V 

• 

n 

1 .0000 

0.0588 

20.0000 


Figure 41. Sample Output (Continued) . 


-52.809 

-35.51 3 

-74.288 

64. 109 

0. 

0. 

-64. 109 

-43.985 

-94.481 

78.514 

0. 

0. 

-78.514 

-55.586 

0. 


FREQUENCY=3150 1 

0. 

0. 

0. 

6.449 

4. 162 

6.919 

-6.449 

-4.162 

-7.4 1 8 

12.938 

8.365 

1 3.382 

-12.938 

-8.365 

-15.383 

19.507 

12.653 

1 9.430 

-19.507 

-12.653 

-23.951 

26.206 

17.076 

25.098 

-26.206 

-1 7.076 

-33. 1 88 

33.098 

2 1 .698 

30.422 

-33.098 

-2 1 . 698 

-43. 1 83 

40.266 

26 . 606 

35.433 

-40.266 

-26.606 

-54.069 

47.837 

0. 

0. 

-47.837 

-31 .927 

-66.084 

56.01 6 

0. 

0. 

-56.016 

-37.869 

-79.743 

65.200 

0. 

0. 

-65.200 

-44.830 

-96.655 

76.385 

0. 

0. 

-76.385 

-53.807 

0. 

95.962 

0. 

0. 

-95.962 

-71.192 

0. 


HRE0UE:NCY=4000 HZ 


0.3556 

0. 1 800 

0.0588 

5 

0. 

1 .0000 

0.0588 

20 

0.3290 

0.1200 

0.0588 

7 

0. 

1 .0000 

0.0588 

20 

0. 1.0000 
TRANSMISSION 

0.0583 
LOSS= 6.37 

20 

0.2808 

0.9600 

0.0435 

4 

0.2871 

0.9100 

0.0435 

4 

0.2551 

1 .0400 

0.0435 

5 

0.2809 

0.8700 

0.0435 

4. 

0.2073 

1 . 1 900 

0.0435 

7. 

0.2704 

0.8300 

0.0435 

5. 

0.0402 

1.7100 

0.0435 

21 . 

0.281 1 

0.7800 

0.0435 

4. 

0. 

-8.8000 

0.0435 

25. 

0.2810 

0.7300 

0.0435 

4. 

0.5850 

-0.3500 

0.0435 

1 . 

0. 

1 .0200 

0.0435 

25. 

0.4374 

0.0700 

0.0435 

2. 

0. 

1 .0000 

0.0435 

20. 

0.378 

0. 1 700 

0.0435 

3. 

0. 

1 .0000 

0.0435 

20. 

0.3443 

0. 1 700 

0.0435 

5. 

r- 

1 .0000 

0.0435 

20. 

0.328b 

0.1100 

0.0435 

8. 

0. 

1 .oooo 

0.0435 

20. 

0. 

1 . 0000 

0.0435 

20. 

0. 

1 . 0000 

0.0435 

20. 

0. 1.0000 
TRANSMISSION 

0.0435 
L()SS= 6.31 

20. 


Figure 41, Sample Output (Concluded). 
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.0964 

.0000 

.9326 

.0000 

.0000 


.7659 
.5451 
,6418 
,7115 
,5142 
,0045 
, 8024 
6035 
0000 
4967 
0000 
0000 
0704 
0000 
9821 
0000 
8353 
0000 
2839 
0000 
0000 
0000 
0000 
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APPENDIX C 


COUPLING OF LINE SOURCE TO DUCT MODES 


To understand the coupling of the sound source with the various duct 
modes possible in an annulus, we consider an idealized problem in which the 
annulus is unwrapped into a rectangular duct and the siren source is modeled as 
a line source (see Figure 42); "r** is the mean radius of the annulus. 

We have to solve an equation for the pressure (denoted by : 

V^p + k^p = 6(x)6(y) 

where k = o)/a, <5(x), <S(y) are delta functions. The solution is for - ttt _< y 
^ TTr and has to be periodic with wavelength ”27ir*\ Also, at x ±<», the 
radiation condition is to be satisfied. Let 


p = A (x) cos (^) , 

m r 

m=0 


Note that 6(y) can be expanded in the even Fourier series as: 

oo 

{- 2 ~ + ^ ^ cos (^)). Then A[J (x) + Aq(x) = 6(x)/27rr and 


A” (x) + [k^ ~ (~)^] A (x) = 6(x)/7rr for m > 1. The solution for p can be 
m r m — 

finally written down as; 


gj k |x| ^ 

47Tj k r ^ ^ 


N 


P = 


i ‘‘m l^'lcos (SE) 


2j k nr 
m 


N+1 


e -"m 1^1 cos (^) 

2k TTr 
m 


= - (7) 


where k =Wk‘‘ - (~) and N is the largest value of m for which kr > m and 
m * “ 


V m 2 2 

(— ) - k for m > N. We are not interested in the terms of the 

r 

series involving in the above since they represented nonpropagating terms. 
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The energy flux associated with each of the propagating terms can be deduced 
by first writing 


cos (^) 
r 


(niy/r)_^^-j (my/r) 


and considering each cos (my/r) term to involve two plane waves (one for +y 
and other for -y and then noting that the energy flux will be proportional to 
cos (niy/r) x (km/k). The square of the cosine term is proportional to the 
power in the wave direction and it is the axial power component that is of 
interest. Hence the product of the cosine squared and the direction cosine is 
considered. Since the cut-off frequency for each mode is = (am/r) , the net 
result is that the line source will excite energy levels in each propagating 
mode proportionaj. to {f - where f is the frequency of excitation and 

the cut-off frequency of the relevant mode. 
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NOMENCLATURE 


a 

A/A* 

A. . 
ij 

®( ) 
BPF 

D. . 


m 


( ) 


H 

-> 

I 


J 

k 

K 

£ 

L 

m 

m 

M 

N 

P( ) 

^ ) 
r 

R 

o 

S 

t 

T 

( ) 
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Speed of sound, m/sec 
Isentropic area ratio 
Matrix element 

Upstream coefficient matrix for n-th blade row 
Amplitude of backward-travelling wave 
Blade passing frequency, Hz 
Matrix elements 

Downstream coefficient matrix for n-th blade row 
Unit vector 

Energy assignment to m-th mode 
Frequency, Hz 

Amplitude of forward-travelling wave 
Aeroacoustic flow function 
Height, cm 
Intensity vector 


Wavenumber, co/a 
(k^ + 

X y 

Axial spacing, cm 
Length, cm 

Circumferential lobe number 
Mass flow rate, kg/sec 
Mach number 
Newtons 

2 

Pressure, N/m 

Amplitude of vorticity wave 
Radius , cm 


Real component of impedance at an interface 

Cross-sectional area 

Time coordinate 

Amplitude of transmitted wave 

Temperature, K 
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NOMENCLATURE (Continued) 


Transmission loss, dB 

Perturbation in axial velocity component, m/sec 
Mean axial velocity component, m/sec 
Perturbation in transverse velocity component, m/sec 
Mean transverse velocity component, m/sec 
Rotor physical speed, m/sec 


Axial Cartesian coordinate, fixed to blade row 
Reactive component of impedani e at an interface 
Transverse Cartesian coordinate, fixed to blade row 
Upstream relative flow angle, degrees 
Downstream relative flow angle, degrees 
Ratio of specific heats 
Increment or decrement 
Acoustic energy density 


Strouhal number (dimensionless frequency) 


Wave propagation angle 


(1 + ^ 

2 abs 

Dimensionless 


Dimensionless 


) 

velocity 
length parameter 


3.14159 

3 

Density, kg/m 
Hub/tip (radius) ratio 


Absolute flow angle relative to axiai direction, degrees; also, 
dimensionless pressure in Section 3.2, <j> = p-/yp 

Angle of Inclination of acoustic wave incident on a shock 
Angle of inclination of refracted wave leaving shock 
Angle of inclination of vorticity wave behind shock 


Circular frequency, 2iTf 



NOMENCLATURE (Concluded) 


Subscripts 

a.bs Absolute 


flow parajeter 


( ) 


c 

f 

F 


( ) 


1 

I 


m 

n 


Q 

R 


( 

( 


s 


) 

) 


Backward- travelling wave parameter 
Cut-on 

Final value in passage problem 
Coordinate frame moving with fluid 

Forward-travelling wave parameter 
Initial value in passage problem 
Incident wave parameter 
Value downstream of n-th blade row 
Value upstream of n-th blade row 
Vorticity wave parameter 

Reflected wave parameter 
Static value 


T Transmitted wave parameter 

X Axial component 

y Transverse component 

0 Stagnation value 

1 Station upstream of turbine 


Superscripts 

( ) * Perturbation quantity 

( )* Root mean square value 
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